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ABSTRACT 

A  general  approach  to  the  elastic  and  inelastic  analysis  of 
thin-walled  members  of  arbitrary  open  cross  section  is  presented  and 
a  finite  element  formulation  is  developed. 

Using  this  formulation,  the  solution  for  critical  loading  as 
well  as  nonlinear  load  deflection  response  for  elastic  members  is 
obtained. 

Numerical  results  are  presented  for  simple  cases  when  local 
and  member  buckling  are  coupled. 

Inelastic  beam-column  problems  are  solved  using  an  iterative 
incremental  technique  based  on  an  equilibrium  balance,  and  inelastic 
buckling  problems  are  solved  using  a  standard  eigenvalue  solution  to 
evaluate  the  critical  length.  Numerical  examples  are  presented  and 
compared  with  available  results. 
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CHAPTER  I 


INTRODUCTION 

1. 1  Introductory  Remarks 

Thin-walled  beams  occupy  an  important  place  in  the  field  of 
modern  structural  engineering.  Rolled,  welded  or  rivetted  metallic 
beams,  columns,  girders  and  elements  of  frames,  are  examples  of  thin- 
walled  members.  In  addition,  many  spatial  engineering  structures, 
such  as  shear  wall  structures,  have  such  proportions  that  they  can  be 
considered  as  thin-walled  structures. 

The  use  of  cold-formed  structural  members  is  increasing,  and 
hence  the  range  of  loading  for  which  elastic  analysis  is  applicable  is 
also  increasing.  The  stability  of  such  structures  has  been  the  subject 
of  many  investigations  and  the  study  of  stability  analysis  is  an 
important  aspect  of  design  since  it  predicts  an  upperbound  for  the 
ultimate  carrying  capacity  of  the  member.  The  general  approach  to 
such  problems  has  been  to  formulate  a  set  of  equations  for  the  specific 
problem  under  consideration.  In  some  cases  the  design  criteria  is 
based  on  deflection  limitations  and  hence  this  aspect  of  analysis 
should  not  be  overlooked. 

Thin-walled  beam  theory  is  based  on  the  assumption  that  the 
cross  section  does  not  distort.  However,  in  some  cases,  failure  can 
be  initiated  by  local  buckling  which  involves  distortion  of  the  cross 
section.  Hence  it  is  necessary  to  design  such  members  so  that,  at 
design  load,  adequate  safety  exists  against  failure  by  local  buckling. 

The  problem  becomes  more  complex  when  nonlinear  material  or 
geometry  is  considered.  A  three  dimensional  structure  is  often 
analyzed  as  a  collection  of  two  dimensional  planar  structures.  While 
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this  idealization  has  usually  resulted  in  safe  design,  it  does  not 
necessarily  represent  the  best  model  for  the  analysis  of  such  struc¬ 
tures.  To  the  author’s  knowledge  no  general  study  of  the  inelastic 
behaviour  of  beam-columns  is  available.  Today,  many  structures  are 
designed  by  plastic  design  methods  and  each  member  must  undergo  con¬ 
siderable  inelastic  deformation  in  order  for  the  whole  frame  to 
develop  its  full  strength.  These  inelastic  deformations  may  cause 
premature  failure  by  various  types  of  instability.  Hence  there  is  a 
need  to  study  inelastic  instability. 

Although  the  general  differential  equations  for  thin-walled 
beams  (68,74)  are  valid  for  complex  geometry,  loading  and  boundary 
conditions ,  the  complexity  of  these  equations  is  such  that  only  a  few 
problems  have  been  solved  using  classical  mathematics.  More  general 
approaches  consist  of  approximating  the  differential  equations  using 
the  finite  difference  technique,  (26) ,  or  solving  the  equations  by 
direct  numerical  integration  (68) .  Even  these  methods  offer  only  a 
partial  improvement  over  classical  methods  of  solution.  On  the  other 
hand,  the  development  of  a  general  treatment  by  finite  element  method 
is  encouraging  and  has  resulted  in  a  number  of  successful  applications 
(7,8). 

1 . 2  Purpose  and  Scope 

The  objectives  of  this  investigation  are: 

1)  to  develop  a  general  set  of  beam  equations  which  are  applicable  to 
both  elastic  and  inelastic  analysis  and  adaptable  to  finite  element 
analysis  ; 

2)  to  develop  general  computational  procedures  for  the  finite  element 
analysis  of  elastic  and  inelastic  members,  which  are  applicable  to 
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stability  problems;  and, 

3)  to  investigate  the  effect  of  distortion  of  the  cross-section  on  the 
member  buckling  strength  for  some  simple  cases. 

1. 3  Outline  of  Contents 

Beam  equilibrium  equations,  in  both  total  and  incremental 
form,  and  without  any  assumptions  relating  stress  resultants  to  deform¬ 
ations,  are  derived  in  two  ways  in  Chapter  2.  By  specializing  to 
elastic  response  the  equations  are  shown  to  be  equivalent  to  Vlasov's 
equations  (74).  A  finite  element  model  is  then  developed  for  the 
elastic  case.  Details  of  derivations  of  equations  are  contained  in 
Appendix  A.  A  comparison  of  nonlinear  terms  in  the  potential  energy 
and  virtual  work  formulations  is  carried  out  in  Appendix  B  where  the 
range  of  applicability  of  the  derived  equations  is  investigated. 

In  Chapter  3,  solution  techniques  to  solve  small-deflection, 
beam-column,  and  bifurcation  problems,  using  the  elastic  finite 
element  equations  of  Chapter  2,  are  discussed.  Numerical  results  for 
a  variety  of  problems  are  compared  with  those  obtained  from  classical 
solutions . 

Chapter  4  introduces  distortion  of  the  cross-section  and 
formulates  the  equilibrium  equations,  in  general  terms,  in  which,  local 
and  member  buckling  effects  are  coupled.  A  finite  element  approxima¬ 
tion  to  these  equations  is  derived  for  the  special  case  of  elastic 
wide  flange  sections.  Numerical  results  for  a  variety  of  local 
buckling  problems  in  wide-flange  sections  are  compared  with  available 
solutions  in  the  literature.  The  coupling  between  local  and  member 
buckling  is  also  investigated.  Details  of  the  derivation  of  the 
matrices  required  in  this  type  of  analysis  are  provided  in  Appendix  C. 
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The  extension  of  the  formulation  of  the  equation  of  Chapter 
2  to  the  case  of  inelastic  response,  is  carried  out  in  Chapter  5. 

The  finite  element  equations  are  formulated  with  respect  to  two  dif¬ 
ferent  systems  of  local  reference  axes,  and  assembly  of  the  'tangent' 
stiffness  and  geometric  stiffness  matrices  with  respect  to  reference 
axes  through  the  original  centroid  and  shear  center  of  the  elastic 
section  is  discussed.  Detailed  derivations  of  the  relevant  section 
properties,  stress  resultants,  submatrices  and  transformation  matrices 
are  contained  in  Appendices  D  and  E. 

Chapter  6  illustrates  the  application  of  the  equations 
derived  in  Chapter  5,  to  the  solution  of  inelastic  small-deflection, 
beam-column  and  bifurcation  problems.  A  short  summary  and  discussion 
is  contained  in  Chapter  7. 
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CHAPTER  II 


FORMULATION  OF  ELASTIC  BEAM  EQUATIONS 
2 . 1  Introduction 

This  chapter  treats  the  formulation  of  the  governing  equa¬ 
tions  for  a  thin-walled  member  of  open  cross-section  by  two  approaches. 

A  set  of  differential  equations,  identical  with  those  of  Vlasov  (74), 
are  derived  from  equilibrium  considerations.  It  is  then  shown  that 
the  same  set  of  equations  may  be  derived  from  virtual  work.  A  finite 
element  formulation,  for  elastic  members,  is  derived  from  the  virtual 
work  equations. 

Original  studies  of  the  flexural  buckling  of  compressed  bars 
were  made  by  Euler  in  the  eighteenth  century.  Derivations  of  funda¬ 
mental  differential  equations  are  associated  with  the  names  of  Chwalla 
(9),  Goodier  (28),  Kappus  (9),  Timoshenko  (68),  Vlasov  (74),  and  others. 
Djanelidze  (18)  derived  the  governing  equations  for  thin-walled  sections 
by  the  energy  method.  The  early  studies  of  torsional  and  flexural 
buckling  have  been  described  by  Bleich  (13). 


2 . 2  Governing  Equations  Derived  From  Equilibrium 

A  thin-walled  open  cross  section  is  shown  in  Fig.  2.1.  Let 

point  ’  C'  be  an  arbitrary  origin  of  coordinates  and  point  ’S',  whose 

coordinates  are  (e  ,e  ) ,  be  an  arbitrary  reference  point.  Let  the 

x  y 

transverse  load  and  support  reactions  pass  through  the  1 S-axis1,  as 


shown  in  Fig.  2.2. 
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Define  the  following  stress  resultants 


a  d  A 
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(2-la) 
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a  y  d  A 
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(2-lb) 
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a  x  d  A 
z 


(2-lc) 


a  d  A 
zx 


(2-ld) 


V 


a  d  A 
zy 


(2-le) 


M 


(a  (x-e)-a  (y  -  e  ) }  d  A 
zy  x  zx  y 


(2-lf) 


where  a  ,  a  ,  a  are  components  of  the  stress  tensor  acting  on  the 
z  zx  zy 

cross  section,  and  A  is  the  area  of  the  section. 

If  the  effect  of  displacements  on  the  equilibrium  equations 
is  neglected  the  normal  small  deflection  equilibrium  equations  result 
from  a  consideration  of  the  equilibrium  of  an  infinitesimal  length  of 
the  beam  (26),  (Fig.  2.2),  and  are 


M"  -Ip  =  0 

x  y 


(2~2a) 


M"  +  p  =  0 

y  rx 


(2-2b) 


M  ’  +  in 
t  t 


0 


(2-2c) 
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where  (T)  indicates  differentiation  with  respect  to  z;  is  the  torque 

about  the  S  axis;  and  p  ,  p  and  m  are  distributed  loads  and  torques 
applied  to  the  beam  (Fig.  2.2).  These  equations  represent  the  sum  of 
the  forces  in  the  x  and  y  directions  and  the  moment  about  the  S  axis, 
respectively,  for  an  infinitesimal  length  of  beam.  However,  if  the 
effect  of  displacements  on  the  equilibrium  equations  is  to  be  con¬ 
sidered,  Equations  2-2  are  inadequate  and  it  is  necessary  to  include 
the  effects  of  the  change  in  orientation  and  location  of  the  stresses 
acting  on  the  cross-section  in  the  summations. 

Let  u  ,  v  denote  the  displacements  of  the  axis  passing 
s  s 

through  S  in  the  x  and  y  directions,  respectively;  denote  the 
(uniform)  axial  displacement  of  point  C  (see  Fig.  2.1);  and  9  denote 
the  angle  of  twist. 

It  is  now  necessary  to  make  some  of  the  assumptions  normally 
associated  with  the  behaviour  of  thin-walled  beams,  (74),  namely 

1)  Because  of  the  bending  and  torsional  flexibility  of  thin-walled 
open  sections,  the  relative  effect  of  shearing  strain  along  the  mid¬ 
surface  of  the  plate  is  extremely  small  and  can  be  neglected. 

2)  The  shape  of  the  projected  cross-section  is  unaltered  during 
deformation. 

Consider  now  a  longitudinal  fibre  of  length  Az  and  area  AA 

with  its  center  at  an  arbitrary  point  A  on  the  midsurface  contour 

(Fig.  2.1).  On  the  application  of  loads,  the  fibre  displaces  uA>  and 

v.  as  shown  in  Fig.  2.3a.  When  it  reaches  an  equilibrium  configuration, 

it  has  developed  a  force  of  a^AA  which  acts  along  the  deformed  axis  of 

the  fibre.  Let  p  and  p  denote  the  total  change  in  transverse  force 

x  y 

per  unit  of  length  resulting  from  the  projection  of  the  normal  and 
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shearing  stresses  due  to  the  deformation.  Assuming  the  section  is 
thin  so  that  shearing  stresses  act  in  the  direction  of  the  contour, 
these  changes  may  he  evaluated  from  Figs.  2.3b  and  2.3c,  which  result 
in  the  rates  of  change  shown  in  Fig.  2.3d.  Summing  over  the  cross- 
section  yields 


p  Az 

y 


a  3v 

3  z  AN  .  . 

~  (a  - — )  AA  - 
az  z  3 z 


— —  (a  d>)  AA  sin  a  }  dz 
3z  zn  ' 


(2-3a) 


p  Az 

rx 


i  3u 

f-  (a  ■—)  AA  - 

o  Z  Z  a  Z 


—  (a  d>)  AA  cos  a  >  dz 

3Z  ZT\  ' 


(2— 3b ) 


The  change  in  torsional  moment  about  S  resulting  from  these  deforma¬ 
tions,  may  be  evaluated  from  Figs.  2.3d  and  2.4,  as 


m^Az  = 


a, 

3z  °z  3z  3z 


(a  <j>)  sin  a\  (x 


zp 


ex} 


AA 


3u 


3z  ^°z  3 z 


A)  „  3 


3z 


(a  <j>)  cos  a  >  (y  -  e  )  AA 


zn 


+ 


3_ 
3  z 


Az 


(2-3c) 


Expressing  the  previously  stated  assumption  for  the  direction  of  the 
shearing  stress,  as 
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-  o  cos  a 

zy  zn 

(2-4a) 

0  „  =  -a  sin  a 

zx  zn 

(2-4b) 

,  expressing  u^  and  by  the  kinematic  relationships  (Fig. 

2.1) 

v  =  v  +  (x  -  e  )d> 

A  s  x 

(2-5a) 

u.  =  u  -(y-e)d) 

A  s  y  T 

(2-5b) 

;  and  carrying  out  the  integration  in  Equations  2-3,  yields 


Px  =  (P  up  +  (P  ey  <*.’)'  -  (Mx  *)" 

(2-6a) 

p  =  (P  v')  -  (P  e  4.  ’  +  (M  «" 

y  o  a  y 

(2-6b) 

—  t  »  ?  ? 

m  =  (M  v’)  -(Pe  v')  -(M  u’)  -f-  (P  e  u') 

t  y  s  x  s  x  s  y  s 

»  »  » 

+  (M  cj>T)  -  (V  v  )  +  (V  u  ) 

p  x  s  y  s 

(2-6c) 

where  it  has  been  necessary  to  define  the  additional  stress  resultant 


» 

M  =  a  { (x  -  e  )  ^  +  (y  ~  e  )  ^}  AA 

p  z  x  y 

A 

(2-7) 

Augmenting  Equations  2-2  by  the  effective  changes  in  distributed  forces 
expressed  in  Equations  2-6,  yields  the  equilibrium  equations  including 


the  effects  of  member  displacements,  as 
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M"  -1  (P  v')'  +  (M  <j))"  -  (P  e  <j>')’  +  p  =  0 

x  s  y  x  y 


(2-8a) 


»*  » 


M"  +  (P  u')  -  (M  <j>)  +  (P  e  4>')  +  p  =  0 

y  s  x  y  x 


(2-8b) 


+  (M  d) ' )  -  (V  v  )  +  (V  u  )  +  m 

P  x  s  y  s  t 


0 


(2~8c) 


It  should  be  noted  that  Equations  2-8  are  valid  for  arbitrary  location 
of  the  reference  points  C  and  S  and  that  it  has  been  unnecessary  to 
specify  the  distribution  of  a  and  a 

Z  ZT) 

2 . 3  Governing  Equations  Derived  By  Virtual  Work 

In  this  section  incremental  equilibrium  equations  and  an 
incremental  stiffness  formulation  are  developed  using  the  principle  of 
virtual  work.  The  principle  of  virtual  work  is  selected  since  it  is 
less  restrictive  than  a  potential  energy  formulation  and  the  total 
and  incremental  equilibrium  equations  developed  by  this  principle  are 
valid  for  both  elastic  and  inelastic  members.  Rectangular  cartesian 
coordinates  and  the  principle  of  virtual  displacements  are  used 
throughout. 

Referring  to  Fig.  2.5,  let  E  represent  the  deformed  equili¬ 
brium  configuration  (the  initial  configuration)  of  the  structure  under 
the  surface  loads  t^,  and  r  -f  AT  represent  the  deformed  equilibrium 

configuration  under  the  loads  T. .  Let  cr. e. and  u.  be  Kirchhoff's 
&  i  ij  ij  i 

tensor,  Green's  strain  tensor,  and  the  displacements  respectively  in 
configuration  T;  and  ,  E_  and  be  Kirchhoff's  stress  tensor, 
Green's  strain  tensor  and  the  displacements  in  configuration  f  +  AF . 
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Green  s  strain  tensor  in  the  equilibrium  configurations  F  and  T  +  AT 
are  given  by. 


2 


e.  . 
ij 


u.  . 


+  u.  . 
3 


(2-9a) 


2  E.  . 
iJ 


U.  . 


+  u.  .  +  u.  .  u. 

J,1  k,r  k,  j 


Let 


(2-9b) 


S . .  =  a. .  +  a. . 

iJ  ij  ij 


E.  .  =  e. .  +  e.  . 

±3  ij  ij 


U.  =  u.  +  u. 
r  li 


(2-10a) 

(2-10b) 

(2-- 10c) 


T.  =  t.  +  t.  (2-10d) 

ill 

where  o...  e..,  u.,  and  t.  are  increments  in  the  corresponding  quan- 

1J  1J  1  1  r  to  n 

tities  going  from  configuration  F  to  F  i  AF. 

Applying  the  principle  of  virtual  displacements  in  position 
T  +  AT ,  and  regarding  the  displacement  increments  u_^  as  variables, 

yields 


S  .  .  6 (E .  . )  dV  = 
ij 


T.  SU.  dS  + 

l  l 


F.  6U.  dV  (2-11) 

l  l 


Substituting  for  S..  and  E..  Equation  2-11  becomes, 
6  ij  iJ 
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ij 
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u , 


J  >i 
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+  "k.i 


\,j  +  \,j 


4  \,i)  dv 


(t.  +  t.)S  u.  dS 

ill 


+ 


F.  6 
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u.  dV 

l 


(2-12) 


The  equilibrium  equation  in  configuration  T;  which  may  be  obtained  by 
letting  cv.  =  0,  t^  =  0,  and  u^  =  0,  is 


n .  .  ( 5  u .  .  +  5  u .  .  +  u.  .  6  ii  .  +  u.  .  6  u.  . )  dV 

ij  i,j  j,i  k,i  dc,j  k,j  hc,i 
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The  incremental  equilibrium  equation  can  be  obtained,  using 
Biot's  (10)  procedure,  by  finding  the  difference  between  the  equili¬ 
brium  equations  for  positions  T  3-  AT  and  T.  Subtracting  Equation  2-13 
from  Equation  2-12  the  incremental  equation  for  position  Y  +  AT  is 
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The  first  term  in  Equation  2-14  represents  the  work  of  the 
"initial  stress"  state  during  the  virtual  variations  of  increments  in 
displacement  gradients  and  is  of  the  same  order  of  magnitude  as  the 
second  term.  This  gives  rise  to  the  "initial  stress"  or  "geometric  1 
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stiffness.  The  second  term  represents  the  work  of  the  increments  in 
the  stress  state  during  virtual  variations  of  increments  in  strain. 
This  gives  rise  to  the  normal  "flexural  stiffness".  The  third  term 
arises  because  of  the  presence  of  the  initial  displacement  gradient 
terms  in  the  virtual  variation  of  increments  in  strain.  This  gives 
rise  to  Marcal's  (44)  so-called  "initial  displacement"  matrix  (see 
Appendix  B) .  The  term  on  the  right  hand  side  represents  the  work 
done  by  the  increments  in  the  external  loading  on  the  virtual  varia¬ 
tion  of  displacement  increments. 

Equation  2-11  represents  the  final  total  equilibrium  equa¬ 
tions.  In  the  following  development  Equation  2-11  will  be  written, 
for  simplicity  of  notation,  by  replacing  with  u_^,  with  a^. 

and  T.  with  t..  Equation  2-11  then  becomes 
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However,  the  reader  must  now  clearly  distinguish  between  the  symbology 
of  Equation  2-15,  which  is  in  terms  of  total  quantities,  and  Equation 
2-14,  which  is  in  terms  of  incremental  quantities.  Equation  2-15  will 
form  the  basis  of  the  remaining  derivations  in  this  section. 

Using  the  symmetry  of  the  stress  tensor,  Equation  2-15  can 

be  written  as 
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where  body  forces  have  been  neglected. 
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This  equation  may  be  written  symbolically  as 


s(Ii  +  q  -  i  )  -  o 


(2-17) 


where 
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For  thin-walled  beams,  the  stress  tensor  is  of  the  form. 
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(2-19) 


where  a  ,  a  are  related  to  a  ,  by  Equations  2-4  (see  Fig.  2.1). 
zx  zy  zq  & 

Substituting  Equation  2-19  for  a^ ^ ,  and  Equation  2-5  for  u^ 
into  Equation  2-17;  and,  applying  the  variation  and  integrating  by 
parts;  yields  the  Euler-Lagrange  equilibrium  equations  and  the  natural 
and  geometric  boundary  conditions.  Details  of  this  procedure  are 
given  in  Appendix  A. 

The  equilibrium  equations  are, 
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where  is  the  distributed  axial  load  per  unit  of  length;  m  and  m 

X 

are  distributed  couples  per  unit  of  length;  W  is  defined  as  the 


bimoment  (see  Appendix  A);  T  is  defined  as  St.  Venant’s  torque  (see 
Appendix  A) ;  and  the  other  stress  resultants  have  been  defined  in 
Section  2-2. 


It  should  be  noted  that  the  St.  Venant  torque,  T  ,  does  not 


sv 


arise  entirely  naturally  from  Equation  2-18a,  but  the  virtual  work  from 
this  effect  has  been  added  as  described  in  Appendix  A. 


Equations  2-20  derived  from  the  virtual  work  principle,  are 


the  same  as  Equations  2-8,  derived  from  the  equilibrium  approach, 
except  for  minor  details  in  notation.  They  apply  to  a  thin-walled 
member  regardless  of  the  constitutive  relationships  and  with  respect 
to  any  pair  of  reference  axes  C  and  S.  All  the  terms  in  which  displace¬ 
ments  appear  are  a  result  of  second  order  geometric  effects.  Discard¬ 
ing  these  terms  yields  the  small  deflection  beam  equations.  Equations 


2-2. 


2 . 4  Displacemen t  E quations  of  Equilibrium 


The  solution  of  Equations  2-20  for  the  general  statically 
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indeterminate  case  is  a  complex  problem.  Since  stress  resultants  are 
themselves  dependent  upon  the  deformations,  it  is  necessary  to  express 
them  in  terms  of  the  dependent  displacement  variables  before  a  solu¬ 
tion  for  displacements  can  be  attempted.  If  "large"  displacements 
were  included  the  expressions  for  the  stress  resultants  would  them¬ 
selves  include  higher  order  terms.  The  introduction  of  these  higher 
order  terms  adds  to  the  complexity  of  the  problem.  Since,  for  the 
usual  type  of  instability  problem  they  have  negligible  effect  on  beam 
behaviour  (see  Appendix  B) ,  it  is  common  practice  to  compute  the  stress 
resultants  using  only  first  order  quantities. 

Because  a  solution  for  displacements  is  required,  it  is 
always  necessary  to  express  at  least  the  first  stress  resultants  in 
each  of  Equations  2-20,  and  T  in  Equation  2~20d,  in  terms  of  the 
displacement  variables.  The  stress  resultants  in  the  remaining  terms 
may  be  expressed  either  in  terms  of  displacements  or  evaluated 
directly  from  statics.  The  selection  of  the  technique  to  be  used  may 
depend  on  the  type  of  problem  under  investigation  but  a  more  efficient 
numerical  solution  is  obtained  when  these  stress  resultants  are 
evaluated  from  statics,  if  this  is  possible. 

The  first  order  expressions  for  stress  resultants  in  terms 
of  displacements  are  derived  in  Appendix  A  and  expressed  by  Equations 
A- 26 .  These  equations  uncouple  when  we  choose  C  as  the  centroid;  S 
as  the  shear  center;  B  as  a  sectorial  centroid;  and  the  x  and  y  axes 
as  principal  axes  of  cross  section  (see  Fig.  2.1).  Under  these  con¬ 
ditions  Equations  A-26  become 

=  E  A  w' 
c 


P 


(2-21a) 


. 


' 

' 


17 


M 

x 


-El  v" 
X  s 


(2--21b) 


M  =  -  E  I  u"  (2-21c) 
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to  to 

where  I  ,  I  and  I  are  the  principal  moments  of  inertia,  and  the 
x  y  to 

sectorial  moment  of  inertia,  respectively. 

Equation  2-16,  which  is  the  variational  equivalent  of  Equa¬ 
tions  2-20,  may  then  be  written  for  elastic  case,  by  substituting 
Equations  2-21  into  Equation  A-49 ,  which  yields  the  specialized 
equilibrium  equation 
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The  coefficients  C  ,  C  ,  C  and  I  are  properties  of  the  section,  as 
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derived  in  Appendix  A  (Equation  A-53) ,  and  d  ,  d  are  the  coordinates 

x  y 

of  the  load  axis  from  the  S-axis. 

It  is  emphasized  that  Equation  2-22  is  valid  only  for  linear 
elastic  response  where  C,  is  the  centroid,  S  is  the  shear  centre  and 
the  axes  are  the  principal  axes  of  the  cross  section. 


2 . 5  Comparison  of  Equations 

A  comparison  of  the  variational  Equation  2-22  with  the  spec¬ 
ialized  case  derived  by  Bleich  (13)  ,  reveals  that  he  omitted  the  terms 

(H  C  )<j>"  and  (M'C  )<j>'.  This  omission  has  already  been  pointed  out  by 
XX  XX 

Masur  and  Milbrandt  (49)  .  The  same  terms  were  omitted  by  Barsoum  and 

Gallagher  (7,8).  Substituting  Equation  2-21  into  Equations  2-20  the 

equilibrium  equations  derived  in  this  chapter  agree  with  Timoshenko 

(68)  and  Vlasov’s  (74)  equations,  except  that  they  have  neglected  the 

normal  stress  due  to  bimoments  in  the  M  term.  Oden  (54)  omitted  the 

P 

change  in  the  transverse  forces  and  torsional  moment  resulting  from  the 
rotation  and  translation  of  shearing  stresses  (as  shown  in  Fig.  2.3c, 

2.3d  and  2.4)  in  his  derivation  of  these  equations.  If  his  equations 
are  specialized  to  principal  axes  they  therefore  differ  from  Vlasov's 
and  those  derived  in  this  chapter  by  the  terms  arising  from  these  effects. 
It  would  be  necessary  to  make  numerical  comparisons  in  order  to  arrive 
at  conclusions  regarding  the  significance  of  these  terms  omitted  by 
other  authors. 

The  conceptual  difficulty  of  including  all  these  factors  in  an 
equilibrium  approach  demonstrates  the  desirability  of  investigating  the 
formulation  from  a  variational  point  of  view. 

2 . 6  Review  of  Finite  Element  Beam  Formulations 

Turner  (72),  et  al,  showed  that  a  'new  class' 
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matrices  had  to  be  introduced  into  the  equilibrium  equations,  when 
large  deflections  and  stability  were  to  be  investigated.  They  pre¬ 
sented  the  derivation  for  this  stiffness  matrix  for  an  axially  loaded 
member.  The  first  published  development  of  beam-column  stiffness 
matrix  was  by  Gallagher  (78)  and  Padlog.  Archer  (2)  also  has  dis¬ 
cussed  structural  stability  and  gives  an  energy  expression  for  the 
beam-column  which  is  virtually  a  verbation  restatement  of  Gallagher  . 
Arghyris  (3)  showed  that,  by  using  a  large  number  of  elements  to 
represent  a  nonuniform  column  excellent  convergence  could  be  achieved 
for  the  buckling  load,  in  spite  of  the  fact  that  nis  geometric  stiff¬ 
ness  matrix  was  inconsistent  .  Martin  (45)  arrived  at  the  initial 
stress  stiffness  matrices  for  axial  force  members  and  for  beam-column 
problems.  He  expressed  the  bending  displacement  as  a  cubic  polynomial 
and  the  axial  displacement  as  a  linear  function. 

Krajcinovic  (40)  developed  a  finite  element  formulation  for 
a  thin-walled  member  based  on  the  use  of  hyperbolic  functions  to 
represent  twist.  These  functions  represent  the  exact  displacement 
field  for  the  stable  formulation  and  yielded  excellent  results.  His 
formulation  cannot  account  for  the  effect  of  moment  gradient  on 
lateral  buckling.  Barsoum  (7,8)  and  Gallagher’s  formulation  is  similar 
to  Krajcinovic  except  that  the  effect  of  moment  gradient  can  be  accounted 
for  and  twist  is  expressed  in  terms  of  a  cubic  polynomial.  They  also 
investigated  fifth  order  Hermitian  polynomials  and  concluded  that  fifth 
order  polynomials  are  superior  to  cubic  polynomials  in  linear  stability 
problems.  They  recommended  the  use  of  lower  order  displacement  junctions 
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for  the  incremental  matrices.  They  also  concluded  that  the  desired 
accuracy  can  be  achieved  by  choosing  a  larger  number  of  elements 
when  dealing  with  nonlinear  problems.  In  addition,  Powell  (59) 
solved  the  problem  of  lateral  buckling  of  I-section  beams  using  the 
elastic  and  geometric  stiffness. 

The  formulation  developed  in  this  dissertation  uses  cubic 
polynomials  for  the  displacements  u,  v,  and  <£ ;  a  linear  function  for 


• 

axial  displacements  in  the  case  of  elastic  members  with  prismatic 
elements,  and  a  second  degree  axial  displacement  function  for  elastic 
members  with  non-prismatic  elements  and  for  all  inelastic  problems. 
Stiffness  matrices  are  derived  for  an  element  with  a  linear  variation 
of  geometric  properties.  Using  this  formulation,  it  is  possible  to 
carry  out  an  approximate  elastic  analysis  of  tapered  thin-walled 
beams  of  general  cross  section  subjected  to  general  loading  and 
support  conditions.  The  effect  of  residual  stress  is  also  included 
in  the  analysis.  The  resulting  equations,  for  the  elastic  case,  are 
similar  to  those  derived  by  Barsoum  and  Gallagher  (7,8)  who  used 
potential  energy,  rather  than  virtual  work,  for  their  formulation. 

2 . 7  Finite  Element  Idealization 

In  the  finite  element  approach  the  displacements  are  approxi¬ 
mated  by  forming  a  linear  combination  of  assumed  shape  functions, 
where  the  parameters  are  appropriate  nodal  displacements.  Let  (f^), 

{ f 2 } ,  and  (f  }  represent  the  cubic,  quadratic  and  linear  interpolating 
functions,  respectively.  Then  we  may  write  the  displacements  of  the 
reference  axes,  in  any  element  as 
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where  ug ,  vg ,  4> ,  w<;  represent  continuous  displacements  and  { u} , 

{v} ,  {£}  and  {wc}  represent  the  corresponding  nodal  displacement 
vectors,  as  defined  in  Fig.  2.6.  Substituting  Equations  2-23  into 
Equations  2--22,  we  obtain  for  one  element. 
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where  the  subscript  £  indicates  that  the  quantity  is  a  variable  with 
respect  to  the  z  coordinate. 

Assuming  all  the  geometric  properties,  the  axial  loads,  and 
the  moments  vary  linearly  along  the  element,  and  integrating  with 
respect  to  3  yields,  (see  Appendix  A,  Equations  A-54  and  A-55) 

[ks]  (rE}  +  [*G]  {r£}  =  (y  (2-25) 

where  (r  }  is  defined  in  Fig.  2.6; 

hi 

[kg]  is  the  usual  element  stiffness  matrix  for  the  linear  elastic  pro¬ 
blem  and  [k  ]  is  called  the  '’initial  stress"  or  "geometric  stiffness" 

G 

matrix. 

Equations  2-24  and  2--25  express  the  total  equilibrium  re¬ 
quirement  implied  in  Equations  2-11  and  2-15.  The  matrix  [k  ]  must  be 

G 

evaluated  using  the  stress  resultants  in  the  final  equilibrium  position, 

which  are  directly  dependent  on  the  displacements,  (r  }. 

h 

The  incremental  equilibrium  Equation  2-14  can  be  written  in 
finite  element  form,  neglecting  the  initial  displacement  matrix  (see 
Appendix  B) ,  as 


[kg]  { ArE>  -I-  [kg]  (Ay  =  (Ay 


(2-26) 


where  [lc  ]  and  [k  ]  have  exactly  the  same  form  as  in  Equation  2-25. 
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However  a  distinction  occurs  in  that,  for  Equation  2-26  the  geometric 
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stiffness  is  evaluated  from  the  stress  resultants  at  the  beginning  of 
the  increment. 

The  expressions  for  [k^]  an<^  [k  ]  are  derived  for  a  typical 
element  in  Appendix  A.  They  may  be  partitioned  in  the  following  form 
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2 . 8  Assembly  of  Finite  Element  Equations 

The  element  stiffness  matrices  of  Section  2-7  have  been 

evaluated  with  respect  to  nodal  displacements  referenced  to  a  'local' 

coordinate  system,  where  u  and  v  are  shear  center  displacements  in 

s  s 

the  directions  of  the  principal  axes;  w  is  the  uniform  axial  dis¬ 
placement  of  the  centroid;  and  (p  is  the  twist.  In  many  cases  it  is 
convenient  to  select  a  different  set  of  reference  axis  for  the  global 
system  of  nodal  displacements. 

Let  the  lateral  loads  be  applied  through  an  axis  passing 
through  an  arbitrary  point  'S'  of  the  cross  section  (not  the  shear 
center) ,  and  the  axial  load  be  applied  through  C  (not  the  centroid) . 
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Let  C  be  the  centroid  of  the  section  and  S  the  shear  center  and  let  the 
principal  axes  of  the  cross  section  be  at  an  angle  of  a  to  the  global 
axes,  as  shown  in  Fig.  2.7.  The  element  displacements  with  respect  to 
the  local  coordinate  system  can  be  related  to  those  in  the  global 
coordinate  system  by  the  transformation 
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where  [T]  is  given  by 
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;  the  submatrices  are  shown  in  Table  2-1; 
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/dU\  P  p  /c^v\  P  9 

dz  W  5 


Ap 

dz 


l  >  L 


(2-31b) 


and 


di  p 

^  >G 


(2-31c) 


;  and  p  and  q  are  nodal  numbers. 


25 


Rewriting  Equation  2-25  in  this  nomenclature  yields 

[ks]  {rE}L  +  [kG]  {rE}L  =  {Re}L  (2-32) 

Suostituting  for  in  terms  of  and  premultiplying  by  [t]^ 

yields , 

[T]T  [ks]  [T]  { rE>G  +  [T]1  [kG]  [T]  {rE)G  =  {Re>g  (2-33a) 
or 

[kS]G  {rE}G  +  [kG]  {rE}G  “  {REJG  (2“33b> 

The  formation  of  the  complete  stiffness  matrix  for  the  entire  structure 
is  obtained  by  direct  addition  for  all  the  interface  nodes  (thus  satis¬ 
fying  the  equilibrium  and  compatibility  requirements  (35) ,  and  finally 
introducing  the  kinematic  constraints.  The  stiffness  matrix  for  the 
entire  structure  can  be  written  as, 

[Kg]  {r}  +  [Kg]  {r}  =  (R)  (2-34) 

Where  [K  ]  is  the  elastic  flexural  stiffness  matrix;  [K  ]  is  the 

b  (j 

geometric  or  initial  stress  stiffness  matrix;  and  {r}  is  the  assembled 
vector  of  global  nodal  displacements. 

Equation  2--34  is  the  finite  element  form  of  the  total  equili¬ 
brium  equations  as  expressed  in  Equation  2-15  and  evaluated  in  Appendix 
A.  The  finite  element  form  of  the  incremental  equilibrium  equations, 
Equation  2-14,  with  the  simplifying  assumptions  of  Appendix  B,  may  be 
written  as 


[K J  {Ar}  +  [K J  {Ar}  -  (AR> 


(2-35) 


. 


' 
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Where  [K  ]  and  [K  ]  are  identical  to  those  of  Equation  2-34  but  [K  ] 
b  b  G 

is  determined  from  initial  stress  prior  to  the  load  increment. 

2 .9  Summary 

These  equations  up  to  and  including  those  in  Section  2.3, 
apply  whether  the  member  is  elastic  or  inelastic.  In  the  remainder 
of  the  chapter,  the  equations  were  specialized  for  the  elastic  case 
and  a  finite  element  formulation  for  this  case  was  developed. 
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TABLE  2-1  TRANSFORMATION  MATRIX 
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FIGURE  2.1  CROSS  SECTION  OF  A  THIN-WALLED  BEAM 


FIGURE  2.2  THIN-WALLED  BEAM  OF  OPEN  SECTION-  A  GENERAL  SYSTEM 
OF  LOADS  IN  A  TYPICAL  BEAM  SEGMENT 
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FIGURE  2.3b  PROJECTIONS  OF  A  DEFORMED  BAR  ON  YZ  &  XZ  PLANE 
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FIGURE  2.3c  INTERNAL  FORCES  INDUCED  DUE  TO  SHEARING  DEFORMATION 
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FIGURE  2.3d  INPLANE  FORCES  INDUCED  BY  NORMAL  AND  SHEARING  STRESSES 
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FIGURE  2.4  ADDITIONAL  TORSIONAL  MOMENT  INDUCED  DUE  TO  SHEAR 
STRESSES  DURING  DEFORMATION 
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are  increments  in  stresses,  strains, 
displacements  and  surface  leads 
respectively 


FIGURE  2.5  STRESS  STRAIN  AND  DISPLACEMENT  NOMENCLATURE  IN 
INITIAL  AND  INCREMENTED  CONFIGURATION 
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<  w  q  ,  wq  >  For  Elastic  prismatic  elements 
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FIGURE  2.6  REFERENCE  AXES,  NODAL  DISPLACEMENTS  AND  STRESS  RESULTANTS 
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FIGURE  2.7  THIN-WALLED  BEAM  WITH  RESPECT  TO  LOCAL  AND  GLOBAL  ORIENTATION 


CHAPTER  III 


SOLUTION  OF  ELASTIC  BEAMS  AND  BEAM-COLUMNS 

3 . 1  Introduction 

In  this  chapter  the  elastic  finite  element  formulation 
developed  in  Chapter  II  is  applied  to  a  number  of  different  types  of 
problems.  The  problems  may  be  separated  into  three  categories,  namely; 
a)  small  deflection  problems  involving  the  solution  of  first  order 
beam  equations;  b)  buckling  or  bifurcation  problems  involving  the 
eigenvalue  solution  of  the  incremental  equilibrium  equations;  and  c) 
beam-column  type  problems  which  include  the  effect  of  "prebuckling" 
deformations  by  satisfying  the  nonlinear  total  equilibrium  equations.  - 

The  study  of  beam-column  response  in  the  elastic  range  is 
motivated  by  a  number  of  considerations.  The  elastic  critical  or 
'buckling*  load  forms  an  upper  bound  on  the  carrying  capacity  of  such 
a  member.  However  if  the  design  criteria  is  controlled  by  a  deflection 
limitation  then  it  is  necessary  to  determine  the  loading  condition  at 
which  this  limitation  is  reached  and  useful  to  relate  this  load  to  the 
upper  bound  on  carrying  capacity.  The  range  of  loading  for  wThich  an 
elastic  analysis  is  applicable  is  limited  by  the  proportional  limit  of 
the  material  and  the  actual  carrying  capacity  is  normally  governed  by 
inelastic  beam-column  response.  It  is  useful,  however,  to  study  the 
elastic  beam-column  response  as  a  prerequisite  to  solving  the  more 
realistic  inelastic  problem. 

It  might  also  be  noted,  that,  with  the  current  development 
of  high  strength  materials,  the  range  of  loading  for  which  elastic 
analysis  is  applicable  is  increasing.  Since  the  buckling  problem  is 
based  on  a  classical  eigenvalue  analysis,  it  ignores  all  deformations 
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prior  to  the  attainment  of  the  critical  load  and  these,  together  with 
geometric  misalignments,  may  make  the  eigenvalue  estimate  an 
unrealistically  high  upper  bound.  In  this  chapter,  the  finite  element 
model  derived  in  Chapter  II  is  used  to  study  both  elastic  beam-column 
behaviour  and  the  classical  eigenvalue  or  buckling  analysis. 

The  method  is  also  applied  to  beam-column  problems  with 
initial  geometric  imperfections.  For  this  case  Southwell  (65)  pro¬ 
posed  an  extrapolation  technique  by  means  of  which  experimental 
elastic  critical  loads  can  be  predicted  for  concentrically  loaded 
pin-ended  columns,  with  sinusoidal  imperfections.  Gregory  (29)  and 
Ariaratnam  (6)  have  extended  the  above  method  in  the  case  of  struts 
which  fail  by  flexural  torsional  buckling.  The  elastic  deformation 
of  I-beams  with  sinusoidal  initial  crookedness  or  twist  under  the 
action  of  equal  and  opposite  end  moments  has  been  carried  out  by 
Flint  (20)  and  Massey  (47)  for  simply  supported  beams,  and  by  Horne 

(36)  and  Trahair  (70) (71)  for  beam-columns.  Recently  Leicester  (42) 
has  shown  the  applicability  of  the  Southwell  plot  to  monosyrametrical 
be am- columns . 

3 . 2  Solution  of  First  Order  Beam  Equations 
3.2.1  Introductory  Remarks 

The  finite  element  model  derived  in  Section  2.7  has  been 
used  to  represent  beam  flexure  by  a  number  of  investigators.  Barsoum 
and  Gallagher  (7,8)  and  Powell  et  al  (59)  have  included  twist  in  their 
model  using  the  same  interpolation  functions.  However,  neither  of 
these  studies  investigated,  in  detail,  the  ability  of  the  model  to 
represent  non-uniform  torsional  response. 
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The  bimoment  which  arises  in  the  derivation  of  the  beam 

equation  has  been  designated  as  W  .  This  bimoment  is  produced  by 

co 

restraint  of  warping  deformations  and  the  rate  of  change  of  bimoment 
produces  a  torque  along  the  beam  called  the  "warping  torque",  M  . 

A  measure  of  magnitude  for  bimoment  is  given  by  the  product  of  the 
distance  between  two  parallel  planes  and  the  moment  on  one  of  them. 

A  bimoment  is  assumed  positive  when  the  sense  of  one  moment  seen  from 
the  plane  of  the  other  moment  is  anticlockwise.  A  rigorous  defini¬ 
tion  is  given  in  Appendix  A,  Equation  A-14.  The  distinctive  feature 
of  the  bimoment  is  that  when  it  is  acting  in  an  elastic  body  it  pro¬ 
duces  internal  stresses  and  strains  which  have  conventional  stress 
resultants  of  zero  magnitude.  When  any  beam  is  cut  by  an  imaginary 
plane  the  bimoment  acting  on  the  section  to  one  side  of  the  plane  is 
equal  and  opposite  to  that  acting  on  the  other  side. 

One  type  of  problem  that  a  model  must  be  capable  of  analys¬ 
ing,  to  produce  reliable  results  for  thin-walled  beams  of  open  section, 
is  the  problem  of  nonuniform  torsion,  where  the  total  torque  is 
resisted  by  a  combination  of  the  warping  torque,  M^,  and  the  St.  Venant 
torque,  .  The  following  examples  test  the  capability  of  the  model 
to  represent  stress  conditions  in  beams  subjected  to  nonuniform  torsion. 

3.2.2  Numerical  Results  for  Nonuniform  Torsion  Problems 
E-l  Nonuniform  Torsion  of  a  Prismatic  Beam 

In  order  to  check  the  ability  of  the  finite  element  model  to 
represent  nonuniform  torsional  behaviour,  an  8WF31  cantilever  beam  of 
50  inches  length  was  analyzed  with  a  5  inch-kip  torque  applied  at  the 
free  end.  Eight  elements  were  used  and  the  finite  element  solution  was 
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compared  with  the  closed  fora  solution  (68).  Fig.  3.1  shows  the  varia¬ 
tion  of  the  angle  of  twist  <j) ,  and  its  derivatives,  along  the  length. 

In  Fig.  3.2  the  torsional  moment  of  resistance  due  to  St.  Venant  and 

warping  are  compared  with  the  closed  fora  solution  (68)  .  Since  a  cubic 

polynomial  was  used  for  twist,  the  variation  of  warping  torsional 
moment  of  resistance  is  constant  along  the  element  whereas  the  varia¬ 
tion  of  St.  Venant  resistance  is  quadratic.  Fig.  3.3  shows  a  compari¬ 
son  of  the  twist  curves  along  the  length  of  the  beam  for  different 
numbers  of  elements.  The  convergence  characteristics  of  the  solution 
are  shown  in  Fig.  3.4. 

E-2  Nonuniform  Torsion  of  a  Tapered  I-Beam 

Fig.  3.5  shows  the  comparison  of  the  finite  element  analysis 
with  Lee's  (41)  power  series  solution  for  nonuniform  torsion  of  tapered 
beams.  The  finite  element  analysis  agrees  with  the  closed  form  solu¬ 
tion  (68)  for  a  uniform  beam.  A,  but  this  solution  does  not  agree  with 
Lee’s  results.  For  the  tapered  beams,  B  and  C,  Lee’s  analysis  includes 
the  torque  about  the  longitudinal  axis  resulting  from  the  projection  of 
the  inclined  flange  moments.  This  effectively  reduces  the  torsional 
stiffness  of  the  beam.  Such  an  effect  is  not  incorporated  in  the 
equilibrium  equations  of  Chapter  II  and  therefore  cannot  be  represented 
by  the  finite  element  model.  When  the  beam  taper  is  considerable,  as 
in  these  examples ,  the  contribution  of  flange  moment  to  the  torsional 
moment  is  significant  and  hence  there  is  a  discrepancy  between  Lee’s 
analysis  and  the  results  from  the  conventional  beam  equations. 

3 . 3  Elastic  Buckling  Problems 
3.3.1  Introductory  Remarks 
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Many  structures  when  loaded  to  a  critical  state  will  undergo 
a  marked  change  in  the  character  of  deformation,  which  is  not  the 
result  of  material  failure  or  alteration  of  mechanical  properties. 

Such  a  change  occurs  because  one  mode  of  deformation  becomes  unstable 
and  the  structure  seeks  a  stable  configuration.  The  change  in  deforma¬ 
tion  which  takes  place  when  the  structure  passes  from  an  unstable 
configuration  to  a  stable  configuration  is  generally  known  as  buckling. 
The  load  at  which  this  occurs  is  known  as  the  critical  or  bifurcation 
load.  The  classical  buckling  problem  ignores  all  deformations  prior 
to  the  attainment  of  the  buckling  load  and  hence  the  eigenvalue 
analysis  always  forms  an  upper  bound  for  the  load. 

3.3.2  Solution  Scheme 

A  buckling  problem  is  formulated  by  imposing  a  variation  of 
displacements  in  the  matrix  equation  2-34  under  constant  load.  The 
result  is  identical  to  the  homogeneous  form  of  Equation  2-35,  i.e.:  - 

[[KS]  {Ar}  +  [KG]  (Ar}J  =  {0}  (3-1) 

The  above  equation  is  satisfied,  for  a  non-zero  displacement  vector, 
only  if  the  stiffness  matrix  is  singular,  and  the  "buckling  criterion" 
can  be  represented  by  the  vanishing  of  the  determinant  of  the  stiff¬ 
ness  matrix,  i.e.:  - 

| | [Ks]  +  [KG] ||=0  (3-2) 

By  assuming  that  throughout  the  loading  range  the  stress  state  can  be 
written  in  terms  of  some  function  whose  intensity  is  governed  by  a 
single  parameter,  A,  the  initial  stress  state  takes  the  form 


. 
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[KG]  =  A[Kg]  (3-3) 

where  [Kg]  represents  the  initial  stress  stiffness  matrix  computed  for 
some  nominal  stress  distribution.  Introducing  Equation  3-3  into 
Equation  3-1  yields 

[Ks]{ Ar}  =  -A[KG]{Ar}  (3-4a) 

or 

X  {Ar}  =  -[Kg]-1  [KG]{Ar>  (3-4b) 

The  critical  buckling  load  is  now  obtained  by  solving  Equation  3-4b 
for  the  unknown  eigenvalue  A.  The  eigenvalue  solution  for  the  follow¬ 
ing  problems  was  determined  using  "NROOT"  of  the  IBM-SSPackage  (1). 

The  buckling  mode  can  be  identified  as  the  eigenvector  corresponding 
to  the  lowest  eigenvalue. 

3.3.3  Numerical  Examples 

In  order  to  verify  the  formulation  a  wide  variety  of  problems 
have  been  solved  and  checked  with  closed  form  solutions.  Table  3-1 
shows  the  comparison  of  the  results  of  a  large  number  of  critical  loads 
computed  by  the  finite  element  analysis  with  the  corresponding  avail¬ 
able  solutions.  The  table  includes  examples  of  Euler,  lateral  and 
lateral  torsional  buckling.  Nomenclature  is  specified  in  Fig.  3.6. 

Some  important  additional  examples  are  described  below. 

E-l  Torsional  Buckling  of  a  Symmetrical  I-Beam 

Table  3-2  shows  a  comparison  of  the  results  of  Krajcinovic 
(40)  and  the  present  formulation.  Krajcinovic  used  a  hyperbolic  func¬ 
tion  for  the  twist  and  claimed  that  the  results  from  a  cubic  polynomial 
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would  not  be  able  to  reflect  the  terra  k  (ratio  of  St.  Venant  torsional 
rigidity  to  warping  rigidity)  no  matter  how  many  elements  were  used. 

In  the  present  formulation  a  four  element  approximation  gives  results 
closer  to  the  exact  solution. 

E-2  Influence  of  Point  of  Load  Application  for  Lateral  Buckling  of 
Beams 


Fig.  3.7  shows  the  solution  for  lateral  buckling  of  a  beam 
with  a  single  vertical  concentrated  central  load  applied  at  a  distance 
dy  above  or  below  the  shear  center.  The  section  is  that  used  by 
Powell  (59),  et.  al. ,  and  is  a  21"  x  10"  I  beam  with  a  uniform  thick¬ 
ness  of  1".  The  span  is  300".  The  solution  is  compared  to  that  of 
Timoshenko  (68) . 

E-3  Elastic  Buckling  of  Indeterminate  Structures 

Table  3-3  shows  the  comparison  of  finite  element  results  with 
those  of  Horne  and  Ajmani  (37)  who  developed  a  series  solution  to 
determine  the  critical  loading  on  columns  supported  laterally  by  "side 
rails".  The  side  rails,  acting  with  the  cladding  and  any  bracing 
system,  are  assumed  to  provide  a  rigid  lateral  support  at  the  points 
of  attachment.  They  also  provide  an  elastic  torsional  restraint.  A 
8WF  28  column,  69"  in  length  and  restrained  laterally  by  side  rails  at 
a  spacing  of  23"  and  at  an  eccentricity  of  3"  from  the  centre  line  of 
the  column,  was  analyzed  when  subjected  to  an  axial  load  of  P  and  a 
uniform  moment  of  M  about  the  major  axis.  At  each  rail  position  the 
torsional  restraint  was  assumed  to  be  1040  in.K/Rad.  As  shown  in 
Table  3-3,  the  load  which  causes  lateral  torsional  buckling  is  1972.19 
kips,  when  the  rotations  are  nonzero  at  lateral  supports  (case  8), 


and  6100  kips  when  the  rotations  are  zero  at  lateral  supports  (case  17) . 


. 
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In  the  first  case  the  column  buckles  in  a  lateral  torsional  buckling 
mode  and  in  the  second  case  the  column  buckles  in  an  Euler  mode.  The 
mode  shapes  for  these  cases  are  shown  in  Fig.  3.8. 

E-4  Euler  Buckling  of  an  Unsymmetrical  Section 

The  unsymmetrical  beam  used  by  Vlasov  (74)  for  unsymmetrical 
bending, shown  in  Fig.  3.9,  was  analyzed  with  an  axial  load  applied  at 
the  centroid  and  with  a  span  of  300".  The  eigenvalue  corresponds  to 
lateral  torsional  buckling.  The  finite  element  critical  load  is 
1163  kips  compared  to  a  value  of  1140  kips  determined  from  the  equations 
of  Timoshenko  (68). 

E- 5  Lateral  Torsional  Buck ling  of  an  Unsymmetrical  Section  Under 
Uniform  Moment 

End  moments  were  applied  about  the  X-axis  for  the  beam  of 

example  E-4  and  the  critical  moment  to  cause  lateral  torsional  buckling 

4 

about  the  minor  axis  was  found  to  be  2.3706  x  10  inch-kips  compared 

4 

to  a  value  of  2.367  x  10  inch-kips  determined  by  the  closed  form 
solution  (68) . 

E-6  Lateral-Torsional  Buckling  of  an  Unsymmetrical  Section  with 

Moment  Gradient 

A  vertical  load  was  applied  at  the  mid-span  of  the  beam  of 
example  E-4  and  E-5  acting  2.42"  above  the  shear  centre  as  shown  in 
Fig.  3.9.  The  critical  load  to  cause  lateral-torsional  buckling  was 
found  to  be  equal  to  263.27  kips.  A  classical  solution  of  this  pro¬ 
blem  involves  the  solution  of  three  coupled  ordinary  differential 
equations  and  no  closed  form  comparison  is  available. 

E-7  Elastic  Buckling  of  a  Tapered  Axially  Loaded  Member 


A  tapered  WF  column  with  the  dimensions  as  shown  in  Fig.  3.10 
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was  analyzed  with  the  small  end  pinned  and  the  larger  end  fixed.  The 
finite  element  analysis  estimates  the  buckling  load  4.93  percent 
higher  than  that  given  by  Culver  (17) ,  who  has  published  a  set  of 
curves  for  this  type  of  problem  using  Lee’s  (41)  nonuniform  torsion 
formulation. 

E - 8  Lateral  Torsional  Buckling  of  a  Tapered  Beam  Under  Moment  Gradient 

The  beam  described  in  E-7  was  loaded  with  a  moment  at  the 
larger  end  and  the  critical  moment  was  found  to  be  1035  inch-kips 
which  is  9.21  percent  lower  than  the  critical  moment  given  by  Culver 
(17) .  The  discrepancy  is  attributed  to  the  component  of  the  flange 
bending  moments  in  the  z-direction,  which  were  considered  by  Lee  and 
Culver,  and  discussed  in  example  E-2  of  Section  3.2.2,  and  which  are 
not  taken  into  account  in  present  analysis. 

E-9  Interaction  Surface  for  Critical  Loading 

Fig.  3.11  shows  the  influence  buckling  surface  for  an  8WF28 
beam,  69"  in  length,  when  it  is  loaded  by  axial  load  and  end  moments. 

3 . 4  Solution  Methods  for  Beam-Column  Problems 
3.4.1  Introductory  Remarks 

The  previous  examples  in  this  chapter  have  dealt  with  the 
small-deflection  beam  equations  and  elastic  buckling  or  bifurcation 
problems.  When  pre-buckling  deformations  are  considered  it  is 
necessary  to  solve  a  ’beam-column’  type  problem.  This  involves  the 
solution  of  Equation  2-34  which,  because  of  the  dependence  of  [K  ] 
on  {r},  is  nonlinear.  There  are  a  number  of  methods  of  solving  this 
set  of  nonlinear  equations.  A  ’Newton-Raphson ’  method  was  used  to 
solve  elastic  beam-column  problems  and  a  ’Modified  Newton ’Raphson' 
was  used  to  solve  inelastic  beam-column  problems  in  this  dissertation. 
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A  description  of  these  two  methods,  applied  to  elastic  problems,  is 
given  in  the  remainder  of  this  section.  An  equilibrium  check  is 
incorporated  into  both  techniques. 

3.4.2  Newton-Raphson  Method 

This  method  has  been  employed  successfully  by  a  number  of 
investigators.  For  example.  Walker  and  Hall  (75)  used  it  to  study 
large  deflections  of  beams,  and  Brebbia  (14)  and  Connor  used  it  to 
study  the  stability  and  geometrically  non-linear  behaviour  of 
arbitrary  shells.  Initially,  a  small  displacement  solution  of  Equa¬ 
tion  2-34  is  obtained,  for  the  first  load  increment,  with  [K  ]  set  to 
zero.  The  element  stress  resultants  are  computed  from  the  known 
displacements  and  the  geometric  stiffness  matrix  for  the  structure  is 
evaluated.  The  total  resisting  force  at  any  node  is  the  sum  of  the 
resisting  forces  obtained  from  the  flexural  stiffness  and  the 
geometric  stiffness.  The  difference  between  the  resisting  forces  and 
the  applied  loads  represents  a  set  of  unbalanced  forces  for  the  given 
configuration  of  the  model.  Knowing  the  set  of  unbalanced  forces  on 
the  model,  in  this  configuration,  one  can  solve  for  the  increments  in 
the  nodal  displacements.  The  iteration  is  repeated  till  the  configura¬ 
tion  of  the  model  maintains  equilibrium  with  the  applied  loads.  The 
following  algorithm  explains  the  procedure  and  is  illustrated  in 
Fig.  3.12a. 

1.  For  any  approximate  {r}^,  compute  [k^,]^  aud  the  unbalanced  force 
vector , 

{AR>n+l  =  {Rl  -  CKS]{r}n  *  [KG]n  {r)n 


(3-5) 


' 


. 
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2.  Compute  the  increment  in  displacements  to  equilibrate  the  unbalanced 
forces , 


{Ar} 

n 


[Ks]  + 


(3-6) 


3.  Update  the  displacements 


{r}  =  {r}  .  +  {Ar}  (3-7) 

n  n-1  n 

4.  Repeat  steps  1  to  3  until  the  unbalanced  force  vector  is  arbitrary 
small. 

5.  Apply  a  new  load  increment  and  repeat  steps  1  to  4. 

The  principal  disadvantage  of  this  technique  is  that  it  is  necessary 
to  triangularize  a  new  stiffness  matrix  for  every  iterate. 

3.4.3  Iterative  Initial  Force  Method  (Modified  Newton-Raphs on  Method) 
This  method  retains  the  linear  elastic  matrix  on  the  left 
hand  side  of  the  equation  and  keeps  all  the  nonlinear  terms  on  the 
right  hand  side.  It  is  illustrated  in  Fig.  3.12b.  For  the  nth  iterate. 
Equation  2-34  can  be  written  as, 

[KS]{r}n+l  -  tR)  “  £KG](r>n  (3~8) 

and  solved  for  the  next  estimate  of  nodal  displacements 

(r,n+l  -  Cqd'1  {R}  -  [Ksrl  [KG]{r>„  (3~9) 

The  above  procedure  is  repeated  till  the  displacement  for  the  (n+l)^1 
iterate  is  the  same  as  for  the  n^1  iterate,  at  which  time  the  unbal¬ 
anced  force  vector  is  arbitrarily  small.  The  matrix  [K^]  is  updated 
only  at  the  beginning  of  each  load  step. 
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To  investigate  the  convergence  of  this  method  one  may  follow 
the  procedure  which  Dupuis  et  al  (79)  used  for  direct  iteration: 


(r)  =  [K  ]-1{R}  -  [K„]  1[kp'J{[kJ  1{R}  -  [Kj  X[k  ]{r} 


n+1 


G  r  S' 


GJt  n-1 


(3-10) 


continuing  this  back  substitution  yields 


{r>n+l  =  ^  +  "  Cg ] 3  . (-l)n  [G]n  ) 


x 


[K q]”1{R}  +  (-l)n+1  [G]n+1  {r  } 
b  o 


(3-11) 


where 


[G]  =  [Kg]"1  [Kg] 


(3-12) 


Assuming  the  displacements  are  initially  zero  the  convergence  depends 
on  the  expression 


(-1)°  [G]n  [Kg]’1  {R} 


tending  to  zero 


If  {u^.}  are  the  eigenvectors  of  the  matrix  [g],  we  may  write 


[Kg]  1  (R)  =  I  C  {u  } 

j=l  J 


n 


(3-13) 


Noting  that  [K„]  and  [K  ]  are  symmetric  matrices  and  [K  ]  is  a  positive 
b  '  G  b 

definitive  matrix,  we  have 


~  {u}  =  [g] {u} 


(3-14) 


[G]n  (u)  =  -  {u} 

A 


(3-15) 


[G]n  [Kg] 1  {R}  =  [G ] n  C-,  {u.,  }  +  [G]“  C9{u9}  + 


=  C1{u1}/X 


r  r 

n 


2  2 


(3-16) 
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Since  X  is  greater  than  1  when  the  current  load  is  less  than  the  buck¬ 
ling  load,  the  method  converges. 


3.4.4  Solution  Scheme  for  Beam-Columns  with  Geometric  Misalignment 
\Then  the  beam  has  an  initial  imperfection  of  {r^}  Equation 
2-34  can  be  written  as 


[Kj{r}  +  [K  ]{ r }  =  {R}  -  |X]{r  } 

S  tj  Go 


(3-17) 


The  unbalanced  load  for  the  Newton-Raplison  method 


{AR} 


n 


{R}  -  [K  ]{r)  -  [K  ]  {r}  -  [K  ]  {r  }  (3-18) 

S  n  Gnn  Gno 


and  the  corrections  {Ar}^^  is  given  by  the  solution  to 


[[KS]  +  [KG]J{Ar}n+l  “  {AR}n  (3“19) 

where  { r }  are  the  displacements  from  the  initially  imperfect  configura¬ 
tion  . 


3 . 5  Elastic  Beam-Column  Response 
3.5.1  Introductory  Remarks 

For  beam-column  type  problems  the  load-deformation  response 
was  obtained  by  applying  the  load  incrementally  and  iterating  to  find 
the  equilibrium  configuration  for  each  loading  condition  as  described 
in  Section  3.4.2.  A  solution  to  this  type  of  problem  can  be  obtained 
for  statically  indeterminate,  as  well  as  statically  determinate  mem¬ 
bers,  subjected  to  any  combination  of  externally  applied  concentrated 
or  distributed  loads,  moments,  or  torques,  with  or  without  elastic 


constraints . 
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3.5.2  Examples  of  Beam-Column  with  Initial  Imperfection 
E-l  Inplane  response  of  a  column  (axial  load) 

Fig.  3.13  compares,  with  the  solution  of  Timoshenko  (68), 
the  column  response  for  an  8WF31  section  of  length  69"  with  an  initial 

7T  Z 

deflection  of  0.01  sin  — — ,  when  subjected  to  an  axial  load.  Deflec¬ 
tion  increases  in  a  nonlinear  manner  and  after  0 . 8P  the  deflection 

cr 

increases  rapidly. 

E-2  In  and  out  of  plane  response  of  a  column 

In  Fig.  3.14  the  load  deflection  response  of  a  beam-column 
is  shown,  for  initial  deflections  of  -0.01  sin  ttz/L  in  u,  v,  and  <j> , 
when  the  member  is  subjected  to  axial  load.  Mien  the  load  reaches 
the  critical  load  corresponding  to  the  Euler  buckling  load  about  the 
minor  axis,  the  sign  of  the  determinant  changes.  For  loads  above  the 
critical  value,  the  sign  of  the  u-def lection  changes.  Mien  the  load 
reaches  the  torsional  buckling  load,  the  sign  of  the  determinant 
changes  again  as  well  as  the  sign  of  the  twist.  The  v-deflection 
continues  to  increase  and  becomes  asymptotic  to  the  line  of  the  major 
axis  Euler  buckling  load. 

E - 3  Out  of  plane  response  of  a  beam  due  to  end  moments 

Fig.  3.15  shows  the  moment  and  out  of  plane  response  for  an 

TT  Z 

8WF31  beam  when  there  is  an  initial  lateral  deflection  of  0.01  sin  — — . 

JL 

Moment  and  deflection  curves  are  drawn,  neglecting  the  effect  of 
moment  about  the  minor  axis  on  the  equilibrium  equation  for  torsional 
moment,  and  compared  with  the  results  obtained  by  Trahair  (71). 

Curves  are  also  drawn  when  this  coupling  effect  is  considered.  Closed 
form  and  finite  element  solutions  are  compared  for  both  cases. 


, 

' 
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E~4  Out  of  plane  response  of  beam-column  due  to  axial  load  and  equal 

end  moments 

A  12WF50  with  a  sinusoidal  imperfection  is  analyzed  with 
the  axial  load  applied  at  different  eccentricities.  The  beam-column 
response  curves  are  shown  in  Fig.  3.16  and  compared  with  corresponding 
eigenvalue  solutions. 

3 . 6  Summary 

In  this  chapter,  some  examples  of  solutions  were  presented 
for  a)  the  first  order  beam  equations,  b)  elastic  buckling  of  both 
statically  determinate  and  indeterminate  structures,  and  c)  the  elastic 
analysis  of  beam-columns  with  imperfections. 
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TABLE  3-1  COMPARISON  OF  FINITE  ELEMENT  ANALYSIS  WITH  TIMOSHENKO  (68) 
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TABLE  3-3  COMPARISON  OF  FINITE  ELEMENT  ANALYSIS  WITH  REPORTED  RESULTS 
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DISTANCE  OF  SECTION  FROM  FIXED  END  INCHES 


FIGURE  3.1  NON-UNIFORM  TORSION  OF  WF  BEAM 


MOMENT  RESISTED  BY  ST. VERNANT  SHEAR 
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FIGURE  3.2  MOMENT  OF  RESISTANCE  BY  ST.VENANT  AND  WARPING  SHEAR 
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DISTANCE  OF  THE  SECTION  FROM  FIXED  END 


FIGURE  3.3  NON-UNIFORM  TORSION  -  COMPARISON  WITH  DIFFERENT  ELEMENTS 
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FIGURE  3.4  CONVERGENCE  CHARACTERISTICS 
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FIGURE  3.5  NON-UNIFORM  TORSION  OF  TAPERED  BEAMS 
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CODE  FOR  BUCKLING  MODE 

EW  -  Euler  buckling  about  weak  axis 
ES  -  Euler  buckling  about  strong  axis 
LW~  Lateral  buckling  about  weak  axis 
LS  -  Lateral  buckling  about  strong  axis 
T  -  Torsional  Buckling 
LTB  Lateral  torsional  buckling 


FIGURE  3.6  NOMENCLATURE  FOR  SUPPORT  CONDITION  AND  BUCKLING  MODE 
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FIGURE  3.7  LATERAL  BUCKLING  OF  A  BEAM  WHEN  LOAD  IS  ACTING  AT  d  FROM  SHEAR-CENTER 
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LATERAL  TORSIONAL  BUCKLING  MODE  WHEN  ROTATIONS 
ARE  NON-ZERO  AT  LATERAL  SUPPORTS 


EULER  BUCKLING  ABOUT  THE  MAJOR  AXIS  WHEN 
ROTATIONS  ARE  ZERO  AT  LATERAL  SUPPORTS 


FIGURE  3.8  COLUMNS  SUPPORTED  LATERALLY  BY  SIDE  RAILS 
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NORMALIZED  WARPING  COORDINATE 
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FIGURE  3.9  ELASTIC  BUCKLING  OF  UNSYMMET RI CAL  SECTION 
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FIGURE  3.10  BUCKLING  OF  TAPERED  BEAM-COLUMN 
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FIGURE  3.11  INFLUENCE  BUCKLING  SURFACE  FOR  WIDE  FLANGE  BEAM 
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FIGURE  3.13  COLUMN  WITH  INITIAL  IMPERFECTION 
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FIGURE  3.14  LOAD  DEFLECTION  RESPONSE  OF  AN  'INITIALLY  IMPERFECT  COLUMN 
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FIGURE  3.15  BEAM  WITH  INITIAL  IMPERFECTION 
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FIGURE  3.16  COLUMN  WITH  SINUSOIDAL  IMPERFECTION 


CHAPTER  IV 


LOCAL  BUCKLING  OF  ELASTIC  SECTIONS 

4 . 1  Introduction 

Thin-walled  beams  are  made  up  of  relatively  thin  plates,  and 
under  certain  conditions,  these  plates  may  experience  local  buckling. 
Local  buckling  in  combination  with  other  buckling  behaviour  is  often 
observed  in  failure  of  members.  Local  buckling  of  the  various  plate 
components,  of  which  heavy  structural  sections  consist,  needs  rarely 
to  be  considered,  because  these  sections  are  proportioned  so  that  they 
will  not  buckle  locally  at  stresses  below  the  yield  point.  However, 
with  the  introduction  of  high  strength  steels,  sections  which  had 
previously  been  immune  to  local  buckling  prior  to  yield  may  no  longer 
be  so.  In  addition  cold  rolled  steel  shapes,  and  aluminum  shapes 
often  have  proportions  such  that  local  buckling  may  occur. 

Buckling  theories  of  thin-walled  members  have  been  separately 
developed  for  member  and  local  buckling.  Member  buckling  is  based  on 
the  assumption  that  the  cross  section  does  not  distort  while  local 
buckling  naturally  involves  a  distortion  of  the  section.  In  practice 
these  phenomena  may  not  be  independent  and  hence  it  is  reasonable  to 
investigate  conditions  where  both  may  occur  simultaneously. 

Bleich  (13)  studied  the  flexural  torsional  buckling  of 
T-shaped  stiffners  including  the  deformation  of  the  web.  Goldberg  (27), 
et.  al.,  presented  a  systematic  buckling  analysis,  for  members  with 
arbitrary  cross-sectional  shape,  considering  the  cross-sectional  de¬ 
formation,  by  starting  with  the  usual  plate  and  membrane  equations, 
and  arriving  at  eight  first  order  simultaneous  partial  differential 
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equations.  Their  procedure  consists  of  selecting  a  trial  value  of 
axial  load  and  moment  and  evaluating  the  determinant  of  the  final 
equations.  The  determinant  will  be  zero  if  the  trial  values  are 
chosen  correctly.  If  not,  a  new  set  of  trial  values  is  chosen.  Re¬ 
sults  of  some  specific  examples  were  presented.  For  the  shapes  of 
sections  which  they  studied,  local  buckling  was  critical  for  the  lower 
slenderness  ratios  and  occurred  at  a  much  lower  stress  than  member 
buckling.  Ko lib runner  (39),  et.  al.,  treated  the  thin-walled  beam 
using  folded  plate  theory  of  structures.  They  arrived  at  the  solution 
by  solving  a  system  of  differential  equations. 

In  this  chapter,  a  general  procedure  for  deriving  virtual 
work  equations  is  employed,  which  is  suitable  for  determining  the 
flexural  and  geometric  stiffnesses,  considering  distortion  of  the 
cross-section.  The  matrices,  however,  are  evaluated  only  for  elastic 
prismatic  symmetric  wide  flange  sections. 

4 . 2  Assumptions 

In  addition  to  the  assumptions  made  in  Chapter  II  the  follow¬ 
ing  assumptions  are  also  made. 

1.  The  original  angle  between  the  adjacent  plates  at  their  common  edge 
is  maintained  during  buckling. 

2.  Kirchhoff's  hypothesis  remains  valid  for  each  plate  segment. 
Assumption  2  of  Section  2.2  is  no  longer  applicable. 

4 . 3  Equilibrium  equations  by  virtual  work  including  distortion  of  the 
section 

4.3.1  General  Equations 


The  virtual  work  form  of  the  total  equilibrium  equations. 


. 
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was  expressed  in  Section  2.3,  Equations  2-17  and  2-18,  as 

6(I1  +  I2  -  I3)  =  0 


(4-1) 


where , 


I,  =  a.  .  u.  .  dV 

1  J 

v 


1 

2 


u 


ij  k,i 


u.  .  dV 

k,J 


(4-2a) 


(4-2b) 


I0  =  t.  u.  dS 

3  i  i 

* 

s 


(4-2c) 


The  displacement  at  any  point  of  the  thin-walled  beam  cross- 
section  may  be  considered  as  the  sum  of  beam  displacements  and  the 
individual  plate  deformations,  i.e. 

u.  =  u.B+u.P  (4-3a) 

l  ii 

Similarly  the  stresses  a.,  can  be  written  as  the  sum  of  the  stresses 

B  .  P 

due  to  beam  action  a..,  and  due  to  plate  action,  a... 

ij  ij 

a..  -  a..B+a..P  (4-3b) 

IJ  IJ  !J 

For  the  model  adopted,  the  stress  and  displacements  due  to  beam  action 
are  assumed  to  be  uniform  across  the  thickness  whereas  tne  stresses 
and  displacements  due  to  plate  action  are  proportional  to  the  distance 
from  the  middle  surface  of  the  plate  in  a  manner  consistent  with 
Kirchhoff’s  hypothesis. 

The  derivation  of  the  equilibrium  equations,  which  couples 
the  effects  of  cross— section  distortion  to  the  beam  equations,  may  now 
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be  carried  out  by  evaluating  the  expressions  in  Equations  4.2.  This 
is  done  below  in  Sections  4.3.2  to  4.3.4. 


4.3.2  Evaluation  of  I. 


The  integral  1^  of  Equation  4-1  may  be  evaluated  by  sub¬ 
stituting  Equations  4-3a  and  4-3b  into  Equation  4-2a. 


a .  . B  u .  B  dV  + 
ij 


P  P 

a. .  u.  .  dV 

ij  i,J 


(4-4a) 


=  IB+IP 
1  1 


(4-4b) 


P  P 

where  because  of  the  linear  variation  of  a. .  and  u 

1  >1 


B  P 
a.  .  u .  . 

ij 


P  B  n 
a. .  u.  .  =  0 

ij 


(4-4c) 


Using  the  definition  of  stress  resultants  from  Equations  2-1  these 


integrals  may  be  written  as 

rl 


B 


(Pw 1  -  M  v"  -  M  u"  +  W  4>"  +  T  (j>')  dz 
c  xs  ys  oo  sv 


(4-5a) 


n 

~  I 

±=1 


(HP  X 
a  a 


+  Xg  +  2MPg  XaB)i  dA 


(4-5b) 


where  A  denotes  the  middle  area  of  the  plate;  the  summation  is  carried 
out  for  all  plates;  and  xa»  Xag  denote  the  curvatures  of  the  middle 

surface  of  the  plate  in  the  cl  and  S  local  plate  co-ordinate  directions, 
shown  in  Fig.  4.1. 


4.3.3  Evaluation  of  ^ 


Integral  of  Equation  4-2b  may  be  evaluated  as  follows. 
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neglecting  the  effect  of 


where 


X2  =  X2  +  X2  +  I2BP 


B 


1 

2 


B  B  B 
a.  .  u.  .  u,  .  dV 
ij  k,i  k,j 


_1 

2 


B  p  P 

a .  .  u.  .  u.  .  dv 
ij  k,i  k,j 


BP 


B  B  P  AXT 

a.  .  u.  .  u.  .  dV 
ij  k,i  k,j 


(4-6) 


(4-7a) 


(4— 7b) 


(4-7c) 


and  the  symmetry  of  the  stress  tensor  has  been  used  in  expression  for 

I  BP 
2 

Using  the  notation  of  Chapter  II,  Equation  4-7a  becomes, 

l 


B 


1 

2 


J  o 


(-2V  d)u’  +  2V  cJ>v'+P(u'2  +  v,2  +  2e  u'  <j)f 
vy's  xTs  s  s  y  s 


-  2e  v'  <b')+  2M  v'  4> *  -  2M  u*  <J>  *  -I-  M  <J)'2)  dz 
xs  7  ys  xs  p 


(4-8) 


Considering  only  the  normal  stresses  in  the  z-direction;  retaining 
only  the  out  of  plane  plate  displacement  gradients  (see  Felippa  (19)) 
and  neglecting  the  effect  of  in-plane  shears;  Equation  4-7b  may  be 
written  as 


n 


I  G 


i=l 


j 

A 


0  B  t 
z 


{(uP  )2  +  (VP)2}  dA], 


(4-9) 


Similarly  Equation  4~7c  becomes 
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BP 
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l  [ 

i=l 
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B 

u 


>Z 


+  V* 


vB  )  dv]. 
z  1. 


(4-10) 


4.3.4  Evaluation  of 

If  surface  loads  are  regarded  as  line  loads  the  integral 
remains  as  in  Chapter  II,  and  is  expressed  as  Equation  A-48  of 
Appendix  A. 


4.3.5  Expression  for  total  virtual  work 

Substituting  the  expressions  for  1^,  I I^  into  Equation 
4-1  yields. 


<5[ 


(Pw*  -  M  v"  -  H  u"  +  W  <j>"  +  T  cj> ' )  dz 
v  c  x  s  y  s  w  sv 
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-  I 

i=l 


A 


(MPy  +  MPx0  +  2lr  x  „)  .  dA 

aAa  3  A3  a3Aap  i 


(4-11) 
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X  s 


n 


+  I  Cy  ^ZB  t{(uP  )2  +  (vp  )2}  dA  + 
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B  ,  p  B  p  B  .  -| 

a  (u  u  +  v  v  )  dV  J . 
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-  t .  u .  dS ]  =  0 

l  l 

s 

Equation  4-11  applies  to  any  thin-walled  beam  of  open  section,  composed 
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of  flat  plate  segments,  for  both  elastic  and  inelastic  response.  The 
terms  in  Equation  4--11  arising  from  Equation  4-5a,  ,  give  rise  to 

the  beam  flexure  tangent  stiffness  matrix,  for  inelastic  response,  and 
the  conventional  beam  flexure  stiffness  matrix  for  elastic  response. 
Those  arising  from  Equation  4-5b  give  rise  to  the  flexural  stiffness 
matrix  for  plate  bending.  Equations  4-7a  and  4-7b  lead  to  the  geo¬ 
metric  stiffness  matrix  for  the  work  done  by  the  initial  stress  state 
on  the  out  of  plane  beam  and  plate  deformations,  respectively,  whereas 
Equation  4-7c  represents  the  coupling  effect  between  plate  and  beam 
stresses  and  deformations. 

4 . 4  Specialization  for  Elastic  WF  Sections 

For  elastic  response  the  conventional  beam  stress  resultants 
appearing  in  Equation  4-11  are  defined  in  Chapter  II.  The  stress 
resultants,  M.p,  Mp,  MP^  for  the  plate  segments  may  be  expressed  as 
(see  Fig.  4.1) 


MP 

a 

=  -  D 

(xa  +  V 

Xg) 

(4-12a) 

Mp 

=  -  D 

(Xg  +  V 

(4-12b) 

Mae 

=  -  D 

a  -  v) 

Xag 

(4-12c) 

where  v  is  poisson' 

s  ration 

and  D  = 

=  E  t3  /  12  (1  -  v2) 

(4-13) 

Assuming  a  symmetric  WF  shape  and  designating  the  out  of  plane  plate 
segment  displacements  by  v^P ,  v^  and  u^P  as  shown  in  Fig.  4.2; 
expressing  the  curvatures  in  Equation  4-11  in  terms  of  the  derivatives 
of  these  displacements;  expressing  the  stress  resultants,  in  the  first 
set  of  braces  in  Equation  4-11,  in  terms  of  displacements;  and. 
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jg 

approximating  a  by  the  normal  beam  equation  (the  first  three  terms 

z 


of  Equation  A-51) ;  the  component  expressions  of  Equation  4-11  becomes 

cl 

•2 


uq8)  =  ql 
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(4-14e) 


Where  the  local  coordinate  system  is  shown  in  Fig.  4.2;  u  and  v  refer 
to  X  and  Y  coordinate  directions;  the  subscripts  b,  t  and  w  refer  to 
’bottom  flange',  'top  flange’  and  ’web’  respectively;  the  superscript 
P  indicates  plate  bending  displacement;  and,  u^ ,  v  ,  <J>  and  w  refer  to 

o  S  C 


beam  bending  displacements. 


4 . 5  Finite  Element  Model 

The  plate  segment  displacements  for  a  wide  flange  beam  are 
assumed  to  be  represented  by  the  shape  functions  shown  in  Fig.  4.3. 
These  functions  imply  that  on  any  section  perpendicular  to  the  longi¬ 
tudinal  axis  the  flange  remains  straight.  With  this  model  there  are 
four  degrees  of  freedom  associated  with  cross  section  distortion  at 
each  node,  in  addition  to  the  seven  degrees  of  freedom  associated  with 
beam  action  at  each  node.  The  sketches  in  Fig.  4.3  show  the  deflected 
shape  of  the  I-section  due  to  unit  values  of  the  generalized  displace¬ 
ment  coordinates  corresponding  to  each  degree  of  freedom  of  the  top 
flange  while  the  others  are  kept  zero.  Shape  functions  for  the  bottom 
flange  are  similar  to  those  for  the  top  flange. 

The  plate  displacements  v^  can  be  written  in  terms  of  the 


nodal  displacements  as 
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(4-15) 


where  £  and  c  and  n  are  the  dimensionless  parameters  illustrated  in 
Fig.  4.2,  and  expressed  as 


5  =  2x/bf ;  C  =  zll;  n  =  2y/dw 


(4-16) 


and  (v^)  ,  £(v^  )  ,  (v^)9,  l(vP  )_  are  the  nodal  values  of  the 

generalized  coordinates  as  shown  in  Fig.  4.3. 

Denoting  the  set  of  shape  functions  by,  <k  >  the  displacements 
in  the  top  and  bottom  flanges  are  expressed  as 
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Where  and  represent  the  vectors  of  generalized  nodal  coordinates. 

Since  the  angle  at  the  junction  of  the  web  and  the  flange 
remains  unchanged 

TP 


9ti=  hi  ■  2'bf{KW 


(4-18) 


Where  0™  is  the  rotation  at  the  top  of  the  web,  and  0  is  the  top 

tl 

flange  rotation  at  node,  1.  Using  a  similar  notation  at  the  bottom 
junction  we  may  write  these  rotations  along  the  element  length  as 
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Where  anti-clockwise  rotation  is  assumed  as  positive  for  the  web,  b^ 

is  the  flange  width  and  d  is  the  depth  of  the  web. 

°  w 

Choosing  the  local  coordinate  system  at  the  centre  of  the 
web,  as  shown  in  Fig.  4.2,  the  web  displacements  are  expressed  as 
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and  using  Equation  4-19 ,  Equation  4-20  may  be  expressed  in  terms  of 
the  generalized  coordinates,  as 
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Substituting  the  distributed  displacements,  in  terms  of  nodal  coor¬ 
dinates,  into  Equations  4-14  and  integrating,  we  may  write  the  element 
equilibrium  relationship  between  the  generalized  forces  and  generalized 
displacements  in  the  symbolic  form 


[[ks]  +  [kG]]  {rE>  =  {ReJ  (4-22) 

where  the  matrices  are  derived  in  Appendix  C  and  have  the  following 


structure. 
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and 


guu 

• 

gu* 

gu  l 

gvv 

gvl 

[kG]  - 

g(j)U 

• 

g<}>v 

» 

• 

g<j>£ 

• 

glu 

glv 

gl<p 

giz 

(4-23b) 


<r  >  =  <<u>  <V>  <_£>  <W>  <£•>> 

E 


(4~23c) 


The  nodal  vectors  u,  v,  and  w  have  been  defined  in  Chapter  II  and 
nodal  vector  l  corresponds  to  the  eight  degrees  of  freedom  associated 
with  local  buckling  i.e.: 


UY 


-  <vlv  < 2 ’  v^2> 


(4-24) 


Upon  assembly,  the  equilibrium  equations  for  the  whole  structure  may 
be  written  as , 


[[Kg]  +  [Kg]]  {r}  =  {R} 


(4-25) 
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The  buckling  problem  arises  from  the  variation  of  the  homogeneous  matrix 
equation  in  the  form 

[[Kg]  +  [kg]]  {Ar>  =  0  (4-26a) 

or 

([Ks]  +  X[KG1)  {Ar}  -  0  (4-26b) 

Solving  the  above  eigenvalue  equation  for  the  unknown  value 
X  and  multiplying  it  by  the  nominal  stress  distribution  yields  the 
critical  stress. 

4 . 6  Illustrative  Solutions 
E-l  Local  Buckling  of  Plates 

To  test  the  plate  buckling  representation  a  number  of  pro¬ 
blems  were  solved  without  coupling  the  local  buckling  to  member  buck¬ 
ling.  In  Table  4.1  the  solutions  for  buckling  of  plates  for  various 
boundary  conditions  are  compared  with  the  closed  form  solutions.  The 
nomenclature  of  Table  4.1  is  indicated  in  Fig.  4.4.  The  Duckling  of 
a  single  plate  can  be  simulated  by  setting  the  thickness  of  the  other 
plates  to  zero.  The  model  yields  good  results  for  flange  plate  buck¬ 
ling  but  is  too  stiff  to  yield  accurate  results  for  web  plate  buckling 
since  only  a  linear  variation  curvature  can  be  represented  across  the 
plate. 

E - 2  Uncoupled  Local  Buckling  of  Columns  subjected  to  an  Axial  Load 
Fig.  4.5  shows  the  results  obtained  when  an  axial  load  is 
applied  to  a  symmetric  wide— flange  column  of  constant  thickness  and 
the  critical  local  buckling  load  (uncoupled  solution)  is  calculated  for 


. 


. 
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various  ratios  of  br/d  .  The  minimum  plate  buckling  coefficient,  k  , 

f  w  w 

was  computed  from  the  critical  stress,  and  the  relationship  of  k  to 

w 

b,/d  is  compared  with  Bulson's  (15)  theoretical  curve  in  Fig.  4.5. 
f  w 

It  can  be  seen  that  the  model  closely  represents  the  flange  buckling 
of  a  wide  flange  beam,  but,  as  is  to  be  expected  does  not  accurately 
reproduce  web  buckling  behaviour  because  the  assumed  shape  functions 
cannot  represent  the  fundamental  mode  of  this  plate. 

E-3  Interaction  of  Local  and  Member  Buckling 

Fig.  4.6  shows  the  buckling  load  for  an  axially  loaded  wide 
flange  column  for  (a)  uncoupled  local  buckling,  (b)  uncoupled  Euler 
buckling,  (c)  interactive  column  and  local  buckling.  The  interaction 
curve  gives  values  for  the  buckling  coefficient  only  a  few  percent 
lower  than  the  Euler  column  buckling  curve.  These  results  agree  with 
those  of  Pfluger  and  Bulson  (16) .  Curves  are  also  shown  for  the  case 
where  the  load  is  acting  at  an  eccentricity  of  5". 

E-4  Interaction  of  Local  and  Lateral  Buckling 

Fig.  4.7  shows  a  comparison  of  the  (uncoupled)  elastic  local 
buckling  load  and  the  (uncoupled)  lateral  buckling  load,  with  the 
critical  load  including  interactive  effects,  for  a  simply  supported 
beam. 

4 . 7  Summary 

The  formulation  developed  in  this  chapter  can  form  the  basis 
for  elastic  or  inelastic  buckling  of  a  thin— walled  beam  with  local  and 
member  buckling  interaction.  The  model  adopted  in  this  chapter  for 
elastic  wide  flange  beams  has  eleven  degrees  of  freedom  at  each  node 
compared  to  eighteen  degrees  of  freedom  for  the  simplest  finite  element 
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plate  representation,  illustrated  in  Fig.  4.8. 

The  model  developed  in  this  chapter  appears  to  yield  good 
results  for  local  flange  buckling  but  should  not  be  expected  to 
represent  web  buckling.  The  results  verify  the  validity  of  the  class¬ 
ical  approach  to  local  buckling  by  showing  that,  for  wide  flange 
sections  in  the  elastic  range,  the  interactive  effect  is  normally 
small,  and  local  and  member  buckling  are  initiated  as  relatively 
detached  phenomenon  over  most  of  the  range  of  loading. 


. 


TABLE  4-1  SOLUTIONS  FOR  UNCOUPLED  LOCAL  BUCKLING  (FIG.  4.4) 
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FIGURE  4.1  LOCAL  PLATE  CO-ORDINATE  DIRECTIONS 
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FIGURE  4.2  LOCAL  COORDINATE  SYSTEM  FOR  PLATE  DEFORMATIONS 
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FIGURE  4.3  DEFLECTED  SHAPE  OF  THE  I-SECTION  DUE  TO  UNIT  VALUE 
OF  THE  GENERALIZED  DISPLACEMENT 
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n  -  number  of  elements 


FIGURE  4.4  GENERAL  NOTATIONS 


PLATE  BUCKLING  COEFFiClENl  'K 
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FIGURE  4.5  LOCAL  BUCKLING  OF  COLUMNS  SUBJECTED  TO  AXIAL  LOAD 


«* 


91 


UJ 


32  g 

o  ° 

H3  V"' 
m  ° 

U.J  fQ 

_  i— 

ro  cd 
o  -w 
o  c; 


A 


in 


o 


ii 

~e- 


n 


co 


o 


CM 


■J  o 
o 


0  iA/l4  l'l^=l)j 
Z  Z 


FIGURE  4.6  COLUMN/LOCAL  BUCKLING  INTERACTION 
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FIGURE  4.7  INTERACTION  OF  LATERAL  AND  LOCAL  BUCKLING 
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FIGURE  4.8 


PLATE  ELEMENTS  WITH  18  DEGREES  OF  FREEDOM  AT 
EACH  NODE  OF  THE  BEAM  ELEMENT 


CHAPTER  V 


FORMULATION  FOR  INELASTIC  MEMBERS 

5 . 1  Introduction 

Nonlinear  formulations  relating  kinematic  and  mechanical 
variables  may  fall  into  the  three  following  categories:  1)  nonlinear 
constitutive  equations  (material  nonlinearity) ,  2)  geometric  nonlinear¬ 
ity,  3)  combined  material  and  geometric  nonlinearity.  The  beam-column 
problems  considered  in  Chapter  III  fall  into  the  second  category .  This 
and  the  following  chapter  are  devoted  to  beam-column  problems  falling 
into  the  third  category.  A  comprehensive  review  of  the  application  of 
finite  element  concepts  to  nonlinear  problems  in  structural  mechanics 
has  been  given  by  Oden  (56) . 

Basically  there  are  three  types  of  finite  element  solution 
methods  which  can  be  adopted  for  nonlinear  problems  (30,34);  1)  incre¬ 
mental  methods  without  equilibrium  checks ;  2)  incremental  methods  with 
equilibrium  checks;  3)  direct  solution  of  the  governing  nonlinear- 
equations.  The  first  type  of  approach  is  computationally  fast  but  has 
the  disadvantage  that  equilibrium  is  not  satisfied  at  each  step  and 
there  is  no  direct  indication  of  the  errors  involved  in  the  solution. 

In  the  second  type,  incremental  solutions  due  to  a  series  of  load 
steps  are  determined,  equilibrium  is  checked,  and  corrections  are  made 
for  each  load  step,  if  necessary.  Unlike  the  first  method  this 
insures  that  the  equilibrium  equations  are  satisfied  throughout  the 
loading  history.  The  beam-column  solutions  presented  here  have  been 
carried  out  using  this  method.  The  third  method  has  been  successfully 
applied  by  Oden  (55) . 
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5 . 2  Formulation  of  Inelastic  Incremental  Equilibrium  Equations 

5.2.1  General  Description  of  Solution  Procedure 

Details  of  the  solution  procedures  for  inelastic  small 
deflection,  beam-column  and  bifurcation  problems  are  discussed  in 
Chapter  VI.  However,  since  the  formulation  is,  to  some  extent, 
dependent  on  the  solution  procedure  adopted,  a  general  description 
of  this  procedure  for  beam-column  problems  is  given  in  this  section. 

The  solution  procedure  for  inelastic  beam-columns  employed 
in  this  study  may  be  described  as  follows.  In  the  formulation, 
strains,  displacements  and  rotations  are  assumed  to  be  ’small’  and 
hence  a  second  order  small  displacement  theory  is  used.  A  linear 
incremental  solution  technique  is  adopted  with  equilibrium  checks. 
This  is  similar  to  the  method  adopted  by  Murray  (50)  and  Wilson 
(51,52),  Sharif i  (63),  and  Hofmeister  and  Greenbaum  (34). 

The  incremental  equilibrium  Equation  2-35,  was  derived  in 
Chapter  II  for  elastic  response.  A  similar  equation  is  derived  in 
this  chapter  and  may  be  expressed  as 


[[Kg]  +  [Kg]]  {Ar}  =  {AR} 


(5-1) 


where  [k  ]  is  the  tangent  stiffness  matrix  for  the  increment  of  loading. 

The  Newton-Raphson  method  for  inelastic  response  is  explained 
in  the  following  algorithm. 

1)  For  any  approximate  {r}^,  tlie  total  strains  may  be  determined  from 
the  kinematic  relationships.  Equation  A-36 .  The  stresses  may  then  be 
computed  directly  from  the  stress-strain  relationship,  shown  in  Appendix 


D  in  Fig.  D.la. 
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2)  Knowing  the  strains  and  stresses,  the  tangent  stiffness  matrix, 
the  geometric  stiffness  matrix,  and  the  unbalanced  forces  may  be 
evaluated,  as  detailed  in  Appendix  D.  Equation  5-1  may  be  solved  to 
determine  the  displacement  increments  to  equilibrate  the  unbalanced 
forces ,  i . e . : 

{Ar}  =  [[kJ]  +  [K  ]  ]_1{AR}  (5-2) 

n  bn  bn  n 

3)  The  displacements  are  updated 


(r> 


n+1 


{r}  +  {Ar} 

n  n 


(5-3) 


4)  Steps  1  to  3  are  repeated  until  the  unbalanced  force  vector  is 
arbitrarily  small. 

5)  A  new  load  increment  is  applied  and  steps  1  to  4  are  repeated. 

This  procedure  is  identical  to  that  of  Section  3.4.2  except 
for  the  detailed  evaluation  of  [K  ]  and  [K  ] .  However ,  a  modified 

b  vj 

Newton-Raphson  technique,  holding  the  tangent  stiffness  constant  for 
several  iterations,  was  actually  used  for  the  numerical  work. 

The  procedure  of  computing  the  error  vector  in  nodal  equil¬ 
ibrium  is  described  as  follows.  The  total  forces  at  a  node  are  the 
sum  of  a)  the  applied  forces,  b)  the  resisting  forces  due  to  stresses, 
which  may  be  computed  from  statics  for  the  stress  resultants  produced 
by  any  stress  distribution,  and  c)  the  forces  resulting  from  changes 
in  geometry  which  may  be  evaluated  as  [K^]{r}  for  any  set  of  stress 
resultants  and  displacements.  The  stress  resultants  are  computed  at 
the  ends  of  each  element  by  a  numerical  integration  of  stresses  which 
have  been  previously  determined  from  the  element  strains  as  detailed  in 


. 
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Appendix  D.  The  difference  between  the  applied  loads  and  the  sum  of 
a)  the  resisting  forces  and  b)  the  ’geometric*  forces,  represents  the 
set  of  unbalanced  forces,  designated  as  {AR}  in  Equation  5-1,  on 
the  given  configuration  of  the  model. 

5.2.2  Other  Formulations 

Other  methods,  such  as  the  initial  strain  method  (38),  and 
the  initial  stress  method  (77)  ,  have  been  adopted  for  the  analysis 
of  elastic-plastic  problems.  In  the  initial  strain  method,  fictitious 
loads  are  applied  to  the  structure  to  produce  elastic  strains  numeri¬ 
cally  equal  to  the  plastic  strains.  Although  this  requires  only  one 
inversion  of  the  stiffness  matrix,  an  iterative  solution  for  the 
fictitious  forces  must  be  carried  out  in  order  to  control  solution 
errors.  The  initial  stress  approach,  developed  by  Zienkiewicz,  et.  al. 
(77),  uses  a  slightly  different  incremental  technique,  treating  dis¬ 
crepancies  between  elastic  and  elastic-plastic  stress  increments  as 
unbalanced  body  forces  and  iterating  until  these  forces  disappear. 

The  main  advantage  of  these  methods  is  that  the  assembled  stiffness 
matrix  remains  unchanged. 

5.2.3  Assumptions  for  Derivation  of  Tangent  Stiffness 

Xn  deriving  the  tangent  stiffness  for  a  partially  plastic 
beam  section  it  is  assumed  that  no  strain  reversal  occurs.  If  one  con¬ 
siders  a  section,  which  has  already  yielded,  and  an  infinitesimal 
increment  in  bending  moment  is  applied,  the  resistance  to  this  moment 
is  a  measure  of  incremental  bending  stiffness.  The  yielded  zones, 
will  be  increased  somewhat,  while  some  strain  reversal  may  take  place. 
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Pn  the  already  yielded  zones.  These  effects  are  neglected  in  this 
formulation  and  the  additional  moment  is  assumed  to  be  resisted  by 
all  parts  of  the  cross  section  with  the  stress  increment  determined 
from  the  tangent  modulus  relationship  for  the  material  at  the  current 
value  of  strain.  The  section  properties  relating  increments  in  stress 
resultants  to  increments  in  displacements  may  then  be  computed  by  a 
transformed  area  concept.  This  concept  may  be  utilized  for  thin-wall 
open  beam  sections  since  the  stress  is  essentially  uniaxial  and  the 
effect  of  shear  stress  on  the  principal  stresses  may  be  neglected. 


5.2.4  Basis  for  Incremental  Equilibrium  Equations 

Since  the  incremental  equilibrium  equations  derived  in 
Chapter  II  are  virtual  work  equations,  they  may  be  used  directly  for 
the  incremental  formulation  of  inelastic  problems.  Equation  2  14  is 
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t.  <5u.ds 
l  l 


(5-4) 


The  physical  interpretation  of  the  terms  in  this  equation 
was  discussed  in  Section  2.3  and  the  initial  displacement  terms  were 
shown  in  Appendix  B  to  have  negligible  effect  for  this  class  of  pro¬ 
blem  in  the  range  of  interest  under  consideration.  They  are  therefore 

neglected  in  the  following  formulation. 

An  examination  of  the  first  term  of  Equation  5-4,  from  which 
the  geometric  stiffness  arises,  indicates  that  it  is  dependent  only  on 
the  initial  stress  and  the  geometry  of  the  deformation  and  therefore 


. 


•  *  * 
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is  independent  of  the  constitutive  relationship.  Hence  the  geometric 
stiffness  terms  derived  in  Chapter  II  remain  unchanged  and  can  be  used 
directly  in  the  solution  of  problems  involving  nonlinear  material 
behaviour . 

For  inelastic  material  response,  the  location  of  the  instant¬ 
aneous  centroid,  shear  centre  and  sectorial  centroid,  and  the  incre¬ 
mental  principal  axes  are  functions  of  the  current  strain  distribution. 
For  elastic  response  it  was  convenient  to  refer  displacement  increments 
to  the  original  shear  centre  axis.  However  for  inelastic  response, 
the  variation  of  increments  in  displacements,  in  Equation  5-4  may  be 
referred  to  either 

a)  the  original  shear  centre,  centroid,  sectorial  centroid,  and 
reference  axes  of  the  elastic  section  (reference  points  S,  C  and  B,  and 


coordinates  x-y  of  Fig.  5.2) 


or 


b)  the  instantaneous  shear  centre,  instantaneous  centroid  and  instant¬ 
aneous  sectorial  centroid  of  the  transformed  section,  in  the  instantan¬ 
eous  principal  directions  (reference  points  S,  C  and  B'  and  coordinates 
£-n  of  Fig.  5.2.) 

The  element  equilibrium  equations  will  be  derived  with  respect 
to  both  of  these  reference  systems.  The  formulation  with  respect  to 
the  original  reference  axes  of  the  element  (the  x-y  coordinates  of 
Fig.  5.2)  will  be  referred  to  as  'Formulation  1'.  The  formulation  with 
respect  to  the  transformed  section  principal  axes  (the  coordinates 

of  Fig.  5.2)  will  be  referred  to  as  'Formulation  2'.  The  transforma¬ 
tion  of  displacements  to  the  global  coordinate  system  (reference  points 
s'  and  C  ,  and  coordinates  X,Y)  will  be  discussed  in  Section  5.4,  and 


. 
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becomes  a  part  of  the  assembly  process.  It  is  shown  in  Appendix  E  that 
these  two  formulations  are  substantially  the  same. 


5.2.5  Formulation  of  Incremental  Equilibrium  Equations  -  ’Formulation  1* 
Neglecting  the  initial  displacement  matrix,  and  taking  the 
variation  with  respect  to  the  original  reference  system.  Equation  5-4 
may  be  written  for  one  element  as 


{AP  6w’  -  AM  6v"  -  AM  6u"  +  AW,  6(f)"  +  V  <f>  Sv' 

cxsys^  xs 

l 

+  V  v1  66—  V  6  6u'  -  V  u!  66  +  P  u’  5u'  +  P  v'  6v’ 
xs  y  sys  ss  ss 


+  P  e  u '  66  ’ 

y  s 


+  P  e  6*  6u'  -  Pe  v’  66’  -  Pn  T  6v' 
y  Y  s  xs  x  s 


-  M  u'  56’  +  M  v’  66’  “  M  6*  5u'  +  M  (j)  ’  6v'  +  M  6* 
x  s  ^  y  s  X  sy  sp 

+  T  66’  -  (Aq  6u  +  Aq  ov  +  Aq  6w  +  Am.  6 <J>  +  Am  6v’ 
sv  y  Hx  S  y  S  Z  Z  t  X  s 

+  Am  6u')}dz  =  0  (5-5) 

y  s 

where  P  M  M  ,  M  ,  V  ,  V  are  the  total  stress  resultants  existing 
5  x 3  y  p  x’  y 

prior  to  the  load  increment;  and  AP ,  AM^,  AM  and  AW^  are  increments 
in  stress  resultants,  with  respect  to  the  x— y  reference  axes  travel¬ 
ling  along  with  the  section,  as  shown  in  Fig.  5.2.  The  total  stress 
resultants  may  be  evaluated  by  direct  integration  of  stresses  as 
detailed  in  Appendix  D. 

If  each  element  of  area  is  'transformed1  such  that  the  pro¬ 
duct  of  the  current  tangent  modulus  times  the  original  element  of  area 
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is  equal  to  the  original  modulus  times  the  element  of  transformed  area; 

the  stress  resultant  increments,  AP,  AM  ,  AM  and  AW  appearing  in 

x  y  m 

Equation  5-5  are  given  by  Equation  A-26,  as 


AP  =  E  AT  w’  -  E  AT  y  vn  -  E  AT  x  u"  +  E  ST  6" 
C  OS  OS  CJ 


(5-6a) 


AM 


E  ATy  w’  -  E  IT  v"  -  E  IT  u"  +  E  ST  6" 
o  c  x  s  xy  s  wx 


(5-6b) 


AM 


E  ATx  w'  -  E  IT  v"  -  E  IT  u"  +  E  ST  6" 
o  c  xy  s  y  s  wy 


(5-6c) 


AW 


03 


E  ST  w'  -EST  v"  -  E  ST  u"  +  E  IT  6" 
a)  c  tox  s  (uy  s  U) 


(5— 6d) 


T  T  T  T  T  T  T  T 

where  A, I, I  ,  I  ,  S  ,  S  ,S  ,1  represent  the  properties  of  the 
x  xy  y  co  ux  my  u> 

transformed  section  with  respect  to  points  C  and  S  and  the  x-y  reference 
axes.  The  distances  x^,  y^,  of  Equation  5-6,  and  e^,  e^  of  Equation 
5-5  refer  to  the  location  of  C  and  S,  as  shown  in  Fig.  5.3. 


5.2.6  Formulation  of  Incremental  Equilibrium  Equations  ’Formulation  2* 
Taking  the  variation  of  displacement  increments  with  respect 
to  the  instantaneous  (transformed)  principal  axis  reference  system. 
Equation  5-4  may  be  written  for  one  element  as 


AP  6w  -  AM  Svm  -  AM  6u"  +  AW  6 "  +  V  v’  56  -1  V  cf>  6v* 

5  bn  n  b  bn  bn 


—  V  d>  6ut  -  V  uT  66  +  Pu’  6u'  +  Pv  6v  4-  Pe  u  66  +  Pe  6  5ur 

n^bnb  bb  nn  nb  n  b 

-Pe  v'  66'  -  Pe_6 1 5v ’  -  Mru!  66’  +  M  v*  66’  -  H  6*  6u*  +  M  6*  5v* 

bn  b  nbb  nn  b  bn  n 


+  M^r6'  56'  +  T  66 '  -  (Aq  6u_  +  Aqn  6vn  4-  Aqr  6w^  +  Am  66 

C'ir  b  b  r* 


Pb 


sv 


n  n 


lb  b 
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+  Am_  Sv*  +  Am  6u') 

5  n  n  5 


dz  =  0 


(5-7) 


where  E,  and  n  are  the  transformed  section  principal  axes  of  Fig.  5.2; 
the  displacement  increments  and  coordinate  relationships  are  shown  in 
Fig.  5.3;  and  the  stress  resultant  increments  are  shown  in  Fig.  5.4. 

The  stress  resultant  increments,  AP ,  AM^ ,  AM^,  and  AW^ 
appearing  in  Equation  5-7  may  be  written  in  the  following  uncoupled 
form  (see  Appendix  A) 

AP  =  E  AT  w^  (5 -8a) 

C 

AM  =  -E  if  v"  (5-8b) 

€  5  n 

AM  =  -E  IT  u"  (5-8c) 

n  n  5 

AW  =  E  IT  <J>"  (5-8d) 

T  T  T 

where  A  is  the  transformed  area;  and  1^,  are  the  principal  moments 

T 

of  inertia  of  the  transformed  section  about  the  £-r)  axes;  I  is  the 
principal  sectorial  moment  of  inertia  of  the  transformed  section  about 
s’*  M  and  M  are  the  stress  resultants  about  the  principal  axes;  and 

5  n 

M  is  the  stress  resultant  defined  in  Equation  2-7,  about  S. 

P4 

5 . 3  Finite  Element  Idealization 

The  finite  element  equations  for  one  element  may  be  obtained 
from  Equations  5—5  and  5—6,  or  Equations  5—7  and  5-8,  by  substituting 
for  the  displacements  in  terms  of  shape  functions  and  nodal  parameters, 
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and  carrying  out  the  required  integration.  The  total  virtual  work  may 
be  evaluated  by  summing  the  integrals  over  each  element,  and  the 
reference  axis  in  each  element  may  be  different  providing  the  displace¬ 
ment  variations  are  compatible  at  the  junction  of  any  two  elements. 

For  sections  with  variable  properties,  a  linear  variation  of 
axial  displacement  is  incapable  of  representing  a  constant  axial  load. 
Therefore,  an  additional  degree  of  freedom  for  axial  displacement  is 
introduced  at  the  midpoint  of  an  element  to  achieve  a  balance  of  axial 
forces  at  the  nodes.  Hence  the  quadratic  displacement  function  of 
Equation  2-23e  is  used  to  represent  the  axial  displacement. 

To  determine  the  element  stiffness  for  'Formulation  1* , 
Equations  5—6  are  substituted  into  Equation  5—5,  and  Equations  2— 23a, 
b,  c  and  e  are  then  substituted  for  displacements.  The  element  tan¬ 
gent  stiffness  matrix  that  arises  in  this  substitution  is  derived  in 
Aopendix  E.  The  geometric  stiffness  that  arises  is  also  discussed 
in  Appendix  E  and  differs  from  that  in  Equation  2-28  only  in  the 
[g  ]  submatrix.  The  element  equilibrium  equations  then  become 

M 

[kg]  { ArE>  +  [kG]  UrE>  =  {AR£}  (5-9) 

T 

where  (Ar  }  is  defined  in  Fig.  2.6;  [k  ]  is  now  a  full  matrix  with 
E  b 


submatrices  denoted  as 


. 

. 
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(Equation 
2-28,  but 


E-5  of  Appendix  E) ;  and  [k  ]  is  of  the  same  form  as  Equation 
[g  ]  is  given  by  Equation  E-7a.  In  the  derivation  of  Equa- 

M 


tions  5-9,  and  5-10,  a  linear  variation  of  stress  resultants  and  section 
properties  has  been  assumed. 

To  determine  element  stiffnesses  for  ’Formulation  2',  Equa¬ 
tions  2-23a,  b,  c  and  e  are  assumed  to  represent  displacement  incre¬ 
ments  referred  to  the  S  and  C  axes,  in  the  E, ,  n  coordinate  system. 

With  this  assumption,  the  evaluation  of  stiffness  matrices  based  on 
Equations  5—7  and  5—8  may  be  carried  out.  This  results  in  a  set  of 
element  equilibrium  equations  of  the  same  form  as  Equation  5-9. 

However ,  £k  J  is  now  given  by  Equation  E— 10  and,  as  in  the  elastic 
case,  is  uncoupled.  The  form  of  is  identical  to  that  in  Formu 

lation  1"  but,  because  of  the  different  reference  axes,  the  stress 
resultants  have  different  magnitudes.  The  displacement  increment 


vector 


{ArE} 


is  shown  in  Fig. 


5.5. 


5 . 4  Transformation  and  Assembly 

5.4.1  Reference  Configurations 

The  element  equilibrium  Equation  5-9  has  been  derived  for 
two  different  local  coordinate  systems.  In  general  the  location  and 
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orientation  of  the  local  reference  axes  may  vary  from  element  to 
element.  Prior  to  assembly  it  is  therefore  necessary  to  transform 
the  nodal  displacements  and  forces  from  the  local  coordinate  system 
to  a  global  coordinate  system  in  order  to  maintain  compatibility  of 
deformations . 

For  inelastic  analysis  it  is  assumed  that  the  global  refer- 
ence  points  C  and  S  (Figs.  5.1  and  5.3),  and  the  X-Y  coordinate 
system,  are  the  original  elastic  centroid,  shear  center,  and  centroidal 
reference  axes,  respectively •  An  arbitrary  section  after  deformation  is 
shown  at  P ,  in  Fig.  5.1.  This  same  section  is  shown  at  different 
stages  of  the  deformation  in  Fig.  5.2.  The  original  position  is  indi¬ 
cated  as  position  I.  Application  of  the  load  causes  the  section  to  be 
displaced  to  position  II.  The  subsequent  addition  of  a  small  load 
increment  causes  the  section  to  be  further  displaced  to  position  IV. 

The  movement  from  position  II  to  position  IV  results  from  a 
simultaneous  translation  and  rotation  of  the  cross  section.  For  con¬ 
venience  the  movement  caused  by  the  increment  of  load  will  be  considered 
to  occur  in  two  parts.  The  first  part  is  the  translational  displace¬ 
ment  from  position  II  to  III.  The  second  part  is  the  movement  from 
position  III  to  IV,  caused  by  twist. 

5.4.2  Displacement  Transformation 

"Formulation  l11 

The  nodal  vector  of  Equation  5-9  for  'Formulation  1'  is 
expressed  in  terms  of  displacement  increments  of  points  C  and  S  of  posi¬ 
tion  II  (Fig.  5.2)  in  the  x  and  y  coordinate  directions  which  are  the 

axes  travelling  with  the  cross  section.  ihese  increments  may 
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be  related  to  the  displacement  increments  of  C  and  S  ,  for  a  typical 
node  by  the  transformation  (Fig.  5.3). 


r%~\ 

cos  $  ,  sin  $ 

rAU  A 

s 

V 

s 

-sin  $  ,  cos  $ 

AV 

s 

w 

c 

1 

AW 

c 

♦ 

= 

1 

< 

A$ 

u'l 

s  1 

Jt  cos  $  ,  t  sin  $ 

AU' 

s 

v*£ 

s 

-l  sin  $  ,  t  cos  $ 

AV' 

s 

l 

A$ ' 

\  y 

v  / 

For  nodes  p  and  q  at  the  ends  of  an  element ,  Equation  5  11  may  be 
written  symbollically  as 

{up}  =  [T  ]  {AUP}  (5-12a) 

-s  p  — 

(Uq)  =  [T  ]  { AUq}  (5-12b) 

— s  q  — 

where  the  transformation  matrices  of  Equations  5-12  are  defined  by 
matrix  of  Equation  5—11,  with  $  taking  the  values  of  $  and  $  in 
[I  1  and  [T  1,  respectively.  Reordering  the  vector  of  element  nodal 

p  q 

displacements,  it  may  be  expressed  in  the  local  and  global  coordinate 
systems,  as 


=  <u?  ,  vh  vP,  <fp ,  l  Ug  .  £  'A  1  *P7.'  wcm> 


U^,  Vq,  Wq,  4>q  J  £  Ug  Z>  t  VS.7.»  ^  ^  > 


s  s 


s ,  z 


(5-13a) 


ID 


107 


and 


fArE}G 


=  <AUP,  AVP ,  AWP,  A$P,  AUP  AVP  ^ ,  A$P^,  AW^, 


s  ,z  s  ,z 


AUq,  AVq,  AWq,  A$q ,  AUq  ,  AVq  ,  A$q  > 
s’  s’  c  *  s,z’  s,z’  ,z 


(5-13b) 


which  are  related,  through  Equations  5-12,  in  the  expression 


^ ArE^L 


T 

4upN 

P 

. 

1 

AW 

1 

— m 

T 

AUq  , 

/ 

q 

(5-14a) 


Symbolically  we  may  write  Equation  5-14a  as 


{ Ar„}T  =  [Tj  {Ar  } 


EL 


EG 


(5-14b) 


5.4.3  Displacement  Transformation  -  "Formulation  2" 

The  nodal  vector  of  Equation  5-9,  for  'Formulation  2',  is 

expressed  in  terms  of  the  displacement  increments  of  points  C  and  S 

of  position  II  of  Fig.  5.2,  in  the  5  and  q  coordinate  directions, 

(which  represent  the  centroid  and  shear  center  of  the  displaced  trans 

formed  section,  and  the  principal  axes  of  tne  displaced  transformed 

section,  respectively).  Let  a  be  the  angle  between  the  instantaneous 

5,  n  principal  axes  and  the  original  x,  y  reference  axes  travelling 

th.  the  section.  The  relationship  between  the  displacement  incre 

ments  of  C*  and  S*,  defined  in  Section  5.4.2,  and  those  in  the  vector 

(Ar  }  of  Equation  5-9,  can  now  be  derived  as  follows 
E 
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From  Fig.  5.2  and  Fig.  5.3 


A  A 

u  =  u  cos  a  +  v  sin  a 


A  A 

v  =  -u  sin  a  +  v  cos  a 

n 


(5-15a) 

(5-15b) 


where  u  ,  v  are  the  displacement  increments  of  the  instantaneous 

5  n 


shear  center  along  the  principal  axes  directions,  and  u,  v  represent 
the  displacement  increments  of  the  instantaneous  shear  center  in  the 


directions  parallel  to  the  x,  y  reference  axes. 

Expressing  u,  v  in  terms  of  u  ,  v  (the  displacements  of  the 

s  s 

original  shear  center  in  the  direction  of  the  reference  axes)  results 
in 


v  =  v  -t-(b  —  e  )  4> 

s  x  x 


u  =  u  -  (b  —  e  )  4> 

s  y  y 


(5-16a) 

(5-16b) 


Substituting  Equation  5-16  into  Equation  5-15  yields 


u  =  u  cos  a  +  v  sin  a  +  ( (b  -  e  )  sin  a  -  (b  ~  e  )  cos  a)4> 

^s  s  xx  y  y 

v  =  -u  sin  a  +  v  cos  a  +  ( (b  -  e  )  sin  a  +  (b  -  e  )cos  ot) 4>  (5-17) 

ns  s  y  y  x  x 


where  (j)  is  the  increment  in 

For  the  small  displacement  problem  it  can  be  assumed  that 


<j)  =  $  (5-18) 

Now  u  and  v  referenced  to  the  configuration  at  position  II  should  be 
s  s 

transformed  to  the  global  system.  If  AU  ,  AV  are  the  increments  of 
displacements  of  the  original  shear  center,  S  ,  in  the  global  system. 


then , 


■ 


. 
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AU  cos  <£>  +  AV  sin  $  =  u 

S  S  s 


-AU  sin  $  +  AV  cos  $>  =  v 

s  s  s 


(5-19a) 

(5-19b) 


Substituting  Equations  5-19  into  Equations  5-17  yields: 


u. 


=  AU  cos  ($  +  a)  +  AV  sin  ($  +  a)  +  { (b  -  e  )  sin  a 

—  ^  X.  A 


(b  -  e  )  cos  a}  A$ 

y  y 


(5-20a) 


n 


-AU  sin  ($  +  a)  +  AV  cos  ($  +  a)  +  {  (b  -  e  )  sin  a 
s  s  y  y 


+  (b  -  e  )  cos  a}  A$ 
x  x 


(5-20b) 


Using  the  similar  transformation  for  rotations 


u. 


=  AU’  cos  ($  +  a)  +  AV’  sin  O  +  a)  +  {  (b  -  e  )  sin  a 

o  S  A  A 


(b  -  e  )  cos  a}  A$’ 

y  y 


(5-21a) 


. »  _ 


-AU'  sin  ($  +  a)  +  AV 1  cos  ($  +  a)  +  {  (b  -  e  )  sin  a 
s  s  y  y 


+  (b  -  e  )  cos  a}  A$ ' 
x  x 


(5-21b) 


The  axial  displacement  of  C  may  be  expressed  in  terms  of  the  axial  dis¬ 
placement  of  C  as 


w  =Aw-y0~x0  +  w  f  4>" 

—  C  ^o  x  oy  BB 

c 


(5-22) 


where  w  . ,  is  the  unit  relative  warping  of  the  transformed  section 
BB 

sectorial  centroid  (B f )  and  the  sectorial  centroid  of  the  original 

section  (ft);  and,  0  ,  0  are  the  increments  in  rotations  about  refer- 

x  y 

ence  axes  x  and  y . 


' 
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Neglecting  w  ,  we  may  write 

w  =  Aw  -  (y  cos  $  +  x  sin  $)A0  -  (x  cos  $  -  y  sin  $)A0  (5-23) 

—  C  o  oxo  o  y 

c 


or 

w  =  Aw  -  (y  cos  $  +  x  sin  $)A  V*  -  (x  cos  $  -  y  sin  $)AU'  (5-24) 
—  C  o  o  so  o  s 

c 


Equations  5-20,  5-21  and  5-24  may  be  evaluated  at  the  end  nodes  p  and 
q  of  the  element  to  establish  the  transformation  of  displacement  incre¬ 
ments.  However  it  is  also  necessary  to  establish  a  transformation  of 
the  axial  displacement  for  the  mid-point  node.  This  may  be  accomplished 
as  follows. 

Using  the  expression  Equation  5-24  for  the  displacement  of 

the  middle  node  w  ,  and  designating  the  location  with  the  subscript 

c  m 

m,  yields 


w 


c  m 


Aw  -  (y  cos  $  +  x  sin  $  )AV' 

cm  om  m  om  m  sm 


-  (x  cos  $  -  y  sin  <1  )  AU  ’ 

om  m  om  m  sm 


(5-25) 


For  the  assumed  shape  function  Equation  5  25  reduces  to 


w 


cm 


Aw  +  i  (y  cos  $  +  x  sin  $  )  (w-  VP 

cm  t  om  m  om  m  2 


+  ±a  (Pi  -  |  a  vq  +  \  a  e£  l)  +  \  <xom  cos  * 


-  y  sin  $  )  Cv  A  UP  +  , ,  A  0P  t 

J  m2  4  y 


ora 


+  Y  A  0q  l) 
4  y 


(5-26) 


Ill 


The  increments  of  the  displacements  of  the  instantaneous  coordinates 

*  *  # 

may  now  be  written  in  terms  of  those  of  C  and  S  by  the  relation 


^ArE^L  ^ArE^G 


(5-27) 


where  [t2]  is  given  in  Table  5-1;  {Ar£}L  is 


{  rEJL 


iP,  vP,  wL  4>P 


l  y  V 

K  n 


:E}G  1 


,  l 

ur  » 

l 

vP  , 

n,z 

,  l 

ur  » 

l 

vq  , 

n,z 

Equation 

5- 

13b. 

cm 


(5-28a) 


5.4.4  Assembly 

The  element  equilibrium  relationship  may  now  be  expressed  in 
terms  of  the  original  reference  axis  nodal  displacements  by  using  the 
appropriate  displacement  transformation  from  Section  5.4.2  or  5.4.3, 
and  premultiplying  by  the  consistent  force  transformation.  Equation 
5-9  then  becomes 


{Ar^}^  +  ^rE^G 


"S  G  E  G 


{ARE>. 


(5-29) 


where 


[k*]Q  =  [T]T  [k£]  [T] 


(5-30a) 


[kQ]G  -  [T]  [kG]  [T] 


(5-30b) 


and 


{areig  -  [t]  (are}l 


(5-30c) 
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Assembly  of  the  incremental  equilibrium  equations  now  pro¬ 
ceeds  by  the  direct  stiffness  method  to  yield 

[£]  (Ar>  +  [K  ]  {Ar}  =  {AR}  (5-31) 

S 

Equation  5-31  is  now  in  the  form  of  Equation  5-1  and  forms  the  basis 
for  the  solution  technique  discussed  in  Section  5.2  and  those  employed 
in  Chapter  IV. 


5 • 5  Tangent  Stiffness  Properties 

T 

The  evaluation  of  the  element  tangent  stiffness  matrix  [k  ] , 

which  arises  in  Equation  5-9,  as  described  in  Section  5.3,  is  carried 

T  T  T 

out  in  Appendix  E.  Appendix  E  presupposes  that  the  values  of  A  ,  1^ ,  I^» 

m 

X1  e  e  x  y  and  the  angle  a  are  known  at  each  end  of  the  element, 
co  ’  £  n  o  '  o 

For  any  given  set  of  total  displacements  {r},  the  strains  may  be  deter¬ 
mined  at  any  point  on  the  cross  section  from  the  kinematic  relation¬ 
ship,  Equation  A-36 .  The  transformed  section  properties  may  then  be 
computed  as  detailed  in  Appendix  D.  This  provides  all  the  information 
required  for  the  evaluation  of  [k  ] . 


5 • 6  Stress  Resultants 

5.6.1  Stress  Resultants  for  [k  ] 

The  evaluation  of  [K_] ,  which  arises  in  Equation  5-9  as  des- 
cribed  in  Section  5.3,  requires  that  the  total  stress  resultants  be 
known.  The  term  which  arise  in  Section  5.3  are  identical  to  the  terms 
depending  on  the  stress  resultants  in  Equation  2—24,  for  the  elastic 
case,  except  that,  for  'Formulation  1*  the  stress  resultants  and  dis¬ 
placement  increments  are  referred  to  the  C,  S,  x-y  reference  system, 


. 
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while,  for  'Formulation  2'  they  are  referred  to  the  C,  S,  £-n  refer¬ 
ence  system.  All  coefficient  matrices  remain  the  same  as  those  of 
section  A- 6 ,  except  that  the  expression  for  [g  ]  arising  from  the 
evaluation  of  M  is  now  expressed  directly  in  terms  of  this  stress 
resultant . 

For  any  set  of  nodal  displacements  the  strain  is  determined 
from  Equation  A-36.  The  stresses  are  determined  in  each  region  of 
the  plate  segment  from  the  stress-strain  relationship.  The  axial 
load  and  bending  moment  stress  resultants  are  calculated  with  respect 
to  the  reference  axis  passing  through  the  centroid,  C,  of  the  beam 
cross  section  by  simple  statics  as  described  in  Appendix  D.3.  The 

stress  resultant  W  is  also  computed  by  simple  statics  about  S. 

0) 

For  Formulation  1,  these  stress  resultants  are  now  used 
directly  in  the  evaluation  of  [k  ]  as  detailed  in  Appendix  E.  The 
shear  stress  resultants  required  for  the  evaluation  are  obtained  from 


the  equations 


V 


dM 

__X 

dz 


(5-32a) 


dM 


V 


dz 


(5-32b) 


For  Formulation  2,  the  required  stress  resultants  are 
obtained  from  the  above  evaluations  as  (Fig*  5.6) 


M 

=  (M 

-  P  y  ) 

cos  a  -  (M 

-  P  x  ) 

sin  a 

(5-33a) 

5 

X 

o 

y 

o 

M 

=  (M 

-  P  X  ) 

cos  a  +  (M 

-  p  y  J 

sin  a 

(5-33b) 

n 

y 

o 

X 

o 

dM 

V, 

n 

(5-33c) 

5  dz 


. 
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dM 

V  =  -f- 

n  dz 


(5-33d) 


5.6.2  Stress  Resultants  for  Evaluating  Unbalanced  Forces 

The  resisting  stress  resultants  M^,  computed  in  Section 
5.6.1  are  with  respect  to  the  local  reference  axes  moving  with  the 
cross  section.  In  order  to  compute  the  unbalanced  residual  force 
vector  in  the  global  system,  the  stress  resultants,  <li^,  M  ,  M^> , 
with  respect  to  the  local  system  may  be  transformed  to  the  global 
system  using  the  following  transformation  matrix 


X 

y 

z 

X 

cos  $ 

sin  $ 

.  dU 
dz 

Y 

-sin  $ 

cos  $ 

d V 
dz 

Z 

dU 

dz 

.  iY. 

dz 

1 

(5-34) 


Since  the  torsional  moment  is  relatively  small,  and  this 

work,  is  primarily  concerned  with  prebuckling  deformations,  the  con 

tributions  of  the  torsional  moment  to  the  moments  about  the  X  and  Y 

axes  (i.e.  -  M_  and  —  Mj  have  been  neglected  in  the  computations, 

dz  T  dz  i 

The  contribution  of  the  element  to  the  unbalanced  load  vector 


of  Equation  5-31  may  now  be  evaluated  by  adding  the  approximate  nodal 
stress  resultants  as  determined  above,  into  the  applied  load  vector 


(AR>. 


5 . 7  Summary 


In  this  chapter  the  finite  incremental  equilibrium  equations. 
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for  inelastic  beam-columns,  have  been  formulated  with  respect  to  two 
reference  systems.  It  should  be  noted  that  the  incremental  stiffness 
matrix  developed  in  Formulation  2  is  not  exact,  since  the  rotation  of 
the  principal  axes  throughout  the  element  length  is  not  taken  into 
account.  Using  the  formulations  developed  in  this  chapter,  inelastic 
beam— column  and  bifurcation  problems  are  solved  in  Chapter  VI. 


TABLE  5-1  TRANSFORMATION  MATRIX 
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Table  5-1  (cont’d) 
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FIGURE  5.1  ISOLATED  BEAM-COLUMN 
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FIGURE  5.2  CROSS  SECTION  OF  THE  COLUMN  IN  DISPLACED  POSITION 
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FIGURE  5.3  COORDINATE  TRANSFORMATION 
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V 


FIGURE  5.4  TRANSFORMED  SECTION  AXES,  DISPLACEMENTS  AND 
STRESS  RESULTANT  INCREMENTS 
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FIGURE  5.5 


ELEMENT,  NODAL  PARAMETERS 


a  n 
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FIGURE  5.6  TRANSFORMATION  FOR  STRESS  RESULTANTS 


CHAPTER  VI 


SOLUTION  FOR  INELASTIC  MEMBERS 


6 . 1  Introduction 

To  date  no  general  study  of  the  inelastic  behaviour  of  beam- 
columns  bent  out  of  plane  is  available,  but  a  number  of  solutions  for 
special  cases  have  been  obtained.  The  ultimate  load  carrying  capacity 
of  H-columns  under  biaxial  loading  has  been  studied  by  Birnstiel 
(11,12,32).  Their  procedure  requires  successive  trials  and  corrections 
and  is  based  on  an  indirect  incremental  load  analysis  in  which  strain 
history  is  considered.  Recently  Smith  (64)  considered  the  biaxially 
loaded  column  in  terms  of  the  equations  of  three  dimensional  elasticity 
taking  into  account  the  nonlinear  effects  of  the  end  tractions. 
Inelastic  lateral-torsional  buckling  problems  have  been  investigated 
by  Galambos  (24,25),  Fukumoto  (22,23)  and  others  (47).  However,  a 
comprehensive  theoretical  solution  is  not  yet  available  for  general 
loading  conditions. 

In  this  chapter  inelastic  beam-column  problems  and  inelastic 
buckling  problems  are  solved  using  the  formulations  developed  in 
Chapter  V.  Inelastic  beam-column  problems  are  solved  by  Formulation 
2  (discussed  in  Section  5.2.6)  using  an  iterative  incremental  tech¬ 
nique,  and  inelastic  buckling  problems  are  solved  as  an  eigenvalue 
solution  for  critical  length  by  both  Formulations  1  and  2. 

6 . 2  In-plane  Behaviour  of  Inelastic  Beams 

The  term  1 in-plane ?  behaviour  is  used  to  describe  behaviour 
in  which  the  deflection  of  the  shear  centre  is  in  a  principal  direction 
of  the  cross  section.  In  order  to  verify  the  model,  and  the  solution 
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technique,  the  in--plane  behaviour  of  simple  beam-columns  has  been 
studied.  The  bending  moment  acting  at  any  section  of  a  beam-column 
is  composed  of  the  primary  moment  Mq  (i.e.  the  static  moment  computed 
from  original  geometry),  and  the  secondary  moment  (i.e.  moments  result¬ 
ing  from  a  change  in  the  equilibrium  equations  due  to  the  deflection 
of  the  member) .  For  in— plane  problems  the  secondary  moment  is  expressed 
as  Pv  where  v  is  the  deflection  in  the  principal  direction.  Assuming 
a  portion  of  the  member  becomes  plastic,  the  deflection  v  increases 
and  the  secondary  moment  becomes  a  proportionately  larger  component 
of  the  resisting  moment.  Eventually  the  beam-column  is  unable  to 
support  an  increase  of  primary  moment  and  thereafter  the  value  of  pri¬ 
mary  moment  decreases  as  the  deflection  increases.  For  this  type  of 
problem,  lateral  and  local  buckling  are  assumed  to  be  prevented,  and 
the  strains  and  displacements  are  assumed  to  be  small. 

The  incremental  finite  element  equation  of  Chapter  V  expressed 
as  Equation  5-1  and  Equation  5-31,  may  now  be  applied  to  beam-column 
problems  by  using  the  procedure  outlined  in  Section  5.2.1. 


For  any  iterate  Equation  5—1  can  be  written  as 


(6-1) 


n 


where  { E }  is  the  vector  of  unbalanced  residual  nodal  forces  for  the 


n 


current  configuration.  This  vector  is  evaluated,  as  described  in 
Section  5.2.1,  as 


{E}n  =  {R}  -  {?}n  -  [KG]n  (r)n 


(6-2) 


where  (R)  is  the  vector  of  total  applied  loads;  (F>n  is  the  vector  of 


. 
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resisting  forces  resulting  from  the  stress  resultants  in  the  current 
configuration,  computed  as  described  in  Section  5.6.2;  and  [K^] 
represents  the  vector  of  forces  resulting  from  the  change  in  geometry 
from  the  original  configuration,  as  discussed  in  Section  5.2.1. 

In  the  case  of  in-plane  beam-column  response,  with  constant 
axial  load,  [K_]  is  constant.  Wien  {R}  is  applied  to  the  structure, 
the  response  is  {r}^  =  {Ar}^  for  tlie  first  iterate,  which  can  be 

computed  from  equation 

[[K g]x  +  [KG]]{A  rL}  =  (E>1  =  {R)1  (6-3) 

After  the  first  and  succeeding  iterates ,  nodal  point  displacements  are 
updated  by  adding  the  incremental  displacements  to  the  previous  dis¬ 
placements.  Strains  and  stresses  are  then  computed  from  the  total 
displacements . 

In  general  the  incremental  equation  is  given  by  Equation  6-1 
and  iteration  is  repeated  until  the  structure  is  in  equilibrium  with 
the  applied  forces.  In  the  solution  scheme  [K  ]  is  kept  unchanged  for 
some  iterations,  as  discussed  in  Section  5.2.1,  until  the  rate  of  con¬ 
vergence  begins  to  deteriorate,  at  which  time  the  matrix  is  updated 
with  properties  based  upon  the  current  displacements.  This  is  a 
modified  Newton-Raphson  procedure.  Since  the  nature  of  the  problem  is 
very  sensitive  to  the  applied  forces  and  the  type  of  stress— strain 
curve,  an  under-relaxation  factor  was  used  on  the  residual  force  vector. 
The  procedure  may  diverge  if  the  matrix  is  not  updated  often  enough. 

This  was  also  observed  by  Haisler  (30)  and  Stricklin. 

As  has  already  been  stated  in  Section  5.3  a  linear  variation 
of  axial  displacements  is  incapable  of  establishing  equilibrium  of 
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axial  forces  with  variable  section  properties. 

To  achieve  a  balance  of  axial  load  at  the  end  nodes  of  the 

elements,  an  additional  degree  of  freedom  at  the  mid-point  of  an 

element  was  introduced.  Equilibrium  is  obtained  as  follows.  Consider 

an  element  k  whose  resisting  axial  forces  at  node  p  and  q  are  denoted 

bv  Pk  and  P^  and  an  adjacent  element  (k+1) ,  whose  axial  forces  are 
P  q 

P^+and  pk+1.  The  unbalanced  force  at  the  k+lttL  node  is  -  P  .  For 

p  q  p  q 

inelastic  problems,  P  ^  P  for  a  linear  variation  of  curvature  and  a 

^  k  k 

constant  axial  strain.  Hence  the  unbalanced  force  P^  -  P^  on  the 

element  k  was  equilibrated  by  a  force  at  the  midpoint  of  the  element  k. 

Iteration  was  carried  out  until  this  unbalanced  axial  force  at  the 

center  node  was  reduced  to  an  arbitrarily  small  quantity  and  the  axial 

loads  at  each  end  node  were  balanced.  It  was  possible  to  achieve  this 

condition  because  of  the  linear  variation  of  axial  strain.  The  axial 

stiffness  with  this  additional  degree  of  freedom  is  derived  in  Appendix 

A. 

6 . 3  Solutions  of  Typical  In-plane  Inelastic  Problems 

E-l  Inelastic  Bending  of  Simple  Beam-Column  of  Rectangular  Cross  Section 
An  inelastic  bending  analysis  was  carried  out  for  a  cantilever 
beam  of  rectangular  cross  section,  with  dimensions  and  loading  as  shown 
in  Fig.  6.1.  The  stress-strain  curve  of  Fig.  D.la  was  used  in  the 
numerical  solution  of  all  problems  with  the  values  of  E,  Eg  shown  in 
Table  6— 3a.  The  element  idealization  is  also  shown  in  the  figure.  For 
a  constant  axial  load  and  a  specified  displacement  at  the  free  end, 
the  reaction  Q  at  the  free  end,  where  the  transverse  displacement  is 
specified,  determines  one  half  the  loading  corresponding  to  the 
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displacement  of  a  simply  supported  beam  column  with  a  concentrated 

load  at  the  center.  The  load  vs  end  deflection  curves  are  shown  in 

Fig.  6.1  for  the  cases  of  P/P^  =  0;  P/Py  =  0.152;  P/Py  =  0.632. 

Moment  curvature  relationships  are  also  shown  and  compared  with 

theoretical  values  in  Fig.  6.2.  The  results  are  compared  with  second 

order  elastic  and  second  order  rigid  plastic  solutions  in  Fig.  6.1. 

E-2  Inelastic  Bending  of  a  Simple  WF  Beam-Column 

As  in  example  E-l,  one  half  of  a  simply  supported  beam- 

column  with  a  concentrated  load  at  the  center  was  simulated  by  the 

cantilever  shown  in  Fig.  6.3.  Displacement  was  again  specified  at 

the  free  end.  The  figure  shows  the  load  vs  deflection  plot  for  both 

strain  hardening  and  elas tic-ideally  plastic  cases,  and  compares  the 

results  with  the  second  order  rigid  plastic  solution. 

E-3  Beam-Column;  Moment  at  one  end 

The  behaviour  of  a  beam— column ,  shown  in  Fig.  6.4,  is 

represented  by  relating  the  applied  moment,  Mq ,  to  the  end  rotation, 

6  .  The  beam  was  subjected  to  an  axial  load  of  0.49  P  ,  which  was 
o  y 

kept  constant.  An  end  moment  Mq  was  applied  about  the  major  axis  of 

the  member  at  one  end.  The  solution  was  obtained  by  specifying  the 

displacement  0q  and  iterating  with  this  as  a  displacement  boundary 

condition.  Instability  is  demonstrated  by  the  decrease  in  moment 

capacity  beyond  the  maximum  value.  The  figure  shows  the  comparison 

of  experimental  and  computed  end  rotation  curves.  The  discrepancy 

between  theoretical  and  experimental  results  is  consistent  with  that 

of  Kauren  (73)  and  Galambos  who  reported  a  computed  maximum  value 

for  M  /M  of  0.55  compared  with  the  experimental  value  of  0.60. 
o  p 
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E - 4  Beam-Column  with  equal  End  Moments 

The  section  analyzed  was  an  8WF31  with  residual  stress 

distribution  as  shown  in  Fig.  6.5.  P  was  kept  constant  at  0.7Py. 

Hence  the  column  is  inelastic  from  the  beginning  of  loading.  Half 

the  beam  is  considered  because  of  symmetry.  End  A  was  free  to  move 

vertically  and  loading  was  applied  by  imposing  an  end  rotation. 

Very  small  increments  in  rotation  were  specified  at  this  end  and 

iteration  was  carried  out  to  find  the  equilibrium  configuration. 

The  end  moment  -  rotation  curve  is  shown  in  Fig.  6.5  and  compared 

with  the  curve  given  in  the  Lehigh  University  design  aids  (57) . 

(The  Lehigh  University  design  aid  is  based  on  a =  36ksi, 

a  =  0.3  a  ,  and  no  strain  hardening.)  Convergence  of  this  problem 
rc  y 

was  slow  and  hence  the  solution  scheme  was  not  carried  out  beyond  the 
peak. 

6.4  Solution  of  General  Beam-Column  Problems 

In  the  examples  of  Section  6.3,  all  displacements  except 
those  in  the  y-z  plane  were  suppressed.  Some  examples  of  the  more 
general  beam— column  problem,  in  which  out— oi— plane  displacements  also 
occur,  are  presented  in  this  section. 

E-l  Twist  of  An  Inelastic  Beam-Column 

In  order  to  test  the  behaviour  of  the  model  when  inelastic 
torsional  response  occurs,  a  cantilever  WF  beam,  as  shown  in  Fig.  6.6, 
was  subjected  to  applied  axial  loads  of  40k  at  each  node  and  a  torque, 
Mt,  applied  at  the  free  end.  The  analysis  of  this  problem  was  prin¬ 
cipally  carried  out  to  test  the  convergence  of  the  bimoment  in  the 
pj^gi^gtic  range.  No  theoretical  solution  is  available  for  compar ison , 
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but  convergence  of  unbalanced  forces  was  achieved. 

E - 2  Inelastic  H-Columns  under  Biaxial  Bending 

A  14WF43  section  with  a  doubly  eccentric  axial  load,  applied 

at  accentricities  of  e  =0.5"  and  e  =  5",  and  L/t  =  H7,  was  analysed 

x  y  y 

and  compared  with  the  results  obtained  by  Harstead  (32)  and  Birnstiel, 
et.  al. ,  (11).  Their  procedure  is  given  in  Table  6-1.  This  problem 
involves  coupling  of  axial  deformations,  bending  deformations  about 
two  axes,  and  nonuniform  torsional  deformation  in  the  inelastic  range 
of  Beam-column  response.  The  load  was  incremented  until  the  displace¬ 
ments  no  longer  converged.  Fig.  6.7  shows  the  load  and  midpoint  dis¬ 
placement  relationships.  Unlike  the  case  of  the  problems  discussed  in 
Section  6.3,  boundary  displacement  increments  cannot  be  adopted  for 
loading,  since  the  adjustment  of  the  end  displacements  in  the  proper 
ratios  is  difficult.  The  procedure  of  incrementing  the  forces  yielded 
results  until  a  maximum  load  was  attained,  but  numerical  difficulties 
did  not  permit  an  evaluation  in  the  unloading  region.  The  results 
agree  well  with  Birnstiel* s  analysis. 

E-3  Inelastic  Lateral  Response  with  Initial  Imperfections 

An  8WF31  beam-column  with  the  residual  stress  distribution 
shown  in  Fig.  6.8  was  analyzed  while  subjected  to  a  constant  axial  load 
of  120  kips  (P/P  =0.4)  and  with  an  initial  lateral  displacement  of 
u  =  0.01  sin  it  z/L.  A  beam-column  solution  was  carried  out  until,  with 
an  end  moment  of  440  inch-kips,  the  residual  force  resultants  no  longer 
converged.  Moment  vs  in-plane  and  out-of-plane  response  is  shown  in 
Fig.  6.8.  For  the  same  problem  the  eigenvalue  solution  yields  a  criti 
cal  length  of  100"  for  end  moments  of  480  inch-kips,  and  this  is  shown 
as  an  upper  bound  in  the  figure. 


. 
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6 . 5  Inelastic  Buckling 
6.5.1  Formulation  of  Problem 

Inelastic  buckling  (or  bifurcation)  analyses  can  be  carried 
out  if  it  is  assumed  that,  as  for  the  beam-column  problems,  a  reason¬ 
able  estimate  of  the  load  carrying  capacity  may  be  obtained  by  employ¬ 
ing  Shanley’s  (62)  tangent  modulus  concepts.  The  basic  assumption  is 
that  no  strain  reversal  occurs,  and  therefore  the  effective  incremental 
stiffness  at  any  section  may  be  determined  from  the  instantaneous  tan¬ 
gent  modulus  values  at  all  parts  of  the  cross  section.  Since  the  stress 
in  the  member  is  essentially  uniaxial, the  effective  strain,  for  the 
determination  of  tangent  moduli,  is  considered  to  be  the  longitudinal 
strain.  The  incremental  section  properties  may  then  be  computed  by  a 
transformed  area  concept,  as  discussed  in  Section  5.2.3  and  5.5. 

In  addition  to  the  above,  the  following  assumptions  are  made 

for  bifurcation  problems: 

1.  The  members  are  initially  straight  and  prismatic. 

2.  The  projection  of  the  cross  section  on  a  plane  normal  to  the 
centroidal  axis  does  not  distort. 

3.  Displacements  of  a  point  may  be  obtained  by  superimposing  warping 
displacements  on  plane  section  displacements,  each  of  which  are 
consistent  with  simple  beam  theory  (54,68,74). 

4.  The  residual  stress  distribution  satisfies  statics. 

5.  Variation  of  the  longitudinal  stress  across  the  plate  thickness 
may  be  neglected. 

6.  The  effect  of  prebuckling  displacements  on  the  equilibrium  equa¬ 
tions  may  be  neglected. 

7.  Boundary  conditions  are  such  that  the  structure  remains  statically 
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determinate. 

8.  The  stress-strain  curve  is  trilinear  (Fig.  D.la). 

It  should  be  noted  that  it  is  due  to  assumptions  1)  and  6) 
that  the  more  realistic  but  numerically  more  difficult,  beam-column 
problem  is  converted  to  a  bifurcation  problem.  The  remaining  assump¬ 
tions  are  the  same  as  those  used  for  beam-columns ,  except  for  (7). 

A  bifurcation  loading  condition  has  been  attained  if  it  is 
possible  to  determine  nontrivial  solutions  to  Equation  5-31,  for 
{AR}  =  0.  The  condition  required  to  satisfy  this  requirement  is  that 
the  determinant  of  the  coefficient  matrix  is  equal  to  zero,  i.e. 


det  {[Kg]  +  [Kg]}  =  0 


(6-4) 


If  [Kg]  were  independent  of  load  and  [Kg]  were  linearly 
dependent  on  load  the  critical  loading  could  be  obtained  from  a  stand¬ 
ard  eigenvalue  analysis.  However,  for  an  arbitrary  cross  section, 

T 

subjected  to  arbitrary  loads  the  matrix  [Kg]  is  highly  sensitive  to  the 
loading  condition  after  inelastic  response  has  been  initiated. 

The  critical  loading  may  now  be  determined  in  a  number  of 
ways.  Fukumoto  (22,23)  has  evaluated  the  determinant  and  extrapolated 
or  interpolated  this  value  to  determine  the  loading  for  which  Equation 
6-4  is  satisfied.  Harris  and  Pifko  (31)  have  used  an  iterative 
approach  on  the  load  and  found  the  load  level  for  which  the  eigenvalue 
was  one.  The  procedure  used  in  this  dissertation  is  to  iterate  for  an 
eigenvalue  to  determine  the  critical  length  rather  than  the  critical 
load . 


6. 5 . 2  Inelastic  Buckline  by  Direct  Method 


, 
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For  any  arbitrary  statically  determinate  loading,  equilibrium 
is  established,  including  the  effects  of  inelastic  material  response 
but  neglecting  the  effect  of  displacements  on  the  equilibrium  equations, 
by  an  iterative  procedure  detailed  in  Table  6.2.  Knowing  the  strain 
throughout  the  member  the  tangent  stiffness  matrix  [K  ]  can  be  esta¬ 
blished  as  outlined  in  Section  5.5  and  detailed  in  Appendix  E.  The 
prebuckling  stress  resultants  can  also  be  established  by  simple  inte¬ 
gration  as  detailed  in  Appendix  D3.  This  permits  the  evaluation  of 

[K  ]  as  detailed  in  Appendix  E. 

G 

An  examination  of  the  stiffness  matrices  in  Appendix  E  shows 
that  Equation  6-4  may  be  written  as 

[[K1]  +  [K2]  +  [K3]]  {Ar}  +  [Kg]  {Ar}  =  {AR}  (6-5) 

2 

where  and  contain  terms  dependent  on  the  inverse  of  L ,  L 

and  l ^  respectively.  Noting  that  [K  ] is  proportional  to  the  inverse 
of  L ,  defining  the  critical  length  as  =  XL  ,  and  factoring  X; 
the  condition  for  the  existence  of  nontrivial  solutions  of  Equation 
6-5  with  {AR}  =  0  may  be  written  as 

~  [[K  ]  +  x[K  ]  +  X2[k  ]]  {Ar}  =  -  [Kg]  Ar  (6-6) 

X 

This  equation  may  be  abbreviated  as 

— y  [[K  (  X.  )]]  Ar  =  -  [Kq3  Ar  (6-7) 

Xi+1 

For  any  general  problem  under  any  general  loading,  assuming  a  value  for 

X. ,  an  eigenvalue  solution  will  yield  X •  Equation  6-7  may  be  iter 

ated  until  X.  =  X.  at  which  time  the  critical  length  has  been 
l+l  l 
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determined.  When  the  equations  corresponding  to  <5u,  5V  and  6^ 

m 

[K  ]  matrix  are  uncoupled  from  that  of  6w,  these  uncoupled  equations 

s 

do  not  depend  on  X.  Hence  the  first  iterate  of  Equation  6-7  gives  the 
exact  solution  for  critical  length.  When  the  tangent  properties  are 
assumed  to  be  the  same  all  along  the  length,  then  also  the  first 
iterate  will  lead  to  an  exact  solution  for  critical  length.  Only  such 
problems  are  considered  in  this  dissertation  since  no  solutions  with 
the  more  general  coupling  are  available  for  comparison. 

It  should  be  noted  that  the  procedure  described  for  the 
solution  of  Equation  6-7  is  only  valid  if  [K  ]  is  proportional  to  the 
inverse  of  L  and  therefore  the  stress  resultants  must  remain  constant 
at  each  node  as  the  member  length  is  scaled.  A  simultaneous  scaling 
of  the  transverse  loading  is  required  to  satisfy  this  condition. 

Since  stress  resultants  are  not  proportional  to  loading  for  statically 
indeterminate  inelastic  structures  the  method  is  applicable  only  to 
statically  determinate  structures. 

6.5.3  Inelastic  Buckling  by  a  Beam-column  Approach 

In  order  to  account  for  prebuckling  deformation,  a  beam- 
column  equilibrium  configuration  can  be  established  for  any  trial 
loading,  by  the  procedure  discussed  in  Section  6.2.  Knowing  the 
deflections,  the  total  moment  at  any  section  may  be  computed  and  the 
tangent  flexural  stiffness  and  geometric  stiffnesses  may  be  evaluated 
as  discussed  in  Section  6.5.2.  Solving  the  eigenvalue  problem,  a 
critical  length  can  be  determined  from  one  iterate  of  Equation  6-7. 

If  this  critical  length  is  the  same  as  the  actual  length  of  the  column 
the  trial  loading  is  the  critical  one.  If  the  critical  length  is 
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greater  (less)  than  the  actual  length,  the  true  loading  level  is  less 
(greater)  than  the  critical  state.  In  this  case  the  load  may  be 
incremented  (decremented)  and  the  resulting  eigenvalue  problem  is 
solved  by  successive  trials  until  X  =  1.  This  procedure  is  time  con¬ 
suming  and  expensive  and  hence  has  not  been  generally  adopted  in  the 
present  investigation.  However,  the  procedure  was  used  for  example  . 
E-2  of  Section  6.4  and  the  results  are  shown  in  Fig.  6.19  where  they 
are  compared  with  those  from  the  methods  of  Section  6.5.2. 


6 . 6  Solution  of  Inelastic  Buckling  Problems 

E-l  Flexural  and  Lateral  Torsional  Buckling  of  Double  Angle  Strut 

Double  angle  struts  are  used  in  building  structures  as  com¬ 
pression  members  in  trusses.  They  are  usually  subjected  to  axial  load 
and  the  secondary  bending  moment  arising  from  the • connections  is  neg¬ 
lected.  The  details  of  the  section  analysed  and  the  residual  strain 
distribution  are  shown  in  Fig.  6.9a  and  Table  6-3b .  The  St-Venant 
torsional  stiffness,  GK^,  was  kept  at  the  elastic  value  in  order  to 
compare  with  the  results  obtained  by  Nuttall(53) .  Nuttal's  solution 
was  obtained  by  imposing  a  uniform  axial  strain,  evaluating  transformed 
section  properties  numerically  and  solving  the  normal  stability  equa¬ 
tions  for  critical  length.  For  this  particular  problem  the  yy  axis 

buckling  strength  virtually  coincided  with  the  lateral  torsional  buck- 

2  2 

ling  strength,  since  (1-e^/r^)  is  approximately  one  throughout  the 

loading  range.  Fig.  6.9b  shows  the  plot  of  a/o^  values  against  the 

slenderness  ratio,  L/r  ,  for  lateral  torsional  buckling  and  xx  axis 

x 

buckling.  The  influence  of  the  residual  stress  and  the  significance 
of  the  inelastic  response  is  apparent.  The  results  obtained  from  the 
finite  element  analysis  agree  with  the  results  of  Nuttall. 
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E-2  Flexural  and  Lateral  Torsional  Buckling  of  Hat-shaped  Sections 

Hat  shaped  sections  are  often  used  as  chord  members  of  open 
web  steel  joists.  The  dimensions  of  the  section  analyzed,  the  normal¬ 
ized  warping  coordinates  and  the  residual  stress  distribution  are  shown 
in  Fig.  6.10a,  b  and  c,  respectively.  Critical  buckling  occurs  about 
the  minor  axis  as  shown  in  Fig.  6.11.  Critical  lengths  for  lateral 
torsional  buckling  are  also  computed  and  compared  with  the  results 
obtained  by  Heaton  (33). 

E-3  Inelastic  Lateral  Buckling  of  a  Doubly  Symmetric  Section 

The  inelastic  beam  buckling  curve  for  an  8WF31  section  sub¬ 
jected  to  equal  end  moments  is  shown  in  Fig.  6.12.  This  curve  consists 
of  three  parts.  The  first  and  last  part  of  the  curve  represent  classi¬ 
cal  elastic  buckling  with  E=E  and  E=E  respectively.  The  intermediate 
portion  of  the  curve  represents  buckling  in  the  inelastic  range  and  is 
applicable  when  some  parts  of  the  cross  section  have  yielded  while 
other  parts  are  elastic.  The  strain  hardening  and  elastic  curves  are 
hyperbolas  which  do  not  intersect.  The  curve  for  inelastic  buckling 
provides  a  transition  between  these  two  idealizations.  The  St.  Venant 
constant  was  assumed  to  vary  linearly  with  the  transformed  area.  A 
solution  of  the  problem  was  obtained  by  Galambos  (25).  The  importance 
of  the  residual  stress  distribution  is  again  demonstrated  and  the 
comparison  is  reasonably  good. 

E-4  Inelastic  Lateral  Buckling  of  Mono-Symmetrical  Sections 

The  double  angle  section  of  E-l  was  analysed  for  inelastic 
lateral  buckling.  The  inelastic  buckling  curve  is  shown  in  Fig.  6.13. 
Only  a  theoretical  elastic  buckling  solution  is  available  for  comparison. 


' 
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E-5  Inelastic  Lateral  Torsional  Buckling  of  an  Unsymmetrical  Section 
(example  1) 

In  the  case  of  mono  or  doubly  symmetric  sections  the  response 
of  the  beam-column  is  in-plane  until  the  buckling  occurs.  But  in  the 
case  of  unsymmetrical  sections  the  beam  will  deflect  laterally  and 
twist  from  the  start  of  the  application  of  the  load.  The  eigenvalue 
solution  of  the  problem  will  give  an  upper  bound  for  the  load  since 
prebuckling  deformations  are  neglected.  Fig.  6.14  shows  the  inelastic 
lateral  torsional  buckling  of  an  angle  section  (67)  subjected  to  end 
moments.  The  experimental  point  of  Thomas  and  Leigh  is  shown  on  the 
plot.  The  discrepancy  is  probably  due  to  the  fact  that  prebuckling 
deformations  have  been  neglected. 

E-6  Inelastic  Lateral-Torsional  Buckling  of  an  Unsymmetrical  Section 
(example  2) 

Results  for  a  similar  type  of  problem  to  example  E-5  are 
shown  in  Fig.  6.15,  for  the  inelastic  lateral  torsional  buckling  of 
an  unsymmetrical  section  when  one  channel  is  placed  over  the  other. 

An  elastic  solution  to  this  problem  was  obtained  (in  E-4  of  Section 
3.3.3).  No  comparitive  solution  is  available. 

E-7  Inelastic  Lateral  Torsional  Buckling  of  Beam-Columns  (Axial  Load- 
Equal  Eccentricity) 

Fig.  6.16  shows  a  comparison  of  the  finite  element  analysis 
for  inelastic  torsional  buckling  of  an  8WF31  beam-column,  under  an 
axial  load  with  constant  eccentricity,  with  that  of  Galambos  (24). 
Since  the  secondary  moment  due  to  prebuckling  deformation  is  neglected 
the  inelastic  properties  are  assumed  to  be  the  same  at  all  sections. 
Hence  properties  for  only  one  section  need  be  calculated.  The  results 
of  the  finite  element  analysis  agree  well  with  those  of  Galambos. 
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E - 8  Inelastic  Lateral  Torsional  Buckling  of  Beam-Column  with  Moment 
Gradient 

Fig.  6.17  compares  the  critical  end  moment  required  to  pro¬ 
duce  lateral  torsional  buckling  of  an  8WF31,  subjected  to  constant 
axial  load,  to  that  obtained  by  Fukumoto  (22,23).  In  contrast  to  the 
preceding  comparative  solutions,  Fukumoto1 s  solution  included  the 
effect  of  in-plane  prebuckling  displacements  in  arriving  at  an  equili¬ 
brium  configuration.  The  critical  load  is  then  determined  as  that  for 
which  the  determinant  vanishes.  In  the  finite  element  analysis  an 
eight  elements  idealization  is  used  and  the  results  are  compared  with 
Fukumoto’ s  analysis.  The  effect  of  prebuckling  displacements  is 
apparently  small. 

E-9  Inelastic  Lateral  Buckling  of  a  Wide  Flange  Beam  Subjected  to  a 
Central  Load 

In  Fig.  6.18  the  critical  moment  vs  slenderness  ratio  curve 
is  shown  for  an  18WF50  beam  subjected  to  a  central  load.  A  large 
reduction  in  the  buckling  strength  results  when  the  presence  of 
residual  stress  initiates  inelastic  behaviour.  The  curves  are  com¬ 
pared  with  Massey's,  et.  al. ,  (48)  recommended  curve  obtained  from  a 
theoretical  analysis.  Massey  used  Mpp/M  =  1.05  to  1.06  and  Mp/M^  - 
1.14  to  1.15,  which  are  applicable  to  most  wide  flange  sections,  where 

M  is  the  bending  moment  at  which  the  flanges  are  fully  yielded.  In 
FP 

Massey’s  analysis,  residual  stress  distribution  is  not  taken  into 
account.  In  the  present  analysis,  the  computed  points  assume  the 
loads  act  through  the  shear  centre.  When  the  load  is  considered  to 
act  at  the  top  or  bottom  flange  the  problem  must  be  solved  by  a  trial 
and  error  solution  for  the  critical  length. 
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E-10  Inelastic  Lateral  Torsional  Buckling  with  Biaxial  Bending 

An  inelastic  buckling  analysis  was  carried  out  for  the  beam- 
column  problem  of  example  E-2  of  Section  6.4.  Fig*  6.19  shows  a  plot 
of  critical  loads  for  biaxial  bending  of  the  section  determined  from  an 
eigenvalue  analysis  and  neglecting  prebuckling  deformations.  The 
curves  are  drawn  for  two  cases:  a)  without  residual  stress  b)  with 
residual  stress.  In  addition  to  the  eigenvalue  analysis  of  Section 
6.5.2,  the  method  of  Section  6.5.3  was  also  applied.  For  a  column 
length  of  221"  (L/r  =  37.6)  and  a  load  of  125  kips  the  critical  length 
obtained  was  283  inches,  showing  that  the  load  level  was  below  the 
ultimate  load.  At  P  =  134  kips  the  critical  length  was  found  to  be 
226  inches  indicating  that  134  kips  was  very  near  to  the  ultimate  load. 
However,  the  critical  load  based  upon  the  eigenvalue  analysis  and 
neglecting  prebuckling  deformation  (i.e.  the  method  of  Section  6.5.2), 
is  27%  greater  than  beam-column  result  obtained  by  either  the  method 
of  Section  6.5.3  or  6.4.  In  such  situations  the  effect  of  prebuckling 
displacement  cannot  be  ignored. 

6 . 7  Summary  and  Conclusions 

In  this  chapter  the  inelastic  equations  derived  in  Chapter  V 
have  been  applied  to  a  variety  of  problems.  Inelastic  beam  column 
problems  were  solved  by  a  modified  Newton— Raphson  technique.  Although 
this  technique  was  shown  to  yield  good  results  for  in-plane  problems 
and  yields  results  for  a  more  general  biaxial  bending  beam-column 
problem,  which  corresponds  to  the  only  known  available  solution  of 
this  problem,  the  technique  is  difficult  to  apply.  Problems  solved  in 
this  way  require  considerable  computing  time  and  it  is  necessary  to 
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input  information  controlling  the  load  step  size  and  under-relaxation 
parameters,  in  addition  to  a  criterion  for  updating  the  coefficient 
matrices . 

The  solution  of  inelastic  stability  problems  neglecting  pre¬ 
buckling  deformation,  was  shown  to  yield  results  which  compare  well 
with  existing  solutions  in  the  literature,  for  a  wide  variety  of  appli¬ 
cations.  The  determination  of  maximum  load  carrying  capacity,  including 
the  effects  of  prebuckling  deformations  by  the  method  of  Section  6.5.3, 
yields  the  same  results  as  the  beam-column  solution  of  Section  6.5.4, 
and  at  the  same  time  yields  an  upper  bound  to  the  carrying  capacity  at 
any  stage  of  the  analysis. 
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Birnstiel 

1. 

2. 

3. 


4. 

5. 

6. 

7. 

8. 
9. 

10. 

11. 

12. 


TABLE  6-1 

s  Procedure  for  Solution  of  H--Columns  Under  Biaxial  Bending 
Assume  trial  values  of  curvatures  corresponding  to  M^}  , 

W  and  axial  strain 
go 

Numerically  integrate  the  curvatures  to  determine  the  dis¬ 
placed  position  of  the  column 

Calculate  the  values  of  internal  bending  moment  and  axial 
load  at  each  station.  Because  the  values  of  axial  force 
differ  (in  the  inelastic  range)  select  the  value  at  mid 
point  as  the  control  value  of  P 

On  the  basis  of  the  displaced  position  and  P  compute  the 
external  bending  moment  at  each  station 

Compare  external  moments  with  internal  moment.  If  they  do 

not  agree  adjust  curvatures  corresponding  to  M  and  M 

x  y 

Repeat  steps  2  to  5  until  convergence  is  attained  in  step  5 
Determine  external  torque  at  each  station 
Compute  internal  torque  at  each  section 

Compare  external  torque  to  internal  torque.  If  they  do  not 
2  2 

agree  adjust  d  <f>/dz  and  integrate  numerically  to  determine 

<j> 

Repeat  7  to  9  until  convergence  is  attained  in  step  9 
Repeat  6  and  10  until  convergence  is  attained  in  5  and  9 
Repeat  steps  1  to  11  until  the  desired  displacement  or 
ultimate  carrying  capacity  is  attained 


' 


' 
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Step 


Note 


Step 


TABLE  6-2 

ALGORITHM  FOR  DETERMINING  INELASTIC  CRITICAL  LENGTHS 


.  Determine  Static  Stress  Resultants 

(a)  Subdivide  the  member  into  elements  of  equal  length. 

(b)  Compute  the  elastic  section  properties  A,  1^  Iy,  Ixy» 


I  and  I  and  I  , 
co  x  y 


where 


I 

x 


I 

X 


/I  I  ) 
x  y 


(6-8a) 


A 

I 

y 


i 

y 


/i  i  ) 

X  y 


(6-8b) 


With  C  and  S  located  at  the  centroid  and  shear  centre 
of  the  section,  (Fig.  2.1). 

(c)  Solve  for  displacements,  with  [K^]  =  0  and[Kg]  deter¬ 
mined  from  elastic  properties. 

(d)  Compute  stress  resultants  P,  Mx,  M^  and  at  each 
node  for  this  small  deflection  solution. 

:  The  stress  resultants  determined  in  this  way  are  those  which 

equilibrate  the  external  forces. 


2.  Determine  Strains  (Iteratively)  Which  Produce  Static  Stress 
Resultants 

(a)  Compute  strains  from  the  initial  approximate  displace¬ 
ments  ,  as 


e  =  w'  -  y  v"  -  x  u"  +  io  <j>"  +  ^ 

t  c  y  s  s  K 

(b)  Calculate  stress  resultants,  P,  Hx>  My  and  W^,  (see 
Appendix  D)  for  the  given  strain  distribution. 


' 


4 
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Table  6 


-2  (cont’d) 


(c)  Compute  the  unbalance  in  stress  resultants,  as 


P  =  P  -  P 


(6-10a) 


M  =  M  -  M 

X  XX 


M  =  M  -  M 

y  y  y 


w  =  w  -  w 

(0  (0  tO 


(6-10b) 

(6-10c) 

(6-10d) 


(d)  Compute  equivalent  unbalanced  stress  resultants,  P  , 

JL  X  % 

M  ",  M  and  W  about  C  and  S,  where 
x  y  to 


P  =  P 


(6-lla) 


*  — 

W  =  W 
to  to 


M  =  M  -  M  I  /I 
x  x  y  xy  y 


(6-llb) 

(6-llc) 


M  =  M  -  M  I  /I 
y  y  x  xy  x 


(e)  Estimate  strain  increments  from  the  equation 


*  M  M 

As.  -  l  +-A-y+J- 

i  A  $  $ 


W 

x  H - = —  to 


to 


(6-12) 


and  a  new  approximate  strain,  as 


ei+l  =  £1  +  Aei 


(6-13) 


(f)  Iterate  on  steps  (2-b)  to  (2-e)  until  the  unbalanced 
stress  resultants  are  negligible. 


. 
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Table  6-2  (cont’d) 

Note:  A  set  of  strains  have  now  been  determined  which  will  provide 

stress  resultants  to  equilibrate  the  external  forces,  assuming 
the  prebuckling  displacements  have  negligible  effect  on  the 
equilibrium  equations. 

Step  3.  Computation  of  Tangent  and  Geometric  Stiffness  Matrices 

T 

(a)  The  tangent  stiffness  matrix,  [K  ] ,  of  Equation  6-4  may 

s 

now  be  evaluated  as  outlined  in  Section  5.3.1  and 
detailed  in  Appendices  D  and  E. 

(b)  The  geometric  stiffness  matrix,  [K  ] ,  of  Equation  6-4 
may  now  be  evaluated  as  outlined  in  Chapter  5 
detailed  in  Appendices  D  and  E. 

Step  4.  Solution  for  Critical  Length 

(a)  Equation  6-7  is  solved  for  X,  to  determine  the  critical 
length  L  =  XL,  as  described  in  the  solution  procedure. 


TABLE  6-3a  PARTICULARS  OF  BEAM-COLUMN  PROBLEMS 
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TABLE  6-3b  PARTICULARS  OF  BUCKLING  PROBLEMS 
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For  the  nomenclature  for  loading  see  Table  6-3c 
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CASES 

LOADING 

A 

P  ka  7$  P 

B 

ft  ^ 

“U  ir 

C 

D 

M .  \M 

K  VVL  ~  -  -  AM  1  r 

E 

^  Jk%’ 

F 

iFy 

G 

|Fy 

r\  K.  V  .  -A,_  p 

H 

.  <J“j  P  <h  P  <J  j  P  <1  j  P 

/  - — -  /v»y 

P 

I 

M  J,  V-J  !  M  y  My  M  X 

VKy  y  yy\K 

TABLE  6.3c  NOMENCLATURE  FOR  LOADING 


2nd  Order  Elastic 


143 


OQ 

UJ 

zn 

o 

o 

i- — 

o 

UJ 

__j 

Li- 

LU 

o 


H 


Sdl)l  Nl  0  QV01 


FIGURE  6.1  LOAD  VS  END  DEFLECTION  FLO' 
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FIGURE  6.2  MOMENT  VS.  CURVATURE  RELATIONSHIP 


Second  Order  Plastic 
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FIGURE.  6.3  LOAD  VS  END  DEFLECTION  PLOT  FOR  WF  BEAM-COLUMN 


NON-DIMENSIONAL  MOMENT 
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FIGURE  6. 4  MOMENT  VS  END  ROTATION 
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TGURE  6.5  MOMENT  VS  END  ROTATION 
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FIGURE  6.6  TORQUE  TWIST  RELATIONSHIP  FOR  AN  INELASTIC 
BEAM-COLUMN 
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FIGURE  6.7  LOAD  VS  MID-SPAN  DIS PLACEMENTS 


500  j- Eigen  Value  Solution  480  Inch  Kips 
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NORMALIZED  WARPING  COORDINATE 


A  =  8.08  sq.  in 

K-J-  =  0.3884  in4 

Y  =  4.624” 

lx  =  41.596  in4 

I  =  20.49  in4 

ey  =  0.754” 

I  =  66.68  in4 
P 

I  =  22.367  in6 

CO 

C  =  -  1.566  in 

y\. 

cy  =  0.0 

,0  =  2.87 
p  =  0.931 


FIGURE  6.9a  CROSS  SECTION  PROPERTIES  OF  A  DOUBLE  ANGLE  SECTION 
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FIGURE  6.9b  INELASTIC  BUCKLING  OF  DOUBLE  ANGLE  STRUT 
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FIGURE  6.10c 


RESIDUAL  STRESS  DISTRIBUTION 


(see  Fig.  6.10) 


159 


FIGURE  6.11  INELASTIC  BUCKLING  OF  A  COLUMN  (HAT  SECTION) 
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FIGURE  6.12  INELASTIC  LATERAL  BUCKLING  OF  A  WIDE  FLANGE  BEAM  (UNIFORM  MOMENT) 
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FIGURE  6.13  INELASTIC  LATERAL  BUCKLING  OF  A  BEAM  (DOUBLE  ANGLE  SECTION  -  UNIFORM  MOMENT) 
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FIGURE  6.14  INELASTIC  LATERAL  TORSIONAL  BUCKLING  OF  A  BEAM  (ANGLE  SECTION-UNIFORM  MOMENT) 
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FIGURE  6.15  INELASTIC  LATERAL  -  TORSIONAL  BUCKLING  OF  A  BEAM  OF  UNSYMMETRICAL  SECTION 

(UNIFORM  MOMENT) 


NONDIMENSIONAL  MOMENT 
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GALAMBOS  (24) 
a  FINITE  ELEMENT 


M/M 


FIGURE  6.16  INELASTIC  LATERAL  TORSIONAL  BUCKLING  OF  A  WIDE  FLANGE 
BEAM-COLUMN 
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FIGURE  6.17  INELASTIC  LATERAL  TORSIONAL  BUCKLING  OF  A  BEAM-COLUMN  WITH  MOMENT  GRADIENT 
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FIGURE  6.18  INELASTIC  LATERAL  BUCKLING  OF  WIDE-FLANGE  BEAM  (CENTRAL  LOAD) 
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FIGURE  6.19  INELASTIC  BUCKLING  OF  BIAXIAL  BENDING  WITH  AXIAL  LOAD 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 

A  general  approach  to  the  elastic  and  inelastic  analysis  of 
thin-walled  members  of  arbitrary  open  cross  section  has  been  presented 
using  finite  element  analysis. 

Both  incremental  and  total  equilibrium  equations  have  been 
derived,  using  the  principle  of  virtual  work,  which  are  applicable  to 
elastic  and  inelastic  analysis.  For  the  elastic  case,  a  special 
application  of  the  formulation  is  the  numerical  determination  of 
critical  loading,  and  the  associated  mode  shapes  which  may  be  deter¬ 
mined  from  an  eigenvalue  solution.  The  eigenvalue  analysis  essentially 
neglects  the  prebuckling  deformations.  However,  initial  imperfections 
have  been  introduced  into  the  model,  and  by  using  an  incremental 
iterative  procedure,  the  nonlinear  load-deflection  response,  rather 
than  critical  load,  has  been  obtained  and  numerical  comparisons  were 
made. 

The  displacement  functions  utilized  in  the  beam  and  beam- 
column  analyses  were  augmented  for  symmetrical  wide-flange  sections 
by  representing  local  buckling  patterns  with  the  shape  functions  shown 
in  Chapter  IV.  By  coupling  the  local  and  member  buckling,  numerical 
results  verified  the  validity  of  the  classical  approach  showing  that 
for  these  sections,  the  interactive  effect  is  small,  and  local  and 
member  buckling  are  initiated  as  relatively  detached  phenomenon  over 

most  of  the  range  of  loading. 

Two  formulations  have  been  developed  for  the  inelastic 
problems.  The  inelastic  beam-column  problems  were  solved  by  an 
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iterative  incremental  technique  based  on  an  equilibrium  balance  of 
applied  loads  and  resisting  forces  in  the  current  configuration  of 
the  structure.  These  problems  were  found  to  be  very  sensitive  to  the 
type  of  stress  strain  curve  and  applied  forces  in  the  inelastic  range 
and  hence  an  under  relaxation  factor  was  used  on  the  residual  force 
vector.  The  convergence  of  these  problems  was  slow  when  the  load  was 
near  the  peak  load.  In  spite  of  the  uae  of  a  modified  Newton- 
Raphson’s  procedure  the  solution  was  found  to  be  most  expensive  from 
the  standpoint  of  computer  expenditure. 

Inelastic  buckling  problems  were  solved  using  standard 
eigenvalue  solutions  to  evaluate  the  critical  length.  This  procedure 
is  applicable  to  thin-walled  prismatic  members,  with  arbitrary  shaped 
open  cross  section,  for  any  statically  determinate  loading  condition, 
with  a  trilinear  stress-strain  curve  and  including  the  effects  of 
residual  stress.  The  analysis  is  capable  of  predicting  lateral, 
torsional  and  lateral-torsional  buckling  modes,  and  yields  reasonable 
results  for  large  classes  of  problems.  The  formulations  which  were 
developed,  were  found  to  yield  the  same  results.  It  is  shown  in 
Appendix  B  that  for  the  class  of  problems  investigated  in  this 
dissertation  it  is  unnecessary  to  include  the  higher  order  terms  in 
the  initial  stress  and  stress-strain  increments  unless  investigating 
the  post-buckling  behaviour,  since  the  maximum  load  carrying  capacity 
will  normally  be  obtained  before  these  higher  order  terms  have  any 
significant  effect. 

Even  though  only  single  member  problems  are  solved  in  this 
work,  the  same  formulation  may  be  applied  to  solve  for  instability  of 
planar  frames  composed  of  arbitrary  open  thin-walled 
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sections.  However,  more  efficient  numerical  techniques  must  be 
developed  before  geometric  and  material  non-linearities  of  the  types 
considered  in  this  dissertation  can  be  practically  incorporated  into 
the  design  process. 
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APPENDIX  A 

DETAILS  OF  ELASTIC  FORMULATION 

A-l  Expressions  For  Torque  and  Blmcment 

The  twisting  stress  resultant,  or  torque,  about  the  arbi¬ 
trary  point  S  of  Fig.  2.1,  was  defined  in  Section  2.1  as 


A 


{o  (x-e)-a  (y-e)}dA 
zy  x  zx  y 


(2— If ) 


Let  p  be  the  distance  from  S  to  A,  and  let  r^  and  r  be  the  projec¬ 
tions  of  p  on  the  tangent  at  A  and  the  perpendicular  to  the  tangent 
at  A,  respectively,  as  shown  in  Fig.  2.1.  Then  r  and  r^  may  be 
written  as 


r 


(A-la) 


n 


=  (y  -  e  )  cos  a  -  (x  -  e  ) 
y '  X 


sm  a 


( A-lb) 


where  a  is  the  angle  between  the  y  axis  and  the  tangent  at  A. 

Substituting  the  expressions  for  components  of  shearing 
stress,  from  Equations  2—4,  into  Equation  2— If ,  and  simplifying  with 
Equation  A-la,  we  may  write 
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a  (r  -  O  dA 

zn 
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(A-2) 


where  E,  is  the  coordinate  orthogonal  to  p.  The  approximation  is 

required  since  Equations  2-4  are  not  exact. 

If  a  is  assumed  to  vary  linearly  through  the  thickness  we 
zq 


may  write 
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zn 
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zn 


a* 

-  2f  E 


(A-3) 


where  a*  is  the  deviation  of  the  stress  at  5  =  t/2  from  that  at  A, 

zn 

as  shown  in  Fig.  A.l.  Substituting  Equation  A-3  into  Equation  A— 2, 
and  treating  A  as  an  arbitrary  point,  yields 
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zn 
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where  s  is  a  coordinate  along  the  midsurface  contour  (Fig.  2.1). 
Equation  A- 4  may  be  written  symbolically  as 


M  -  M  +  T  /2 
t  0)  sv 


(A-5) 


where  M  is  known  as  the  warping  torque  and  T  is  known  as  the  St. 

03  sv 

Venant  torque.  A  study  of  torsion  of  general  cross-section  (69) 

shows  that  the  torque  resisted  by  shear  stresses  perpendicular  to 

the  contour  is  also  T  / 2  so  that  we  may  write  the  total  torque  as 

sv 


M  =  M  +  T 
t  0)  sv 


(A-6) 


If  it  is  understood  in  the  following  that  c  is  the  shear 

Z  Tj 

stress  along  the  midsurface  contour,  may  be  evaluated  as  follows. 
By  definition 
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Integrating  by  parts,  with  B  as  the  origin  of  the  coordinate  s,  yields 
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Recognizing  that  a 


zn 


must  be  zero  at  E  and  F,  and  using  the  equilibrium 
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relationship 
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Equation  A-8  may  be  written  as 
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where  we  define 
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This  may  also  be  written  as 
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where 


r  ds  =  oo 


(A-13) 


is 


the  sectorial  coordinate  of  point  A  with  respect  to  the  origin  E 


and  the  pole  S. 

The  bimoment  stress  resultant,  is  defined  as  the  quantity 
being  differentiated  in  Equation  A-10,  i.e. 

•F 
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Equation  A-6  may  now  be  written  as 
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A- 2  Derivation  of  Stress  Resultants  in  Terms  of  Displacements 

To  convert  the  equilibrium  Equations  2-2,  2-8  or  2-20  to 

displacement  equations  of  equilibrium,  it  is  necessary  to  express  the 

stress  resultants,  P ,  M  ,  M  and  M  ,  in  terms  of  displacements.  This 

x  y  t 

may  be  accomplished  by  deriving  an  expression  for  strain  ,  as  follows. 

The  in-plane  displacements  at  any  point  A  on  the  midsurface 
contour  are  illustrated  in  Fig.  2.1,  and  may  be  expressed  in  terms  of 
the  displacements  at  the  arbitrary  reference  point  S,  as 

v  =  v  +  (x  -  e  )  <J)  (2-5a) 

A  s  x' 

u  =  u  -  (y  -  e  )  <j>  (2-5b) 

A  s  y  y 

If  nA  is  the  displacement  at  A  in  the  (positive)  tangential  direction 
A 

then 

n.  =  v.  cos  a  -  u.  sin  a  (A-16) 

A  A  A 

Substituting  Equations  2-5  into  Equation  A-16,  and  simplifying  with 
Equation  A-la,  yields 

n  =  v  cos  a  -  u  sin  a  +  r 4>  (A-17) 

As  s 

Using  assumption  1,  of  Section  2.2,  we  may  now  write 

3n  +  =  o  (A-18) 

3  z  3s 

Substituting  Equation  A-17  into  Equation  A-18,  and  using  the  relations 
(see  Fig.  2.1) 
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results  in 
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Integrating,  from  B  to  an  arbitrary  point  A,  yields 
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where  lo  has  been  defined  in  Equations  A-ll  and  A-12.  Equation  A-21 
SB 

may  also  be  written  as 
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where 
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(A-23) 


The  longitudinal  strain  may  now  be  obtained  from  the  kinematic 
relationship  in  Equation  A-22  by  differentiating  with  respect  to  z. 

For  an  elastic  material  the  resulting  stress  is 


a  =  E  w '  -  E  y  v"  -  E  x  u"  +  E  co  4>"  (A  24) 

z  c  s  s  OD 

The  primary  stress  resultants,  required  in  Equations  2-20,  may  now  be 
evaluated,  for  an  elastic  response,  by  substituting  Equation  A-24  into 
the  definitions  of  Equations  2-la,  2-lb,  2-lc  and  A-14.  This  results 


in 
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Each  integral  in  the  matrix  defines  a  property  of  the  cross-section 
which  is  designated  by  the  notation  below. 


(A-26) 


The  expressions  for  the  stress  resultants  contained  in  Equation  A  26 
may  now  be  substituted  directly  into  the  equilibrium  Equations  2-20, 
although  it  leads  to  considerable  simplification  if  the  off-diagonal 
terms  are  made  to  vanish  by  properly  selecting  the  reference  points. 


Selecting  C  as  the  centroid,  and  the  axes  as  principal  axes, 
results,  by  definition,  in 


A  x  =  A  y  =  I  =  0 

3  xy 

The  sectorial  linear  static  moments,  and  S ^  vanish  if  the  point  S 

is  located  at  the  shear  center.  To  establish  the  location  of  this 
point,  the  distance  r  in  Fig.  2.1  may  be  expressed  as 
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r  =  r  -  e  cos  a  -  e  sin  a 
ox  y 


(A-27) 


Using  the  definition  A-ll  and  the  relations  A- 19  we  may  then  write 


•A 

-  r0  dS  +  ex  -  ey  (x  -  XB) 

B 


(A-28) 


Substituting  into  the  definitions  of  Equation  A~26  and  assuming  the 

centroidal  location  of  C  and  principal  orientation  of  the  axes,  the 

requirements  that  S  and  S  vanish  are 
H  oox  my 
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CB  y  xx  J 

“CE  y  ^ 


+  e  I  =  0 

x  x 


Solving  for  the  location  of  the  shear  center  yields 


(A- 29b) 


(A-30a) 


e  ,  S^-c-  (A-30b) 

x  I 

x 

Where  S  and  S  are  identified  with  the  corresponding  integrals  in 
ojxc  o)yc 

Equations  A-29 . 

The  sectorial  static  moment,  S  ,  will  vanish  if  the  point  B 
is  selected  at  a  sectorial  centroid.  To  locate  such  a  point  we  use 
the  definition  from  Equations  A~26  and  Equation  A-12  to  write 


' 


. 
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(A-31) 


Solving,  the  point  B  must  be  located  such  that 


co 


B 


co  dA 


(A-32) 


In  accordance  with  the  accepted  nomenclature,  when  S  is  the  shear  center 

and  B  is  a  sectorial  centroid,  coOT)  is  designated  as  co  ,  the  normalized 

bo  n 

unit  warping. 

Since  all  the  off-diagonal  terms  in  Equation  A-25  and  A-26 
have  now  been  made  to  vanish,  the  displacement  derivatives  are  easily 
expressed  in  terms  of  stress  resultants  and  Equation  A-24  may  be  written 
as 
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(A-33) 


For  thin-walled  elastic  sections  the  stress  resultant  may 

be  written  as  (68) 


T  =  G  K  <J>’  (A-34) 

sv 


where 


K  =  l  bt3/3 


(A-35) 


A- 3  Residual  Stress  and  Strain 

If  residual  strain  exists  in  the  section  the  total  strain  may 

be  written,  using  Equation  A-22,  as 
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e  =  w'  -  y  v"  -  x  u"  +  10  <J>"  +  e 

™  ~  °  °  bB 


(A-36) 


where  e  is  the  residual  strain.  The  stress,  for  elastic  response  is 
R 

then 


a  =  E  w’  -  E  y  v"  -  E  x  u"  +  E  to  (j>M  +  a 


(A-37) 


The  residual  stress  distribution  must  satisfy  the  following  equili¬ 
brium  requirements 


°R  ^  = 
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(A-38) 


and 
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(A-39) 


where 
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(o  dA 
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(A-40) 


and  T  is  the  residual  St.  Venant  Torque. 
svR 

A- 4  Virtual  Work  Derivation  of  Equilibrium  Equations 

The  virtual  work  equilibrium  equations  for  members  of  thin- 
walled  open  cross-sections  may  be  derived  from  Equation  2-16;  on 
substitution  of  the  relevant  form  of  displacements,  from  Equations 
2-5  and  A-22;  using  the  appropriate  assumptions  for  the  stress  field, 
expressed  by  Equations  2—19  and  2—4;  integrating  over  the  area  of  the 
section;  and  making  use  of  the  definitions  of  stress  resultants. 

Equation  2-16  has  been  written  symbolically  as  Equation  2-17 
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6  (1^  +  ~  I^)  ~  0 

where  I  ,  I  ,  and  I„  are  defined  in  Equations  2-18.  The  operation  of 
1  2  >3 

taking  the  variation  is  carried  out  separately  for  each  of  the  inte 
grals  as  follows. 

Expanding  I-  for  the  nonzero  stress  components  yields 
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(A-41) 


substituting  for  the  displacements  from  Equations  2  5  and  A  22,  yields 
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Using  Equations  2-4  this  may  be  written  as 
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+  a  (w’  -  y  v"  -  x  u"  +  m  dA  dz  (A.  43) 

z  c  s  s  SB 

The  terms  in  a  add  out  by  virtue  of  Equations  A-la  and  A-ll.  Inte- 
zq 

grating  the  remaining  a ^  terms  and  using  the  definitions  of  the  stress 
resultants,  yields 
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Expanding  for  the  nonzero  stress  components  yields 
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The  last  integral  in  this  equation  contains  product  terms  in  w  displace¬ 
ment  gradients  and  is  considered  of  higher  order. 

Substituting  the  displacement  assumptions  into  the  first  two 
integrals  of  Equation  A-45  yields 
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where  the  stress  resultant  T  is  defined  as 

rn 


rn 


a  r  dA 
zq  n 


(A-47) 


This  stress  resultant  has  not  been  included  in  any  set  of  beam  equa¬ 
tions  to  the  author’s  knowledge.  It  is  presumed  small  and  omitted  from 
the  final  equations. 

The  integral  1^  is  evaluated  by  assuming  q^,  q^.,  and  m^  act 

along  the  S  axis  and  m  ,  m  ,  and  q  act  along  the  C  axis.  Then 

x  y  z 
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Adding  Equations  A-44,  A-46  and  A-48:  dropping  the  T^  term;  and  aug¬ 
menting  the  equation  by  the  virtual  work  due  to  St.  Venant  torsion; 
Equation  2-17  becomes 
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These  are  the  equilibrium  equations  referred  to  any  arbitrary  origin 

of  coordinates  C  and  reference  axis  S  (see  Fig.  2.1)  and  make  no 

assumption  as  to  the  distribution  of  a  and  a  other  than  those 
r  z  zr\ 

required  for  St.  Venant  torsion.  They  are  therefore  suitable  to  use 
as  a  basis  for  a  finite  element  idealization. 

Integrating  by  parts  yields  the  Euler-Lagrange  equilibrium 
Equations  2-20  and  the  natural  and  geometric  boundary  conditions 
shown  in  Table  A-l.  The  force  quantities  associated  with  the  natural 
boundary  conditions  give  the  generalized  forces  corresponding  to  the 
displacement  coordinates. 


A- 5  Evaluation  of  M  for  Elastic  Response 

P  


The  stress  resultant,  M  ,  defined  as 
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arises  in  the  equilibrium  derivation  (Section  2.2)  and  the  virtual 

work  derivation  (Section  A-4)  of  the  governing  equations.  ror  elastic 

response  o  is  expressed  by  Equation  A— 37 .  If  x  and  y  are  principal 
z 

centroidal  coordinates  and  S  is  the  shear  center  Equation  A-37  may  be 
written  as  (see  Equation  A-33) 
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Substituting  Equation  A-51  into  Equation  2-7  yields 
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where 


c 

X 

-  H  /I  - 

X  X 

2  e 

y 

(A-53a) 

C 

y 

-  H  /I  - 

y  y 

2  e 

X 

(A-53b) 

c 

w 

=  H  /I 

0)  to 

(A-53c) 

H 

X 

• 

P 

y(x2  + 

i 

y2)  dA 

(A-53d) 

H 

y 

* 

P 

f  (  2 

x(x  + 

L 

y2)  dA 

(A-53e) 

H 

0) 

. 

2  2 

to  (x  +  y  )  dA 
n 

(A-53f) 

A 

I 

[  +  I 

2  2 
+  A(e  +  e  ) 

(A-53g) 

x  y 


X 


Equation  A- 52  is  valid  only  when  C  is  the  centroid ,  S  is  the  shear 
center  and  the  axes  are  principal  axes  of  the  Section.  Expressions  for 
the  numerical  evaluation  of  all  the  necessary  section  properties ,  for 
any  arbitrary  thin-walled  section,  are  given  in  Appendix  D. 

A- 6  Evaluation  of  Element  Flexural  Stiffness  Matrix  [k  ]  and 


Geometric  Stiffness  Matrix  [k  ] 


It  was  assumed  in  Section  2.7  that  continuous  displacements 
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are  given  in  terms  of  nodal  displacements  by  Equations  2-23.  The 
interpolating  functions  may  be  expressed  as 
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<fL>  =  <(1  -  3),  3> 


(A-54b) 


<f2>  =  <(23  -  1)  (3-D,  -4  3(3  -  1),  3(2  3  -  1)>  (A-54c) 

where  the  nondimensional  coordinate,  3  and  the  nodal  vectors  {u_} ,  ( v} , 

{j£_}  and  (w)  are  defined  in  Fig.  2.6. 

To  obtain  Equation  2-25  from  Equation  2-24  it  is  necessary 

to  carry  out  the  integration  in  Equation  2-24,  assuming  A,  1^,  I  ,  1^, 

P ,  M  ,  M  ,  e  ,  e  all  vary  linearly  in  the  form 
*  x*  y  x*  y 


A  =  A  +  (A  -  A  )  3  (A-55) 

4  P  9  P 


where  p  and  q  are  adjacent  nodal  numbers,  the  integrated  matrices  may 
all  be  written  in  terms  of  coefficient  matrices  of  the  form 
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o 

^  =  the  degree  of  the  interpolation  vector  (see  Equations 
and  t  indicate  the  order  of  differentiation  of  the  respective 
and  j  is  the  exponent  of  3*  The  numerical  values  of  these 
for  specific  values  of  s,  t,  m,  n  and  j  are  given  in  Table  A-2. 
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The  matrices  arising  from  Equation  2-24,  which  are 
symbolically  in  Equations  2-27  and  2-28,  then  become 
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In  Equation  A-57c 
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and  K  is  St.  Venant’s  constant. 

If  the  thin-walled  beam  is  supported  on  an  elastic  medium 
suitable  spring  constants  may  be  added  into  the  flexural  stiffness 


matrix. 
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Natural  Boundary  Conditions 

Geometric  Boundary 
Conditions 
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TABLE  A-l 


BOUNDARY  CONDITIONS 
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TABLE  A- 2 


COEFFICIENT  MATRICES 
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o-A  Warping  Shear 

z  rj 

o-*,j  St.  Venani  Shear 


FIGURE  A1  SHEAR  STRESSES  IN  A  THIN -WALLED  BEAM 


APPENDIX  B 


COMPARISON  OF  ELASTIC  FORMULATIONS 
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APPENDIX  B 


COMPARISON  OF  ELASTIC  FORMULATIONS 


B-l  Introduction 

The  governing  equations  of  this  dissertation  were  derived 
from  the  principle  of  virtual  work.  It  has  been  shown  in  Chapter  II  . 
and  Appendix  A  that  the  total  equilibrium  equations,  Equations  2-20, 
are  essentially  the  same  as  those  of  Vlassov  (74).  However,  the 
incremental  equilibrium  equations,  derived  from  virtual  work  following 
the  procedure  of  Biot  (10)  and  developed  by  Felippa  (19),  Murray  (50), 
Hofmeister  et  al  (34),  and  others,  differ,  from  the  potential  energy 
derivation,  in  the  higher  order  terms. 

The  most  consistent  derivation  of  nonlinear  terms  in  finite 
element  formulations  is  based  upon  the  potential  energy  approach  and 
due  to  Mallet  and  Marcal  (43).  The  objectives  of  this  appendix  are 

a)  to  compare  the  present  method  of  formulation  with  that  of  Mallet 
and  Marcal ; 

b)  to  indicate  the  modifications  required  in  the  virtual  work  formu¬ 
lation  to  make  it  correspond  to  the  incremental  formulation  of 
Mallet  and  Marcal; 

c)  to  demonstrate  by  numerical  comparison  that,  for  members  of  thin- 
walled  open  sections,  the  formulation  adopted  in  this  dissertation 
yields  reliable  results  in  the  range  of  interest.  The  beam-column 
formulation  of  Mallet  and  Marcal  is  used  for  this  comparison. 

B-2  Formulation  of  Beam-Column  Equations 


3-2-1  Virtual  Work  Formulation 
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The  incremental  equilibrium  equation  described  in  Chapter  II, 
is  given  by  Equation  2-14  which  may  be  written  in  the  following 
notation 


a  .  .  6 (e .  . )  dv  + 
ij  iJ 


a  .  .  6 (e .  . )  dv 
ij  ij 


t.  5u.  ds 
l  l 


(B-l) 


where 

e. .  =  e? .  +  e}.  +  e? .  (B-2) 

ij  ij  3.3  13 

0  9  1 

and  e  ,  e  ,  e  .  contain  linear  terms  in  displacement  increments  and 
13  13  i3 

product  terms  of  initial  and  incremental  displacements ,  respectively , 

i. e.  : 


0 

e.  . 
13 


+  u, 


.!> 


(B-3a) 


1 

e.  . 
ij 


u. 


k,i 


(B-3b) 


2 

e.  . 
ij 


^“k.i  uk,j +  uk,j 


(B-3c) 


In  Table  B-l  linear  incremental  formulations  are  compared  for  the  in¬ 
plane  beam-column  problem 

a)  when  only  first  order  terms  are  included  in  a  and  e_ 

b)  when  the  higher  order  terms  are  included  in  and  e_ 

B-2-2  Mallet  and  Mar cal* s  Formulation 

The  strain  energy  for  the  beam  element  is  written  as 


U  = 


E(w  -  Y  v  _  +  v  / 2)  dv 


zz 


(B-4a) 


v 
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’  <d>  [f^  [kg]  +  [nj  +  If  Cn2m d}  dz 

l 


where 


10  0 


0  0  0 


0  0  I/A_ 


0  v  0 

,z 

v  w  0 

,z  ,z 

0  0  0 


0  0  0 

0  3/2(v2  )  0 

,z 

0  0  0 


and 


<  d  > 


<  w 

>z 


V  ,  V 
,z  ,zz 


For  equilibrium 


SU  +  SW  =  0 


(B-4b) 


(B-4c) 


(B-4d) 


(B-4e) 


(B-4f ) 


(B-5a) 


Hence  equilibrium  equation  is 

<«d>  [EA[ks]  +  [nx]  +  [n2]]ld)  dz  +  6W  =  0  (B-5b) 

l 

Retaining  only  linear  terms  in  the  increments  in  displacements ,  the 
incremental  equiliorium  equation  may  be  written  as 


6 Ad  [EA  [kg]  +  EA  [n^  +  EA  [n^]^  (Ad)  dz  +  {6AW>  = 

l 


0 


(B-6) 
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Equation  B™6  is  the  same  as  the  incremental  equilibrium  equation 
developed  by  using  the  virtual  work  formulation  for  case  b  (Table  B-l) . 

Expressing  [d]  in  terms  of  the  generalized  coordinates  (3) 
through  a  linear  transformation  of  the  form 

«}  =  [D]  {6}  (B-7) 


the  potential  energy  is 

U  =  <3>  [“  [kg]  +  |  [nx]  +  [n2]l  (B)  (B-8a) 


The  equilibrium  equation  is 

[[kg]  +  |[^x]  +  ^[n2]]  (3)  +  (Q>  =  0  (B-8b) 

where  { Q }  is  the  vector  of  generalized  forces  associated  with  the  (3) 
displacement  coordinates. 

The  incremental  equilibrium  equation  is 


[[kg]  +  [n1]  +  [n2]]^  +  (aQ}  =  0 


where 


rl 


[ts] 


EA  [D]  [kg]  [D]  dz 


rl 


[ni] 


EA  [D]1  [nx]  [D]  dz 


o 


[n2] 


l 


EA  [D]  [n2]  [D]  dz 


(B-8c) 


(B-9a) 


(B-9b) 


(B-9c) 


B-2-3  Comparison  of  Formulations 


The  derivations  in  Chapter  II  and  Appendix  A  show  that  by 
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ignoring  the  initial  displacement  matrix,  and  using  only  first  order 

terms  in  a...  a.,  and  e.#  (as  in  Table  B-l,  case  a)  the  virtual  work 
ij  ij  iJ 

formulation  is  consistent  with  the  formulation  of  the  classical  beam 
equations  normally  applied  to  stability  problems.  These  equations 
may  be  regarded  as  examples  of  the  lowest  level  of  nonlinear  equations 
in  the  heirarchy  of  nonlinear  formulations  discussed  by  Mallet  and 
Marcal.  In  this  work  the  term  "virtual  work  equations"  refers  to  sets 
of  equations  derived  by  assumptions  consistent  with  those  for  case  a, 
of  Table  B-l.  The  term  "potential  energy  equations"  refers  to  sets 
of  equations  derived  in  a  manner  consistent  with  the  derivation  of 
Mallet  and  Marcal,  or  case  b.  Table  B-l. 

B-3  Nonlinear  Analysis  of  Beam-Columns 

In  Section  B-2  a  comparison  was  made  of  the  nonlinear  terms 
arising  in  the  virtual  work  and  potential  energy  formulations.  The 
question  arises  as  to  the  range  of  reliability  of  the  virtual  work 
formulation.  Some  indication  of  this  is  achieved  by  comparing  the 
numerical  result  for  a  typical  in-plane  beam-column  problem.  Such  a 
comparison  is  carried  out  in  the  following. 

B - 3 - 1  Mallet  and  Marcal' s  Linear  Incremental  Formulation 

Using  the  nomenclature  of  Mallet  and  Marcal,  the  linear 
incremental  potential  energy  formulation  is  written  as 

[ [k  ]  +  [N  ]  +  [n  ]]-  {Ar}  =  {AR}  (B-10) 

b  -J*  Z.  1C 

where  element  matrices  are  given  by 


[ks]  -  [t]T  [ks]  [T] 


(B-lla) 
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[nx]  -• 

[T]T 

[nj 

[T] 

(B-llb) 

[n2]  - 

[t]t 

[n2] 

[T] 

(B-llc) 

and  [k  ]  ,  [n  ]  ,  [n  ]  and  [t]  are  specified  in  Figs.  B.l,  B.2,  B.3,  and 

u  -L 

B.4  respectively.  The  vector  of  generalized  coordinates,  (3)>  specifies 
the  centroidal  displacements  through  the  relations 


w  =  3 +  3 2  z  (B-12a) 

and 

v  =  +  3^  z  +  5^  z2  +  3^  z3  (B-12b) 

The  generalized  coordinates  are  related  to  the  nodal  displacements,  {r} 
through  the  transformation 

(3)  =  [T]  {r}  (B-13) 


where 


VP’ 


w  > 

q 


(B-14) 


The  total  equilibrium  equation  may  be  written  as 

[[Ks]  +  |  [Nx]  +  i  [N2]]  {r}  =  {R>  (B-15) 

The  solution  to  this  set  of  equations  was  obtained  by  a  Newton-Raphson 
procedure  where  iteration  was  carried  out  until  the  unbalanced  nodal 
forces,  { E } ,  evaluated  as 

{E}  =  {R}  -  [[Ks]+  \  [tq]  +  i  [N2]]  {r} 


(B-16) 


■ 
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were  negligible.  A  flow  chart  of  this  procedure  is  shown  in  Fig .  B.5. 

B-3-2  Incremental  Virtual  Work  Formulation 

The  incremental  virtual  work  equations  may  be  written  as 

[ [K  ]  +  [K  1]  {Ar}  =  {AR}  (B-17) 

b  ^ 

The  element  geometric  stiffness  matrix,  [kG]  for  the  generalized  coor¬ 
dinates  is  obtained  by  summing  [n^]  and  [n^]  of  the  previous  formulation, 

with  8.,  3  ,  and  6,  equal  to  zero  and  =  P./AE.  This  is  equivalent  to 
4  5b  ^ 

discarding  the  matrix  [n2]  and  retaining  only  the  3x3  main  diagonal 
submatrix  in  the  definition  of  [n^] . 

The  total  equilibrium  equation  is  written  as 

[[Ks]  +  [KG]]  {r}  =  {R}  (B-18) 

The  solution  to  this  set  of  equations  was  also  obtained  by  a  Newton- 
Raphson  procedure  where  iteration  was  carried  out  until  the  unbalanced 
nodal  forces,  {E}  evaluated  as 

{E>  =  {R>  -  [[Kg]  +  [Kg]]  {r}  (B-19) 


were  negligible. 

B-4  Numerical  Results 

The  beam-column  shown  in  Fig.  B.6  was  subjected  to  an 

increasing  axial  load  and  a  constant  transverse  load  of  5  kips  applied 

at  the  central  point,  A.  The  contribution  of  the  matrices  [Kg]  and 

[K  ]  to  the  equilibrium  of  point  A  is  shown  in  Fig.  B.6.  The  stiff- 
G 

ness  of  the  structure  approaches  zero  as  P/Pcr  approaches  1. 


■ 


. 


* 
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The  same  beam-column  was  analyzed  by  the  method  of  Mallet 
and  Mar cal.  The  contribution  of  the  matrices  to  the  equilibrium  of 
the  transverse  force  at  A  are  shown  in  Table  B-2.  The  results  are 
also  plotted  in  Fig.  B.7. 

A  comparison  of  the  load  deflection  plots  by  the  two  methods 
is  shown  in  Fig.  B.8.  While  it  is  apparent  that  the  second  order 
incremental  stiffness  matrix  [N^]  or  a  suitable  equivalent  must  be 
included  to  analyze  post-buckling  behaviour,  it  is  also  apparent  that 
the  pre— buckling  behaviour  may  be  closely  approximated  by  omitting 
these  terms.  In  addition,  inelastic  material  response  may  be  expected 
to  occur  prior  to  any  significant  stiffening  effect  from  the  difference 
in  the  large  displacement  terms  between  the  matrices  [N^  and  [N2]. 
qt  3L3  therefore  concluded  that  the  inclusion  of  inelastic  material 
response  is  more  significant  for  this  class  of  problem  than  the 
inclusion  of  higher  order  geometric  effects. 

B-5  Summary  and  Conclusion 

It  has  been  shown  that  the  virtual  work  formulation  can  be 
made  to  correspond  to  the  linear  incremental  potential  energy  formula 
tion  by  including  the  higher  order  terms  in  the  expression  for  initial 

stress,  and  stress  and  strain  increments. 

It  is  seen  from  numerical  results  that  it  is  unnecessary  to 

include  the  higher  order  terms  unless  investigating  the  post-buckling 
behaviour.  Because  of  the  influence  of  inelastic  material  response 
the  maximum  load  carrying  capacity  will  normally  be  attained  before 
these  higher  order  terms  have  any  significant  effect. 

It  has  also  been  shown  that  the  virtual  x^ork  formulation 
employed  is  consistent  with  the  classical  differential  equation 


N 
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formulation  of  the  problem. 


TABLE  B-l  LINEAR  INCREMENTAL  FORMULATION 


211 


Td 

3 


0) 


to 


•r-> 

to 

‘H 

•H 

CM 

0k 

0) 

> 

to 

a 

•H 

to 

CM 

to 

0k 

CO 

to 

1  > 

a 

CM 

G 

to 

CM 

+ 

cu 

to 

0k 

+J 

to 

> 

CM 

a 

N 

to 

to 

>s 

Csl 

0k 

to 

o 

CM 

•k 

> 

Td 

to 

|  \ 

> 

to 

rH 

o 

to¬ 

CM  CM 

CM 

co 

»>  CM  •* 

0t 

1 

Td 

•H 

W 

V 

1  > 

to 

Td 

Csl 

o 

0k 

o 

i — i 

£ 

to 

to 

c/o 

•H 

to 

•rH 

^ — V 

to 

-to 

•rH 

T1 

•H 


0) 

to 

•H 


•H 

to 

to 

to 

to 

0k 

e 

•r-) 

to 

o 

o 

to 

rH 

to 

p- 

*H 

CO 

•H 

CO 

'to 

6 

to 

rH 

to 

to 

to 

•H 

to 

to 

•H 

to 

to 

id 

•H 

to 

o 

bO 

to 

to 

•H 

CO 

to 

to 

O 

*H 

to 

fH 

i — 1 

bO 

to 

s — s 

to 

to 

CM 

CM 


> 

I 


CM 


CM 


CM 

CM  * 

> 


'to 

to 

to 

o 

to 

rH 

bO 

to 

to 


O  *rt 
to 


H  *H 
to 


Si 

CM 

#» 

> 

>s 


CM 


5 

II 

to 


CM 


CM  -H 
to 

+ 


to 

+ 

•r~) 
O  -H 
to 


•o 

CM  -H 
to 


II 


to 


o 

o 

o 

tod 

to 

Csl 

/-—S 

to 

0k 

CM 

0k 

0k 

to 

> 

to 

CM  CM 

rH 

CM 

CM  M 

i  i  ^ 

b£> 

0k 

0k 

1  >  1 

to 

I  > 

1  > 

o 

H* 

O 

to 

CM 

+ 

+ 

«** 

CO 

1  £ 

to 

CM 

CM 

to 

Csl 

CM 

to 

0k 

0k 

0k 

0k 

0k 

6 

> 

1  > 

O 

O 

o 

to 

1 

1 

to 

>s 

to 

H, - 

- 

- N 

rH 

| 

1 

f 

CM 

p- 

CM 

CM 

CM 

CO 

Csl 

CM 

0k 

0k 

0k 

•H 

0k 

0k 

> 

> 

> 

'to 

1  ^ 

°o 

<0 

°o 

l — 1 

w 

w 

V  , 

- V — - - 

_ J 

to 

' - 

to 

to 

<i 

to 

w 

S 

to 

II 

to 

o 

to 

•H 

CM 

'to 

in 

t 

- ^ 

- N 

CM 

CO 

CM 

CM 

CM 

g 

0k 

0k 

0k 

to 

S 

> 

> 

to 

J 

to 

V 

- v" 

to 

r 

- 

o 

c 

o 

O 

to 

■to 

/-N 

o 

Csl 

CM 

CM 

to 

CM 

CM 

o 

0k 

o 

Cto 

0k 

0k 

1  S 

> 

1  > 

>» 

1° 

O 

Ti 

N 

CM 

0k 

0k 

H 

js 

13 

/* - 

- \ 

tS3 

w 

w 

CM 

CM 

CM 

0k 

0k 

0k 

n 

II 

Es 

> 

> 

<o 

<0 

csl 

CM 

V. 

- v - 

_ J 

to 

1  to 

* - 

< 

w 


II 


<o 


•r-)  'I-) 

•H  *H 

O  D 


O 


> 


CM 

'to 


CM  CM 


£  >  > 


zz 


212 


213 


214 


Csl 

I 

pq 

w 

pq 

<1 

H 


1900 


215 


0  0  0  0  0  0 


0 


0  0  0  0 


0 


0  0 

0  0 


0  0 

0  0 


0 

0 


0  0 

0  0 


0  0 

0  0 


41 l  61 l2 

A  A 

6ll2  12ll3 

A  A 


Fig.  B . 1  FLEXURAL  STIFFNESS  (GENERALIZED  COORDINATES) 
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Fig.  33.4  TRANSFORMATION  MATRIX  RELATING  GENERALIZED 


COORDINATES  TO  ELEMENT  DISPLACEMENTS 
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Fig.  B.5  FLOW  CHART  EXPLAINING  THE  SOLUTION  SCHEME 


FOR  BEAM-COLUMN  PROBLEMS 
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FIGURE  B . 6  EFFECT 


OF  INCREMENTAL  MATRICES  IN  SMALL  DISPLACEMENT  RANGE 
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B.7  effect  of  incremental  matrices  in  large  displacement  range 
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FIGURE  B . 8  PRE  AND  POST  BUCKLING  BEHAVIOUR  OF  BEAM-COLUMN 
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APPENDIX  C 


DERIVATION  OF  MATRICES  FOR  LOCAL  BUCKLING 


The  element  stiffness  matrices  appearing  in  Equation  4-22  are 
obtained  by  substituting  Equation  2-23,  4-17  and  4-21  into  4-14;  integrating 
carrying  out  the  variation;  and  grouping  terms  associated  with  each  set 
of  virtual  nodal  displacements.  Assuming  the  moment  is  constant  in  an 
element  length  for  the  evaluation  of  the  geometric  stiffness  of  a  plate 
segment,  this  results  in  the  following  set  of  equilibrium  equations. 
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The  matrices  [k^],  [k^]  ,  [k^qL  [k^]  are  the  saine  as  [K^°l> 

330  330  ^31 

[K22  ] ,  [K10  ] ,  [K21  ]  derived  in  Appendix  A  in  Table  A-2.  The  matrices 
[Kt]L  [Kt2L  ,  [K^],  [Gt],  [Gb],  [a],  [b],  [c],  [d],  [f],  [h] 


are  shown  in  Tables  C-l,  C-2,  C-3,  and  C-4.  The  notation  A  designates 


the  area  of  the  web  and  I  is  the  moment  of  inertia  of  the  web. 
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TABLE  C-l 
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-3 
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TABLE  02 
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32 
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88 

16 
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216 
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APPENDIX  D 

CALCULATION  OF  SECTION  PROPERTIES  AND  STRESS  RESULTANTS 

FOR  THIN-WALLED  SECTIONS 

D-l  Introduction 

In  the  evaluation  of  the  stiffness  matrices  for  elastic 
sections,  carried  out  in  Appendix  A,  Section  A-6 ,  it  is  necessary 
to  know  the  section  properties.  In  the  evaluation  of  the  stiffness 
matrices  for  inelastic  analysis,  carried  out  in  Appendix  E,  it  is 
necessary  to  know  the  section  properties  of  the  transformed  section. 
Once  the  transformed  section  has  been  determined,  the  evaluation  of 
section  properties  is  identical  to  that  of  any  arbitrary  elastic 
section.  These  section  properties  are  always  referred  to  the 
centroid,  shear  center  and  principal  axes  of  the  effective  elastic 
(transformed)  section. 

In  this  appendix  the  method  of  determining  the  transformed 
section  is  discussed  in  Section  D-2,  and  the  equations  for  numerically 
evaluating  section  properties  are  specified  in  Section  D-3.  Equations 
for  numerically  evaluating  stress  resultants,  which  are  required  in 
the  calculation  of  [K  ]  and  in  the  evaluation  of  the  unbalanced  force 
vector  for  inelastic  analysis,  are  given  in  Section  D-4. 

D-2  Determination  of  Transformed  Section 

The  strain  at  any  point  of  a  thin-walled  section  is  due  to 
axial,  bending  and  warping  displacements,  together  with  any  residual 
strain,  and  for  any  set  of  nodal  displacements  may  be  evaluated  from 
Equation  A-36.  If  the  residual  strain  distribution  varies  in  some 
arbitrary  manner,  as  shown  in  Fig.  D.lb,  then  it  can  be  approximated 


. 


by  linear  segments  as  shown  in  Fig.  D.lc.  For  equilibrium 
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a.  dA  =  0 

K 

A 

* 

a  y  dA  =  0 

K 

A 

a  x  dA  =  0 
K 

A 


(D-la) 

(D-lb) 

(D-lc) 


Considering  any  plate  segment  which  has  a  linear  variation  of 
residual  strain,  by  superposition,  the  total  strain  distribution  will 
also  be  linear  in  the  segment.  The  stress  strain  curve  used  is  shown 
in  Fig.  D.la.  It  is  tri-linear,  where  the  slopes  of  the  three  seg¬ 
ments  are  designated  as  E  =  elastic  modulus;  E^  =  plastic  modulus; 

and  E  =  strain  hardening  modulus.  To  approximate  the  idealized 
s 

stress  strain  curve  for  mild  steel,  E^  may  be  set  equal  to  zero. 

Considering  a  segment  under  a  linear  variation  of  strain,  as 
shown  in  Fig.  D.2,  when  £^  and  e^  have  the  same  sign,  and  e^  may  be 
in  any  of  the  three  strain  ranges. 
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(D-2a) 

(D-2b) 

(D-2c) 


Hence  there  are  nine  combinations  of  strain  distribution  when  e^  and 

e  are  of  the  same  sign  in  a  segment,  and  there  will  be  nine  other 
B 

combinations  of  the  strain  distribution  when  e^  and  £g  are  of  opposite 
signs.  In  each  segment  the  incremental  moduli  are  determined  from  the 
strain  distribution,  assuming  no  strain  reversal.  This  divides  the 
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plate  segment  into,  at  the  most,  five  regions  as  shown  in  Fig.  D.2. 

It  is  now  possible  to  transform  the  section  and  evaluate  section 
properties  of  the  transformed  area. 

D-3  Evaluation  of  Section  Properties 

To  determine  the  section  properties,  consider  the  segment  of 
a  plate  of  a  thin-walled  section  under  linear  strain  distribution  as 
shown  in  Fig.  D.2.  As  stated  above,  for  any  arbitrary  segment,  there 
are  at  most  five  possible  regions  in  which  modulii  are  different.  The 
moment  of  inertia  of  the  thin-wall  about  its  y-axis  is  very  small  and 
hence  in  order  to  formulate  tangent  properties  it  is  sufficiently 
accurate  to  modify  the  thickness  between  the  corresponding  regions  by 
the  modular  ratio.  Areas  and  moments  of  inertia  of  each  region  of  the 
plate  segment  are  calculated  about  the  centroidal  axis  of  the  trans¬ 
formed  segment.  Moments  of  inertia  of  each  plate  segment  are  then 
transformed  to  the  global  axis  orientation.  A  similar  computation 
may  now  be  performed  to  determine  the  centroid  and  moments  of  inertia 
of  the  total  cross  section  by  summing  the  properties  for  all  the 
segments.  Let  C  be  the  centroid  of  the  transformed  area  and  0  an 
arbitrary  origin  of  coordinates  X  and  Y.  Referring  to  Fig.  D . 3b ,  we 
may  write  for  the  centroid 
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(D-3a) 
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n 


(D-3b) 


where  A^,  X^,  and  are  the  transformed  area  and  distances  to  the 
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til 

transformed  centroid  of  the  k  plate  segment,  and  there  are  n  plate 
segments.  The  section  properties  with  respect  to  the  x,  y  axes  of 
Fig.  D.3b  are 
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tan  2  a  =  2  I-,-  /  (I-  -  I-) 

xy  y  x 


(D-4d) 


2  2 

I_  =  I-  cos  a  +  I-  sin  a  -  sin  2  a 
5  x  y  xy 


(D-4e) 


2  2 

I  =1-  cos  a  +  I-  sin  a  4-  sin  2  a 
ri  y  x  xy 


(D-4f) 


where  a  is  the  inclination  of  the  principal  axes,  (£  and  n) ,  with 
reference  to  the  original  axes  and  is  positive  if  it  is  anticlockwise 

A 

from  the  x  axis. 

In  order  to  locate  the  shear  center,  choose  C  as  the  origin 
and  CE  as  the  initial  radius.  For  straight  segments,  sectorial  area 
coordinates  at  the  beginning  and  the  end  points  of  each  region  of  the 
plate  may  be  computed,  from  the  definition  of  the  sectorial  coordinate 
in  Equation  A-13,  as  twice  the  sectorial  area  with  respect  to  the 
initial  radius,  CE.  This  is  evaluated,  see  Fig.  D.3b,  as 


/V  A  A  A 


co.  =  co.  +  x.  y.  -  y.  x. 
J  i  J  i  J  * 


(D-5) 
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The  sectorial  static  moment  and  the  sectorial  linear  statical  moment 


about  the  x  and  y  axes,  with  C  as  the  origin  and  CE  as  the  initial 
radius,  are  defined  (see  Equations  A-29 ,  A-30  and  A-31)  as 
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S  = 
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co  y  dA 
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S  = 
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(D-6c) 


The  coordinates  of  the  shear  center,  S,  of  the  transformed  section  may 
then  be  determined  by  the  following  equations,  which  are  Equations 
A~30  for  non-principal  axes  directions. 


xt  —  (S  ^  _  X  ~  S  ~r,  I  ~  ^  )  /  (1-^  X  ^  I  ~  ~  ) 

S  ooxC  y  ooyC  xy  x  y  xy 
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—  (S  I *  —  S  I~~)  /  (I  ~  I ~  -  I  ~  ) 
ooyC  x  coxC  xy  x  y  xy 


(D-7b) 


Sectorial  area  coordinates,  co,  are  now  calculated  with  respect  to  the 
shear  center  S  and  with  SE  as  the  initial  radius,  by  a  similar  com¬ 
putation  to  that  in  Equation  D-5. 

A  principal  radius,  S  B’ ,  is  defined  as  one  for  which  the 
sectorial  coordinates  satisfy  the  following  equation  (see  Equation  A-31) 
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where 


W  SB' 
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r  ds 
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(D-8b) 


becomes  the  principal  sectorial  coordinate  of  point  A.  Choosing  to—  i 

SB 

as  the  principal  sectorial  coordinate  with  S  as  the  origin  and  SB'  as 

the  principal  radius  the  principal  sectorial  static  moment,  S  ,  will 

to 

therefore  be  zero.  The  quantity  to^-  ,  is  usually  called  the  normalized 
unit  warping  or  the  normalized  warping  coordinate,  and  designated  as  co^. 
Now 
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Hence  the  normalized  warping  coordinate  may  be  computed  by  subtracting 
the  constant  C  from  the  to  values  and  changing  sign. 

Equation  D-8a  requires  C  to  be 
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The  principal  sectorial  moment  or  simply  sectorial  moment  of  Inertia 

I  is  now  given  by 
con  -c 
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which  for  the  type  of  section  shown  in  Fig.  D.3b  can  be  expressed  as 
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v  — O  O 

I  =  I  I  b  t  (  w  +  w  w  +  to  )  /3  (D-llb) 

11  k-^l  r=l  L  L  ^  J  J 

where  cu_^  represents  the  normalized  warping  coordinate  at  junction  i; 

r  indicates  the  region  of  the  segment;  and  n  is  the  number  of  segments. 

If  it  is  assumed  that  G  is  proportional  to  E,  St.  Venant’s 

constant  may  be  calculated  by  modifying  the  width.  It  is  therefore 
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where 

b*  =  E  b  /E  (D-12b) 

r  t  r 

Other  properties  which  are  necessary  for  the  elastic  analysis  of  the 
thin-walled  structures,  are  those  appearing  in  Equation  2-24  and  are 
derived  for  elastic  sections  in  Appendix  A,  Section  A-2. 

In  these  equations 
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where  x  and  y_  are  the  coordinates  of  the  shear  centre  with  reference 

¥  ¥ 

to  centroidal  axes.  Fig.  D.3;  and. 
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By  referring  to  Fig.  D.4  the  constant  Ih,  H-  and  may  be 
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evaluated  as 


H' 


n 


yv/^2  ~2.  ,  r* 

x(x  +  y  )  dA  =  £ 


k=l 


x  .Lt 
12 


(2  pf  - 


2  2 
LZ  +  4  ($Z) 


XiLt  2  2  2 

+  i2~  ^  P-  -  L  +  2  pZ) 


(D-15a) 


H, 


n 


y  (x2  +  y2)  dA  =  l 


k=l 


YiLt  .2 
12  (2  pi 


2  2 
LZ  +  4  pZ) 


y.Lt  r.  -  _ 

+  i2~  (4  pi  -  L  +  2  Pj) 


(D-15b) 


H 


w 


n 


a)  (x  +  y  )  dA  =  £ 

k=l 


w  .Lt  ~ 

1 _  ,0  ^2 

12  (2  Pi 


2  2 
LZ  +  4  pZ) 


w.Lt  ?  ?  9  - 

+  -ir  (4  4  - L  +  2  V 


(D-15c) 


For  elastic  sections  the  coordinates  x,  y  are  identical  to  x  and  y 
and  C,  S  are  the  same  as  C  and  S  (Fig.  D.3a).  The  reference  axes 
x  and  y  through  C  of  elastic  section,  the  shear  center  S  of  the  elastic 
section,  the  principal  axes  £  and  n  through  C  of  transformed  section, 
and  the  shear  center  S  of  the  transformed  section,  are  shown  in 
Fig.  D.3c. 


D-4  Evaluation  of  Stress  Resultants 

Assuming  the  strain  is  known,  the  stresses  are  calculated  in 
each  region  of  the  plate  segments.  The  stress  resultants  may  then  be 
evaluated  numerically,  in  the  C,  S  and  x-y  reference  system,  from  their 
definition,  as  follows  (see  Fig.  D.5) 
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where  r  =  plate  region;  n  =  number  of  plate  segments;  k  =  plate  segment 
index;  t  =  plate  segment  thickness;  b  =  plate  region  length;  and 
m  =  the  normalized  warping  coordinate  with  respect  to  S. 

The  quantity  M  may  be  evaluated  about  S  as 
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where  p  and  p  are  the  distances  of  the  l  and  j  nodal  points  of 
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STRAIN,  £ 


FIGURE  D. la  TRI-LINE AR  STRESS-STRAIN  DIAGRAM 


FIGURE  l.b  ACTUAL  RESIDUAL  STRAIN 

DISTRIBUTION 


FIGURE  D.lc  ASSUMED  RESIDUAL  STRAIN 
DISTRIBUTION 
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FIGURE  D . 2  TRANSFORMED  SECTION  OF  A  PLATE  SEGMENT 
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c  Centroid  of  Elastic  Section 
s  Shear  Centre  of  Elastic  Section 
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FIGURE  D . 5  STRESSES  IN  A  TYPICAL  PLATE  SEGMENT 
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APPENDIX  E 

DETAILS  OF  INELASTIC  FORMULATION 

E-l  Introduction 

In  this  appendix  the  tangent  and  geometric  stiffness  matrices 
are  derived,  in  Section  E-2,  for  'Formulation  1'  of  Section  5.2.5,  and 
in  Section  E-3,  for  'Formulation  2'  of  Section  5.2.6.  The  section 
properties  and  stress  resultants  for  the  evaluation  of  tangent  and 
geometric  stiffness  matrices  have  been  determined  as  detailed  in 
Appendix  D. 


E-2  Derivation  of  Tangent  Stiffness  and  Geometric  Stiffness  Matrices 

For  Formulation  1  (Section  5.2.5) 

Since  variations  in  incremental  displacements  are  taken  with 
respect  to  the  reference  points  in  the  original  undeformed  configuration 
the  tangent  stiffness  matrix  contains  coupled  terms.  The  terms  con¬ 
tributing  to  the  tangent  stiffness  in  the  virtual  work  equation  Equation 


5-5  are  given  by 
rl 
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(E-2b) 

(E-2c) 

(E-2d) 


(E-2e) 


y2  dA 


(E-2f ) 


x2  dA 


x  y  dA 


(E-2g) 

(E-2h) 


and  all  the  integrals  are  computed  over  the  transformed  area.  The 

A 

quantity  oo  is  the  normalized  warping  coordinate  of  the  original 
elastic  section  with  pole  at  S. 

The  same  coupled  terms  in  the  tangent  stiffness  may  be 
arrived  at  by  using  ’Formulation  2',  Equations  5-7  and  5-8,  and  the 
transformation  equations,  Equations  5-17  and  5-22,  where  u_ ,  v^  and 
w  are  expressed  as 


u  =  u  C  +  v  S  +  { (b  -  e)S  -(b  ~e)C  }cj)(5-17a) 
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Then  we  can  identify  the  following  relationships 
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(E-3a) 
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(E-3f ) 


where 


x  =  (b  -  e  )C  +(b  -  e  )S 

x  x  a  y  y  a 


(E-4a) 


y  =  (b  -  e)C  -  (b  -  e)S 
y  y  a  x  xa 


(E-4b) 


and  C  =  cos  a ;  S  =  sin  a 
a  a 

By  either  of  the  above  procedures,  the  increraental  tangent 
stiffness  matrix  for  an  element  may  be  written  as 
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where 
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The  coefficient  matrices  in  Equations  5 --6  have  been  defined  in 
Appendix  A,  Section  A-6. 

The  matrix  [k  ]  is  the  same  as  the  geometric  stiffness 
matrix  of  an  element  for  the  elastic  case,  i.e.  Equation  2-28,  except 
for  the  submatrix  [g^],  which  may  now  be  expressed  as 

f  -  %  [k”°]  +  (I  -  I)  [k^1]  (E-7a) 

where 

K  =  K  +  (M  )p/G  (E-7b) 

P  P  ps 

Hence  the  incremental  equilibrium  equation  for  an  element  may  be 
written  as 

[Od]  +  [kG]]  {ArE)  =  {ARe)  (E-8) 

which  is  Equation  5-9. 

After  transformation  and  assembly.  Equation  E-8  becomes 

[[Kg]  +  [kg1]  (Ar)  =  {AR}  (E-9a) 

which  is  Equation  5-31.  Grouping  terms  from  these  matrices  we  may 


write 
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[[Kx]  +  [K2]  +  [K3]]  {Ar}  +  [Kg]  {Ar}  =  {AR}  (E-9b) 

where  [K^]  ,  C^],  and  [K3]  contain  terms ,  dependent  on  the  inverse  of 

2  3 

Z,  Z  and  Z  respectively  and  [K  ]  is  linearly  dependent  on  the 
inverse  of  Z. 


E-3  Derivation  of  Tangent  and  Geometric  Stiffness  Matrices  for 
’Formulation  2*  (Section  5.2.6) 

When  the  incremental  equilibrium  equation,  Equation  5-4,  is 
expressed  with  respect  to  the  instantaneous  shear  centre  and  sectorial 
centroid  of  the  transformed  area.  Equation  5-7  will  apply.  Because 
of  the  form  of  Equations  5-8,  the  incremental  tangent  stiffness  for 
an  element  may  be  written  as 
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(E-llc) 


I  ft1  ]  =  AT  [k??°]  +  (AT  -  AT)  [k??1] 
E  ww  p  11  q  p  11 


(E-lld) 


The  matrix  [k  ]  is  the  same  as  in  Equation  2-28  provided  M  , 

X 

M  are  replaced  by  Mr,  M  and  e  ,  e  by  er  and  e  .  However,  the 
y  3  5  n  x*  y  J  z  n 

sub matrix  [g  ]  is  now  given  by 
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(E-12a) 


where 


K  =  K  +  (M  — )  /G  (E-12b) 

p  P  ps  p 


The  relationship  between  M  in  Equation  E-7b  and  M  —  in 

r  p  s  p  s 

Equation  E-12b  is  given  by 


M  =M-+2M(b  -  e)+2M(b 

ps  ps  y  x  x  x  y 
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+  (b  -  e  )2  -  2(b  -  e  )  (x  +  x-)  -  2(b  -  e  )  (y  +  y~) }  (E-13) 

y  y  x  x  o  s  y  y  s 

where  x  ,  y  are  the  coordinates  of  the  instantaneous  centroid  with 
o  o 

/V  A 

respect  to  the  original  centroidal  reference  axes,  and  x— ,  y-  are  the 

s  s 

coordinates  of  the  instantaneous  shear  centre  with  respect  to  the 
instantaneous  centroidal  axes  parallel  to  the  original  reference 
axes  (see  Fig.  5.3  and  Fig.  D.3b). 

Hence  the  incremental  equilibrium  equation  for  an  element  is 


again  written  as 


Using  the  transformation  matrix,  derived  in  Section  5.4.3  for 
’Formulation  2’,  which  is  shown  in  Table  5-1,  Equation  E-14a  becomes 


[[t]T  [A]  [T]  +  [t]T  [kG]  [T]]  {ArE}Q  =  {ARe}g 


After  assembly.  Equation  E-14b  becomes  Equation  5-31,  namely 
[[KT]  +  [K  ]]  {Ar}  =  {AR} 

S  U 


(E-14b) 


(E-15) 


which  is  the  incremental  equilibrium  equation  for  beam-column  problems. 


